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Abstract 

An  increased  need  for  higher  resolution  space-based  optics  creates  requirements  for  enlarged 
apertures  on  the  primary  optical  surface.  Traditional  space-based  optical  systems  are  limited  in 
aperture  size  by  weight  and  the  payload  capacity  of  current  launch  vehicles.  Membrane-like  optical 
structures  which  could  be  rolled  for  space-lift  and  deployed  in  orbit  have  the  potential  to  greatly 
reduce  weight  and  increase  the  aperture  size  of  next  generation  mirrors  in  space.  However,  thin, 
deformable  membrane-like  structures  require  active  control  for  stabilization  and  shaping  of  the  optical 
surface.  The  primary  goal  of  this  research  is  to  obtain  a  methodology  for  the  characterization  of 
piezo-actuated,  thin,  deformable  optical  structures.  The  response  of  this  membrane-like  structure  to 
dynamic  disturbances  is  given  particular  focus.  Modal  analysis  and  active  wavefront  measurements 
were  carried  out  for  a  five  inch  membrane-like  optical  structure  and  the  results  are  compared  to 
theoretical  analysis.  Control  patches  etched  into  the  piezoelectric  polymer  material  are  used  to 
characterize  the  optical  surface  response  subject  to  static  and  sinusoidal  forcing  functions.  Various 
fabrication  techniques  are  also  investigated  for  a  polyvinylidene  fluoride  (PVDF)  membrane-like 
structure  with  a  silicone  substrate  surface  coating. 
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Dynamic  Characterization  of  Thin  Deformable  PVDF  Mirror 


I.  Introduction 

1 . 1  Overview 

Satellite  imaging  systems  are  critical  to  the  information  gathering  needs  of  the  U.S.  Department 
of  Defense.  The  aperture  size  of  the  primary  optical  surface  determines  the  detail  level  and  distance 
at  which  the  system  is  able  to  gather  useful  data.  An  increase  in  aperture  size  translates  to  an 
increase  in  observed  detail,  referred  to  as  resolution.  An  increase  in  aperture  size  also  allows  the 
optical  system  to  take  quality  images  from  a  larger  distance,  which  equates  to  a  larger  coverage  area 
on  the  surface  of  the  earth,  reducing  the  need  for  large  satellite  constellations  in  low  earth  orbit. 

An  estimate  of  the  resolution  of  an  optical  system  is  given  by  Equation  1.  This  is  known  as 
Rayleigh’s  resolution  criterion  [12]. 


Or  =  1-22^  (1) 

The  variables  A  and  d  are  the  wavelength  and  aperture  diameter  of  the  primary  mirror,  re¬ 
spectively.  Or  represents  the  minimum  angular  separation  that  two  point  sources  can  be  resolved  as 
separate  entities.  Notice  that  d  is  inversely  proportional  to  Or,  which  means  increasing  the  diameter 
of  a  mirror  will  enable  a  smaller  angular  separation  between  two  given  objects  to  be  resolved. 

Figure  1  shows  a  typical  scenario  for  satellite  observation.  From  these  dimensions,  Equation  2 
can  be  used  to  calculate  the  minimum  distance  that  two  objects  on  earth  can  be  discerned  as  separate 
entities. 


L  =  2h  *  tan 


(2) 


Let  L  be  the  minimum  resolvable  separation  distance  from  an  assumed  geosynchronous  orbit 
of  h=  35,  800  km.  Since  the  minimum  angular  separation  increases  as  wavelength  is  increased,  an 
approximation  for  the  largest  wavelength  in  the  visible  spectrum  is  used,  A  =  700  nm.  This  will 
give  a  ‘worst  case’  estimate  for  the  entire  visible  spectrum.  Table  1  provides  a  summary  of  these 
calculations. 

These  calculations  clearly  illustrate  the  fact  that  the  larger  the  lens,  the  finer  the  details  that 
can  be  collected.  There  are  numerous  reasons  that  space-based  optical  systems  created  to  date  have 
not  incorporated  the  large  aperture  sizes  needed  for  a  very  high  level  of  resolution. 
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Satellite 


Figure  1:  Typical  scenario  for  satellite  observation 


Table  1:  Calculated  Separation  Distances 


Diameter  of  Primary  Mirror  (m) 

Or  (rad) 

L  (m) 

3 

2.85e-7 

10.2 

10 

8.54e-8 

3.06 

100 

8.54e-9 

0.31 

Traditional  space-based  optical  systems  are  constructed  of  highly  polished  monolithic  glass. 
These  mirrors  are  expensive  to  produce  in  large  sizes  and  the  polishing  that  is  required  over  the 
entire  surface  is  very  time  consuming.  Additionally,  the  weight  of  traditional  mirrors  makes  them 
costly  to  lift  into  space  with  the  launch  craft  that  are  presently  available.  Another  major  limitation  is 
that  glass  mirrors  cannot  be  stored  in  any  compact  configuration  for  launch,  therefore  the  maximum 
diameter  is  limited  to  what  will  fit  in  the  cargo  bay  of  a  launch  vehicle,  presently  about  3  meters. 

One  obvious  solution  is  a  change  of  material.  If  a  lightweight  material  could  be  substituted  for 
the  heavy  glass  of  traditional  optical  systems,  the  weight  problem  could  be  eliminated.  A  further 
advancement  may  involve  the  use  of  a  thin,  deformable,  lightweight  material.  Such  a  material  would 
allow  lightweight,  compact  storage  of  the  mirror  for  the  space-lift  phase.  The  mirror  could  then 
be  deployed  to  its  final  size  once  positioned  in  its  predetermined  orbit.  Of  course,  such  a  flexible 
material  would  most  likely  need  active  control  to  create  the  desired  surface  shape  and  maintain  a 
high  quality  optical  surface. 

An  optical  quality  surface  must  have  a  surface  flatness  below  the  order  of  magnitude  of  the 
light  one  hopes  to  reflect.  Visible  light  has  a  spectrum  of  400  nm  -  700  nm.  With  advances  in 
compensated  optical  systems  and  adaptive  optics,  this  requirement  can  be  reduced  by  one  order  of 


2 


magnitude.  Therefore,  high  optical  quality  is  defined  to  be  on  the  order  of  0.5  pm  and  compensated 
optical  quality  (or  simply  optical  quality)  is  assumed  to  be  approximately  5  pm. 

The  solutions  described  here  require  advancements  in  a  diverse  range  of  engineering  disciplines. 
Materials  and  manufacturing  researchers  are  finding  new  and  innovative  ways  to  use  lightweight 
materials  as  a  replacement  for  traditional  mirrors.  Electrical  and  control  systems  engineers  are 
researching  the  necessary  stabilization  and  surface  shaping  methods  which  will  be  incorporated  into 
this  new  class  of  mirrors.  Optical  engineers  are  attempting  to  tie  all  of  these  advancements  together 
to  create  a  working  optical  system  which  will  have  unprecedented  imaging  properties. 

1.2  Scope 

Lightweight,  thin  deformable  optical  structures  are  the  key  to  achieving  the  goal  of  large 
aperture  imaging  systems  in  space.  One  of  the  largest  hurdles  preventing  this  advancement  is  the 
problem  of  stabilization  and  control  of  the  membrane  surface.  Previous  efforts  at  the  Air  Force 
Institute  of  Technology  (AFIT)  have  demonstrated  the  validity  of  the  use  of  piezo-actuated  control 
patches  etched  onto  the  under-surface  of  the  mirror  [21].  This  research  involved  testing  the  static 
deformation  pattern  of  a  5  inch  silicone  coated  PVDF  mirror.  Up  to  seven  control  patches  were 
etched  into  the  PVDF  layer  as  actuators.  The  scope  of  the  present  research  seeks  to  extend  the 
previous  effort  by  developing  a  data  acquisition  system  capable  of  sampling  the  optical  wavefront 
at  rates  as  high  as  30  Hz  and  using  it  to  characterize  static  and  dynamic  actuation  of  the  mirror 
surface.  The  scope  of  the  project  is  also  expanded  to  include  the  development  of  a  methodology  for 
a  calculation  of  membrane  tension,  verification  of  previously  developed  static  actuation  results  and 
improvement  in  membrane  fabrication  methodology. 

A  data  acquisition  system  was  developed  to  work  with  existing  wavefront  measurement  equip¬ 
ment.  This  acquisition  system  was  necessary  in  order  to  establish  a  link  between  the  motion  of 
the  optical  surface,  the  acquisition  equipment  and  the  analysis  and  control  equipment.  This  new 
system  enabled  detailed  measurement  of  the  optical  surface  subject  to  a  sinusoidal  forcing  function. 
The  improved  acquisition  system  is  also  capable  of  quickly  taking  data,  which  is  used  for  statistical 
analysis  of  static  deflection  results. 

A  scanning  laser  vibrometer  was  used  to  characterize  the  membrane-like  optical  structures  over 
a  full  range  of  frequencies.  This  data  was  instrumental  in  a  prediction  of  the  tension  in  the  PVDF 
membrane  layer. 

The  methodology  of  fabrication  was  also  improved  to  decrease  the  areal  density  of  the  membrane¬ 
like  optical  structures,  this  lower  weight  is  a  result  of  a  decreased  thickness  of  the  silicone  surface 
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coating.  A  decrease  in  thickness  creates  a  more  membrane-like  structure,  which  allows  standard 
membrane  theory  to  better  predict  the  dynamics. 

1.3  Problem 

The  primary  concern  with  previous  work  involves  the  data  acquisition  system.  Data  could  be 
taken  only  in  single  frames,  and  had  to  be  post-processed  for  an  accurate  description  of  the  surface. 
An  active  control  system  requires  that  the  data  acquisition  system  be  able  to  grab  data  and  provide 
feedback  for  the  controller  in  a  timely  manner,  typically  many  times  a  second.  Therefore  the  primary 
focus  of  this  research  is  to  create  a  data  acquisition  system  which  can  be  used  with  Matlab® .  This 
will  enable  future  researchers  to  use  Simulink®  and  dSpace®  control  hardware  to  develop  wavefront 
control  algorithms  for  the  optical  surface. 

At  present,  the  data  acquisition  system  is  used  to  investigate  the  behavior  of  the  optical  surface 
subject  to  static  and  dynamic  loading.  Multiple  data  samples  on  a  statically  loaded  membrane  could 
provide  statistical  information  about  the  mean  and  variance  of  the  optical  surface.  Additionally,  the 
optical  surface  could  be  observed  in  ‘real’  time,  subject  to  sinusoidal  forcing  functions. 

A  secondary  concern  involves  the  tension  on  the  membrane.  An  accurate  tension  is  essential 
to  the  development  of  a  finite  element  model  for  this  optical  system.  Currently,  the  tensioning  of  the 
membrane  is  done  on  a  tensioning  apparatus  where  the  load  is  applied  via  uncalibrated  mechanical 
torsion.  Modal  analysis  methods  were  utilized  to  estimate  the  tension  of  the  membrane,  which 
provides  critical  information  for  use  by  the  numerical  modeling  researchers. 

Taking  only  one  measurement  at  a  time  will  give  a  general  idea  of  the  results  that  can  be 
obtained,  but  a  proper  statistical  analysis  of  the  system  requires  that  data  be  collected  quickly.  This 
will  enable  a  mean  and  standard  deviation  calculation  and  this  information  clearly  gives  a  more 
accurate  description  of  the  system.  The  newly  developed  data  acquisition  system  was  also  used  to 
perform  this  detailed  analysis. 

1-4  Summary  of  Thesis 

The  first  chapter  is  devoted  to  a  general  introduction  to  the  subject  of  this  research  and  an 
outline  of  the  problems  that  were  approached  here  in.  In  the  second  chapter,  a  review  of  the  relevant 
literature  is  provided.  This  chapter  is  intended  to  further  develop  an  understanding  of  the  problem 
and  quantify  some  of  the  research  efforts  that  have  been  made  regarding  this  issue.  As  this  is  a  very 
active,  multidisciplinary  area  of  research,  this  review  is  by  no  means  exhaustive.  Simply  the  most 
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relevant  information  to  this  project  is  provided  and  the  reader  is  encouraged  to  investigate  other 
references  at  his  leisure. 

The  third  chapter  describes  the  design  and  fabrication  of  the  membrane-like  optical  structures. 
This  summarizes  the  iterative  process  that  was  used  to  produce  the  optical  test  subjects. 

Chapter  four  provides  a  complete  description  of  the  modal  analysis  of  the  test  articles,  which 
was  performed  with  the  scanning  laser  vibrometer.  Included  in  this  chapter  is  a  detailed  outline  of 
the  theory  of  vibrating  circular  membranes,  an  analysis  and  prediction  of  the  membrane  tension  and 
an  analysis  and  comparison  of  the  primary  modes  of  the  membrane  test  subjects. 

The  fifth  chapter  includes  a  complete  description  of  the  wavefront  imaging  system.  The  optical 
test  setup  and  data  acquisition  system  are  described  and  test  results  for  the  static  and  dynamic 
loading  of  the  5  inch  optical  structure  are  provided. 

The  final  chapter  outlines  some  conclusions  and  recommendations  for  further  study.  Great 
care  have  been  taken  to  create  usable  procedures  for  the  test  equipment.  This  information  along 
with  much  of  the  analysis  code  is  provided  to  future  researchers  in  the  Appendix. 
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II.  Review  of  the  Relevant  Literature 


2. 1  Overview 

There  have  been  many  recent  research  advances  which  have  made  the  daunting  task  of  placing 
a  large  optical  quality  mirror  in  space  and  maintaining  a  precise  shape  across  its  entire  surface  seem 
more  attainable.  Two  areas  of  research  are  currently  producing  the  most  promising  potential  for 
meeting  this  ambitious  goal.  Advancements  in  materials  and  manufacturing  techniques  have  greatly 
reduced  the  areal  density  of  rigid,  lightweight  optical  structures.  Similarly,  improvements  in  piezo- 
actuated  polymer  control  and  boundary  control  of  membrane-like  structures  has  paved  the  way  for 
an  order  of  magnitude  decrease  in  areal  densities  by  allowing  the  use  of  flexible,  thin,  deployable 
structures  to  create  new  lightweight  mirror  designs. 

The  areal  density,  p areal,  is  useful  as  a  comparison  of  mass  per  surface  area  for  different  optical 
structures.  It  is  defined  in  Equation  3.  Here  p,  to,  h  and  V  are  defined  as  density,  mass,  height  and 
volume,  respectively. 


Pareal  —  ph  y  h  (3) 

These  large  decreases  in  areal  density  are  necessary  if  a  large  aperture  space-based  optical 
system  is  to  be  realized.  Various  U.S.  Government  agencies  have  proposed  mission  objectives  which 
require  optical  surfaces  up  to  100  meters  in  diameter  to  be  constructed  [15]. 

Decreasing  the  areal  density  produces  many  design  problems  that  have  not  been  addressed  with 
traditional  optics.  Structural  vibration,  durability,  deployability  and  robustness  are  all  important  to 
consider.  Active  and  adaptive  control  algorithms  are  necessary  to  constantly  maintain  the  surface 
shape  and  quality. 

This  section  outlines  some  of  the  recent  advances,  which  have  made  this  design  problem  more 
tractable.  Two  areas  of  research  are  considered:  lightweight,  rigid  optical  structures  and  membrane¬ 
like  structures. 


2.2  Lightweight,  Rigid  Optical  Structures 

Currently,  the  standard  material  used  to  produce  space-based  optical  systems  is  monolithic 
glass,  which  can  be  polished  to  a  very  suitable  optical  surface.  This  polishing  process  is  time 
consuming  and  costly.  The  systems  are  also  heavy,  often  exhibiting  areal  densities  in  the  hundreds 
(approximately  200  for  the  primary  mirror  on  the  Hubble  Space  Telescope).  These  traditional 
systems  are  relatively  robust  and  do  not  require  the  amount  of  active  stabilization  that  a  lighter 
weight,  flexible  membrane-like  system  does. 
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The  space  shuttle  serves  as  an  appropriate  benchmark  for  the  size  and  lifting  capacity  of  a 
typical  modern  day  launch  vehicle  used  to  lift  a  telescope  into  orbit.  The  maximum  cargo  size  for 
the  space  shuttle  is  approximately  3.25  meters  and  it  can  accommodate  approximately  28,800  kg. 
These  numbers  will  be  assumed  for  future  comparisons  of  areal  density.  As  an  example,  a  monolithic 
glass  mirror  which  maximized  both  the  weight  and  diameter  specification  would  have  an  areal  density 
of  approximately  3,472  ^4,  as  shown  by  Equation  4. 


Pareal  — 


28800 


3472^ 

mz 


(4) 


Similarly,  an  optical  system  with  a  primary  mirror  aperture  of  100  meters  would  need  an  areal 
density  of  about  3.67  ^4  to  be  lifted  by  the  space  shuttle. 

One  approach  is  to  combine  the  positive  aspects  of  monolithic  glass  mirrors  with  advanced 
material  processing  and  fabrication  methods  to  create  thinner,  light  weight,  rigid  optical  structures. 
These  structures  help  to  lower  the  launch  weight  of  the  optical  system,  but  do  nothing  to  increase 
the  aperture  size  of  the  mirror. 

Researchers  at  AFRL  are  looking  into  ways  to  fabricate  a  mirror  using  a  spin  casting  technology 
to  form  the  reflective  surface  in  a  shape  close  to  that  which  it  will  be  utilized.  These  methods 
emphasize  the  use  of  stress  coatings  applied  to  the  optical  surface  coupled  with  boundary  control 
techniques  to  ensure  that  an  optical  quality  surface  is  maintained.  Other  approaches  considered  by 
this  group  include  fabrication  of  the  rigid  support  structure  or  mirror  surface  once  deployed  on  orbit. 
This  process  may  involve  the  use  of  shape  memory  alloys  or  electro-active  materials  (such  as  PVDF) 
to  deploy  the  membrane-like  structure  to  its  usable  state  [4]. 

Another  approach  is  to  simply  replace  the  heavy  support  structure  of  traditional  monolithic 
glass  mirrors  with  light  weight  alternatives.  Researchers  at  Eastman  Kodak  Co.  have  produced  a 
1.4  meter  mirror  with  an  areal  density  of  around  12  -^4.  This  mirror  uses  a  carbon  fiber  reinforced 
polymer  matrix  composite  reaction  structure  to  support  the  optical  surface  [15].  It  is  constructed 
in  the  shape  of  a  hexagon,  so  that  it  can  be  fit  with  other  similar  pieces  to  form  a  larger  primary 
mirror. 


AFRL  Materials  Directorate  researchers  are  also  pursuing  glass  or  ceramic  syntactic  foam 
[15].  These  are  structures  which  are  made  from  glass,  ceramic  or  polymer  microballoons.  The 
microballoons  are  suspended  in  a  polymer  resin  matrix.  These  structures  are  stiff,  very  lightweight 
and  can  be  coated  with  reflective  material  to  form  an  optical  quality  surface. 
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2.3  Lightweight,  Deformable  Membrane-like  Structures 

While  the  diameter  of  lightweight,  rigid  mirrors  are  restricted  by  the  cargo  bay  dimensions 
of  today’s  launch  vehicles,  an  optical  surface  which  could  be  compressed  in  some  way  could  be 
substantially  larger.  For  this  reason,  several  ideas  are  being  considered  to  fold,  roll  or  otherwise 
store  flexible,  membrane-like  structures  for  launch.  These  structures  are  then  deployed  in  orbit. 

The  flexible  nature  of  these  structures  necessitates  some  means  of  tensioning  and  control  in 
order  to  maintain  an  optical  quality  surface.  The  membrane- like  structure  can  either  be  controlled 
with  subtle  deformations  around  the  boundary  or  by  control  actuation  applied  to  the  interior  of  the 
membrane  surface. 

2.3.1  Boundary  Manipulation  Control  Methods.  One  method  of  deploying  a  lightweight 
optical  structure  which  has  been  studied,  is  to  enclose  it  with  an  inflatable  structure  and  simply 
pump  positive  pressure  into  the  system.  The  pressure  of  the  gas,  shape  of  the  structure  and  other 
actuation  methods  could  be  used  to  control  the  flexible  surface.  Figure  2  gives  an  illustration  of 
such  a  system.  One  major  drawback  of  this  system  is  that  the  transparent  canopy  would  be  difficult 
to  manufacture.  To  accommodate  much  of  the  visible  spectrum,  the  total  thickness  and  thickness 
variation  across  the  structure  would  have  to  be  kept  to  very  tight  tolerances  [4]. 
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Figure  2:  Lenticular  Optics  System  [4] 

Another  design  approach  is  to  use  a  deformable  boundary  to  enclose  the  membrane-like  struc¬ 
ture.  This  approach  has  been  validated  by  the  inflatable  antenna  experiment,  which  was  carried  out 


by  NASA  in  May,  1996.  Systems  similar  to  this,  with  piezoelectric  control  patches  applied  to  the 
outer  ring  have  been  studied  [10]. 

Work  has  also  been  done  to  develop  analytical  models  to  represent  these  boundary  actuated 
flexible  membrane  optical  structures  [9]  and  [6] . 


S77E5011  1998:05:2008:03:03 

Figure  3:  Inflatable  Optical  System  [8] 

Experiments  have  been  conducted  to  study  the  effects  of  stretching  a  membrane  structure  over 
a  curved  rail  system.  The  rails  can  be  actuated  by  forces  and  moments  applied  to  their  ends.  The 
deformed  rails  alter  the  shape  of  the  membrane  and  offer  some  measure  of  control  [1] . 

2.3.2  Interior  Manipulation  Control  Methods.  Results  from  a  test  conducted  at  the  Uni¬ 
versity  of  Kentucky  indicate  that  a  non-contact  electron  gun  can  feasibly  be  used  to  deform  a  piezo¬ 
electric  mirror  [5].  The  test  subject  for  these  experiments  was  a  piezoelectric  bimorph  (PZT-5H)  2 
mm  thick  and  2  inches  in  diameter.  Deflections  on  the  order  of  3  /iiri  positive  deflection  and  2  /jm 
negative  deflection  were  observed.  This  method  of  actuation  requires  complicated  equipment  in  a 
controlled  environment,  therefore  questions  regarding  the  scaling  of  this  approach  are  legitimate. 

Another  approach  is  to  etch  control  patches  directly  onto  the  piezoelectric  material.  Research 
along  these  lines  has  been  done  at  the  Air  Force  Institute  of  Technology.  It  has  been  shown  that 
static  deformations  on  the  order  of  3.3  /jm  are  possible  using  5  inch  membranes  [22].  For  these  tests, 
membranes  with  a  thickness  of  4  to  6  mm  were  produced.  Membranes  were  constructed  of  etched 
PVDF  material.  Due  to  the  control  patch  placement  directly  on  the  optical  membrane,  layers  of 
silicone  substrate  are  needed  to  smooth  the  surface  and  create  and  optical  quality  mirror. 
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Analytical  modeling  of  piezo-actuated  material  has  also  been  an  area  of  research.  Researchers 
at  NASA  have  developed  an  analytical  solution  to  deform  a  flat  piece  of  electro-active  polymer  into 
a  predefined  deformed  state  and  have  verified  the  results  computationally  [3] . 

Mirrors  with  actuation  patches  require  less  complicated  equipment  than  non-contact  methods, 
however  etching  the  control  pattern  without  imprinting  onto  the  optical  surface  can  be  a  problem. 
No  one  method  is  a  clear  best  alternative,  each  possesses  new  problems  that  will  require  additional 
research  efforts  to  solve. 

Many  potential  solutions  are  currently  under  investigation  within  the  scientific  community. 
Previous  work  completed  at  AFIT  demonstrated  great  promise  for  an  advancement  in  the  area  of 
control  of  a  tin  deformable  PVDF  optical  surface.  Thus,  the  current  research  seeks  to  validate  and 
extend  this  valuable  effort  by  creating  an  acquisition  system  which  will  enable  data  to  be  processed 
in  a  more  timely  manner.  This  development  is  critical  if  any  control  algorithms  are  to  be  developed 
for  control  or  stabilization  of  the  mirror  surface  in  the  future. 
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III.  Membrane  Design  and  Fabrication 

3. 1  Overview 

The  design  and  fabrication  process  for  the  test  articles  is  a  very  important  phase  of  the  project. 
Each  mirror  must  be  produced  by  hand,  which  is  a  time  consuming  process.  There  are  many  ways 
that  each  mirror  could  be  produced,  but  simply  changing  one  small  technique  during  fabrication  may 
produce  separate  mirrors  with  different  dynamics.  Great  care  was  taken  to  follow  the  procedures 
outlined  below  for  every  mirror  that  was  fabricated.  Any  deviation  from  this  is  noted  and  explained. 
By  closely  following  the  outlined  procedures,  wasted  time  and  resources  was  minimized  and  mirrors 
with  similar  dynamics  are  produced  for  testing. 

3.2  Control  Pattern  Design 

Three  control  patterns  are  used  for  this  research,  1-patch,  5-patch  and  7-patch.  The  following 
sections  describe  the  reasons  that  each  control  pattern  was  selected.  A  description  of  the  dimensions 
of  each  control  patch  is  also  provided. 

The  width  of  all  control  leads  was  set  at  |  inch  (3.2e-3  m).  This  dimension  was  chosen  to 
minimize  the  area  of  the  control  lead,  while  still  enabling  fabrication  by  hand. 

Design  considerations  for  the  5  and  7-patch  mirrors  included  the  ability  to  actuate  both  sym¬ 
metric  and  antisymmetric  portions  of  the  surface,  patch  spacing  and  ease  of  fabrication.  The  final 
design  satisfied  all  design  considerations.  The  mirror  surface  can  be  actuated  symmetrically  with 
only  the  central  control  patch,  all  outer  patches  or  all  patches  acting  in  unison.  The  mirror  can  be 
actuated  anti-symmetrically  through  the  use  of  any  combination  of  the  outer  control  patches. 

The  control  patches  were  carefully  spaced,  to  provide  the  greatest  actuation  area,  while  still 
maintaining  a  minimal  patch  spacing.  The  minimum  patch  spacing  was  determined  by  fabrication 
constraints  (|  in)  and  serves  to  provide  insulation  between  charge  on  different  patches. 

The  control  patterns  were  designed  with  only  one  set  of  concentric  patches  because  this  is 
sufficiently  complex  to  demonstrate  both  symmetric  and  antisymmetric  motion,  while  it  is  simple 
enough  to  fabricate  by  hand. 

The  outer  most  edge  of  the  patch  design  is  located  |  inch  (1.9e-2  m)  from  the  aluminum  ring 
support,  which  serves  two  purposes.  By  using  this  inner  portion,  the  slightly  distorting  effects  of 
the  meniscus  can  be  minimized.  The  meniscus  describes  the  slightly  raised  level  of  the  silicone  near 
the  aluminum  wall,  which  is  formed  during  the  cure  process.  Additionally,  the  wavefront  sensing 
equipment  has  an  optimal  illumination  of  approximately  4  inches,  so  by  localizing  the  control  patches 
to  this  diameter,  the  maximum  actuation  effect  can  be  observed. 
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The  control  patch  pattern  is  designed  only  to  be  etched  onto  one  side  of  the  PVDF  material. 
By  placing  the  pattern  onto  the  underside  of  the  mirror  surface  and  mounting  the  aluminum  ring 
onto  the  other,  a  layer  of  insulation  is  maintained.  This  prevents  shorting  of  one  patch  to  another 
through  the  aluminum  ring  or  by  other  means.  It  is  important  to  note  that  PVDF  can  be  actuated 
on  both  the  top  and  the  bottom,  which  may  increase  control  input  options  for  future  researchers. 
Double  sided  actuation  would  require  the  development  of  an  improved  membrane  support  structure, 
which  would  maintain  an  insulated  boundary  between  control  patches. 


Aluminum  Tensioning  Ring 


Figure  4:  Schematic  of  1-Patch  Design 

3.2.1  1-Patch,  12  inch  Membrane.  The  1-patch  etch  design  (see  Figure  4)  was  used 
primarily  as  a  proof  of  concept.  The  diameter  of  the  central  control  circle  is  3  inches  (07.62e-2  m) 
and  the  length  of  the  control  lead  is  7  inches.  The  length  of  the  lead  is  determined  by  the  radius 
of  the  tensioning  apparatus,  which  is  6  inches.  The  diameter  of  the  control  patch  was  decided  upon 
in  a  somewhat  arbitrary  manner.  The  patch  is  large  enough  to  actuate  the  tensioned  membrane  for 
the  vibration  phase  of  testing,  yet  small  enough  not  to  be  cumbersome  during  fabrication. 

This  simple  design  was  chosen  as  a  starting  point  in  order  to  serve  as  a  learning  device  for 
some  of  the  more  complex  control  patterns  which  needed  to  be  fabricated.  This  control  pattern 
demonstrated  that  it  is  possible  to  etch  thin,  clean,  straight  lines  for  the  control  leads  by  hand. 

3.2.2  5-Patch,  5  inch  Membrane.  The  5-patch  control  pattern  was  developed  primarily  for 
use  as  a  final  test  mirror.  The  5-patch  design  has  larger  control  patches  than  the  7-patch,  therefore 
it  offers  a  larger  actuation  area  as  a  test  subject. 
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The  dimensions  of  the  5-patch  test  subject  are  outlined  in  Figure  5.  For  identification,  each 
control  patch  is  given  a  number  label  1  to  5  as  shown  in  this  diagram.  The  area  of  the  central 
control  patch  is  0.7854  in2  (5.067e-4  m2),  this  represents  approximately  4.0%  of  the  mirror  surface 
area.  The  area  of  the  outer  control  patches  is  1.19  in2  (7.6826e-4  m2),  which  is  approximately  6.0% 
of  the  mirror  surface  area. 

Two  5-Patch  mirrors  were  etched  for  this  project,  however  only  one  was  given  a  silicone  surface 
coating.  These  structures  will  be  referred  to  as  5-Patch  A  and  5-Patch  B,  where  only  5-Patch  A  was 
coated. 


Figure  5:  Schematic  of  5-Patch  Design 

3.2.3  7-Patch,  5  inch  Membrane.  The  7-patch  control  pattern  was  developed  primarily  for 
use  as  a  final  test  mirror.  The  7-patch  design  has  more  control  patches  than  the  5-patch,  therefore 
it  offers  more  actuation  options  as  a  test  subject. 

The  dimensions  of  the  7-patch  test  subject  are  outlined  in  Figure  6.  For  identification,  each 
control  patch  is  given  a  number  label  1  to  7  as  shown  in  this  diagram.  The  area  of  the  central 
control  patch  is  0.7854  in2  (5.067e-4  m2),  this  represents  approximately  4.0%  of  the  mirror  surface 
area.  The  area  of  the  outer  control  patches  is  0.866  in2  (5.587e-4  m2),  which  is  approximately  4.4% 
of  the  mirror  surface  area. 
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Two  7-Patch  mirrors  were  etched  for  this  project  and  both  membranes  were  coated  with  silicone 
for  wavefront  testing.  These  structures  will  be  referred  to  as  7-Patch  A  and  7-Patch  B. 


Figure  6:  Schematic  of  7-Patch  Design 


3.3  PVDF  Material 

The  PVDF  material  used  for  this  research  is  manufactured  by  Measurement  Specialties,  Inc. 
It  is  52 fj,m  thick  and  can  be  purchased  in  an  unprocessed  14  inch  wide  roll. 

The  material  is  typically  rolled  during  processing  and  a  usable  center  patch  is  cut  out  before 
delivery  to  the  client.  It  was  determined  that  the  post-processed  PVDF  material  is  too  small  to  fit 
AFIT’s  tensioning  apparatus.  The  supplier  agreed  to  provide  the  roll  of  PVDF  material  without 
this  post-processing  step.  Only  the  center  of  the  PVDF  material  is  utilized  once  it  is  etched  and 
epoxied  to  the  aluminum  ring,  thus  insuring  a  high  quality  piece  of  PVDF  material  for  testing. 

Each  piece  of  PVDF  has  a  thin  layer  of  conductive  copper  and  nickel  coated  on  each  side,  see 
the  illustration  in  Figure  7. 

PVDF  material  is  a  piezoelectric  polymer  which  responds  by  expanding  when  a  voltage  is 
applied  to  it.  Conversely,  if  the  material  is  deformed,  it  will  produce  a  voltage.  PVDF  also  exhibits 
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Figure  7:  Schematic  of  PVDF  Material  Layers 

a  pyro-electric  effect,  where  upon  it  will  produce  a  voltage  due  to  a  change  in  temperature.  For  this 
reason  the  material  is  widely  used  in  the  sensors  community. 

The  deformation  of  PVDF  material  is  a  result  of  the  attraction  (or  repulsion)  of  internal  dipoles 
due  to  an  applied  electric  field.  The  PVDF  material  is  anisotropic,  therefore  it  will  not  expand  the 
same  in  all  directions.  A  positive  polarity  will  cause  the  PVDF  material  to  expand  in  length  and 
contract  in  width  and  height.  A  negative  polarity  will  cause  the  PVDF  material  to  contract  in  length 
and  expand  in  width  and  height.  The  coordinate  system  for  the  PVDF  material  has  been  defined  in 
Figure  8. 


Figure  8:  Definition  of  PVDF  Coordinate  System 


A  L  = 


Ld31v 

t 


(5) 


In  Equation  5,  L,  A L,  d3i,  v  and  h.  are  length,  length  change,  piezoelectric  coefficient  for  length 
direction,  applied  voltage  and  thickness,  respectively. 


Aw  =  ^ 
t 


(6) 
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Where  w,  Aw  and  d32  are  width,  width  change  and  piezoelectric  coefficient  for  width  direction. 


At  = 


td33v 

t 


d33v 


(7) 


In  Equation  7  d33  is  the  piezoelectric  coefficient  for  the  3-direction  (height). 

PVDF  material  cannot  generally  be  used  to  create  large  displacements,  but  it  can  be  used  over 
a  large  range  of  frequencies  (lelO-3  Hz  to  lel09ffz).  Table  2  provides  some  of  the  useful  coefficients 
that  apply  to  the  PVDF  material.  This  data  is  taken  directly  from  the  Measurement  Specialties, 
Inc.,  Piezo  Film  Sensors  Users  Manual  [11]. 


Table  2:  PVDF  Coefficients  and  Material  Properties 


Symbol 

Parameter 

Value 

Units 

h 

Thickness 

52 

fim  (micron,  lO-9  m) 

d3i 

Piezo  Strain  Constant  (length) 

23 

1(T12  f 

d33 

Piezo  Strain  Constant  (height) 

-33 

lO-12  m 

Pm 

Mass  Density 

1.78 

103  — 
m 

3-4  Mirror  Etching 

The  fabrication  process  utilized  for  the  optical  membrane-like  structures  is  an  extension  of  that 
used  in  previous  research  [21].  The  clean  PVDF  material  is  stored  on  a  thick  cardboard  roller  (to 
prevent  creases  in  the  material)  with  the  positive  side  rolled  on  the  inside. 

The  positive  side  of  the  PVDF  material  was  marked  with  permanent  masking  marker  to  avoid 
confusion.  This  will  help  to  make  sure  that  the  aluminum  ring  is  epoxied  to  the  correct  side  of  the 
PVDF  sheet.  It  will  also  help  to  predict  the  positive  or  negative  displacement  of  the  final  mirror 
surface  under  an  applied  voltage. 

To  track  the  orientation  of  the  material  during  fabrication,  three  parallel  lines  are  marked  in 
the  corner  near  the  positive  side  marking.  The  lines  are  drawn  parallel  to  the  length  of  the  material, 
perpendicular  to  the  cardboard  roller  and  correspond  to  the  1  direction  in  Figure  8. 

Computer  generated  masking  patterns  were  created  and  used  to  form  the  mask  for  the  etching 
process.  The  specifics  of  these  masking  patterns  is  covered  in  Section  3.2.1  to  3.2.3.  A  razor  was  used 
to  cut  the  control  patches  out  of  the  computer  generated  stencil.  This  was  done  extremely  carefully 
to  maintain  the  exact  spacing  of  the  control  pattern  design.  The  area  that  is  to  be  the  control 
pattern  and  therefore  must  be  masked  with  a  permanent  marker  must  be  completely  removed. 

A  14  inch  piece  of  PVDF  material  is  cut  from  the  roll  of  clean  material  with  a  straight  edge 
and  the  razor  blade.  The  resulting  PVDF  piece  is  roughly  14  inches  square.  It  is  important  to 
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maintain  this  size,  because  smaller  dimensions  will  not  fit  on  the  tensioning  stand  and  cannot  be 
evenly  tensioned. 

The  center  of  the  PVDF  piece  is  found  by  touching  two  pieces  of  string  to  each  diagonal  corner. 
This  will  form  and  ‘X’,  which  will  mark  the  spot  to  center  the  membrane.  A  small  mark  is  made  at 
the  center  with  the  permanent  marker. 

The  computer  generated  template  is  then  taped  down  to  the  positive  side  of  the  PVDF  material 
with  Scotch  tape.  The  lead  for  the  center  patch  is  always  oriented  perpendicular  to  the  1  direction. 
Eventually,  the  excess  PVDF  material  will  be  trimmed  away  and  this  orientation  will  replace  the 
positive  side  and  parallel  line  markings  which  were  made  earlier.  The  middle  control  patch  is  then 
matched  to  the  center  of  the  PVDF  material. 

The  PVDF  material  is  masked  with  the  permanent  marker  by  tracing  out  the  control  patches 
and  filling  in  the  visible  regions  of  the  PVDF.  Great  care  is  taken  not  to  smear  the  permanent  marker 
while  it  dries.  If  the  pattern  is  deformed  in  some  way,  it  can  be  cleaned  with  isopropyl  alcohol  and 
reapplied. 

The  paper  stencil  is  then  carefully  removed.  The  control  pattern  is  then  cleaned  up  by  hand 
with  the  straight  edge  and  the  permanent  marker.  Typically,  there  are  leads  that  need  to  be  straight¬ 
ened  and  rounded  edges  that  need  smoothing  by  hand.  Occasionally,  a  permanent  mark  was  removed 
with  the  alcohol  from  a  region  not  supposed  to  be  masked.  Great  care  was  taken  not  to  change  the 
dimensions  of  the  control  pattern. 

The  PVDF  material  was  then  etched  with  Ferric  Chloride  by  simply  dipping  a  cotton  ball  in 
the  etchant  solution  and  wiping  away  the  exposed  copper  and  nickel  on  the  PVDF  surface.  The 
reaction  is  stopped  by  applying  water  to  the  etched  region.  The  PVDF  is  then  dried  by  patting  with 
a  soft  paper  towel  and  hanging  to  dry. 

The  PVDF  sheet  was  then  flipped  upside  down,  and  small  patches  of  the  material  were  etched. 
These  patches  serve  as  windows  (see  Figure  9).  Each  window  is  strategically  selected  to  allow  the 
placement  of  the  6  inch  aluminum  ring  onto  the  PVDF  during  the  epoxy  phase.  This  method  is  a 
crucial  step  in  the  centering  of  the  control  pattern  in  the  middle  of  the  ring.  If  the  control  etching  was 
to  be  off  center,  the  dynamics  of  the  membrane  would  suffer  accordingly  and  many  of  the  theoretical 
symmetry  assumptions  which  were  made  would  not  be  valid. 

Finally,  the  permanent  marker  is  removed  with  the  alcohol,  leaving  the  etched  control  patch 
design  on  the  sheet  of  PVDF.  Since  deformations  of  the  sheet  will  cause  charge  generation,  the  sheet 
is  stored  loosely  rolled  and  covered  in  poster  board  material. 
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Figure  9:  Windows  Etched  into  PVDF  Material 
3.5  Mirror  Tensioning  and  Epoxy 

Even  tension  around  the  circumference  of  the  membrane  is  important  to  the  dynamics  of  the 
final  optical  surface.  This  section  will  describe  the  methodology  which  was  used  to  tension  and 
epoxy  the  PVDF  membrane  to  the  6  inch  aluminum  ring.  There  are,  of  course,  many  other  options 
and  variations  for  these  tensioning  methods,  some  of  which  will  also  be  addressed  in  this  section  and 
some  of  which  are  beyond  the  scope  of  this  research  effort. 

The  etched  PVDF  membrane  is  placed  over  a  12  inch  tensioning  ring.  The  etched,  positive 
side  of  the  PVDF  membrane  is  placed  face  down,  so  that  the  aluminum  ring  can  be  epoxied  to  the 
smooth  negative  side  of  the  PVDF.  The  etched  windows  should  provide  the  ability  to  center  the 
aluminum  ring  over  the  control  patches,  without  the  need  to  see  the  control  patches  directly.  The 
tensioning  ring  has  a  rubber  o-ring  set  into  a  groove,  which  will  be  deformed  when  the  PVDF  is 
clamped.  The  deformation  provides  even  tensioning  around  the  circumference  of  the  membrane.  An 
aluminum  top  ring,  which  is  flat,  is  placed  on  top  of  the  PVDF.  Four  clamps  are  then  spaced  evenly 
around  the  circumference  of  the  tensioning  apparatus.  Figure  10  illustrates  this  setup. 

The  PVDF  material  is  carefully  stretched  by  hand  in  all  directions  to  remove  any  wrinkles 
which  may  be  present.  Great  care  is  also  taken  to  insure  that  the  PVDF  material  is  resting  evenly 
on  the  o-ring  around  the  entire  perimeter. 

The  clamps  are  tightened  slowly,  two  at  a  time.  Opposing  corners  are  tightened  together  and 
only  one  eighth  of  a  turn  is  given  to  the  clamp  before  switching  sets.  The  clamps  are  tightened  in 
this  manner  until  it  is  difficult  to  tighten  any  more  with  the  thumb  and  forefinger  only. 
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Figure  10:  Tension  Apparatus  Setup  Diagram 


Other  methods  of  tensioning  which  were  considered  include  measuring  the  rotation  angle  of 
the  clamp  screw,  measuring  the  torque  of  the  clamp  screw  and  placing  force  transducers  between 
the  PVDF  and  the  aluminum  top  plate  to  measure  the  compression  force.  In  the  end  this  method  of 
tensioning  was  decided  upon  because  it  had  been  used  successfully  in  the  past  [21]  and  the  tensioning 
apparatus  was  readily  available. 

One  additional  method  of  tensioning,  which  has  been  successfully  used  for  fabrication  of  optical 
membranes  by  AFRL  Researchers,  is  the  use  of  a  vacuum  system  to  pull  the  PVDF  material  down 
around  the  12  inch  ring  [10].  This  method  could  potentially  provide  very  even  tensioning,  however 
working  in  a  vacuum  introduces  complications  which  are  beyond  the  scope  of  this  project. 

An  aluminum  ring,  which  was  specifically  designed  to  allow  excess  epoxy  to  flow  away  from  the 
membrane  surface  was  used  for  the  membrane  support  structure.  Figure  11  shows  the  design  of  this 
ring.  Notice  the  shallow  epoxy  groove,  which  is  surrounded  on  either  side  by  epoxy  drain  grooves. 


,  Black  ->  Aluminum  Structure 

View^  ^ray  Shallow  Epoxy  Groove 


Figure  11:  Aluminum  Ring  Diagram 

The  grooved  area  on  the  aluminum  ring  is  cleaned  with  a  compressed  air  canister  to  remove 
any  particles  which  may  prevent  a  solid  bond  between  the  membrane  and  the  aluminum  ring.  Some 
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of  the  aluminum  rings  have  been  machined  with  wire  support  plates  on  the  top  surface  (side  opposite 
the  grooves).  These  plates  must  be  removed  with  a  hex  wrench  so  that  the  ring  will  lay  flat  on  a 
surface.  If  the  ring  is  tilted  at  all  during  cure,  the  epoxy  will  be  unevenly  distributed  around  the 
circumference  of  the  mirror,  which  may  jeopardize  the  bond  strength. 

M-bond  610  is  a  two  part  epoxy  system.  It  was  selected  for  the  project  based  on  its  ability  to 
bond  securely  to  both  aluminum  and  the  PVDF  material.  This  epoxy  was  also  selected  because  of 
its  low  viscosity,  which  allowed  the  epoxy  to  easily  fill  the  machined  grooves  of  the  aluminum  ring. 
This  is  not  a  conductive  epoxy,  so  a  grounded  lead  was  attached  to  the  negative  side  of  the  PVDF 
sheet  to  bleed  off  any  stored  charge  accumulated  during  testing. 

A  plastic  dropper  was  used  to  carefully  fill  the  inner  groove  of  the  aluminum  ring.  Once  the 
groove  is  filled,  the  epoxy  will  naturally  flow  up  to  the  shallower  middle  groove.  Epoxy  is  added  until 
the  liquid  has  completely  filled  both  the  inner  and  the  middle  grooves.  The  outer  groove  is  simply 
used  for  runoff.  If  too  much  epoxy  is  applied,  it  will  flow  into  the  outer  groove  and  to  the  outside  of 
the  aluminum  ring  through  the  slit.  This  prevents  the  epoxy  from  running  to  the  inside  of  the  ring, 
which  would  destroy  the  surface  of  the  mirror  and  distort  the  dynamics  of  the  membrane. 

The  epoxy  is  allowed  to  set  for  a  half  of  an  hour.  This  process  is  important,  because  when  the 
ring  is  turned  right-side-up  for  placement  on  the  PVDF  membrane,  the  epoxy  cannot  be  allowed  to 
run  out  of  the  grooves  onto  the  mirror. 

After  the  epoxy  has  started  to  set  up  and  the  ring  can  be  tilted  without  disrupting  the  epoxy 
in  the  aluminum  grooves,  the  PVDF  membrane  is  brushed  with  epoxy  in  the  areas  where  the  ring 
will  be  attached.  The  aluminum  ring  is  then  carefully  turned  right-side-up  and  placed  in  the  center 
of  the  control  patches.  The  windows  etched  on  the  PVDF  membrane  will  help  with  the  placement. 

The  epoxy  between  the  aluminum  ring  and  the  PVDF  is  allowed  to  cure  under  its  own  weight 
for  two  or  more  days.  The  mirror  can  then  be  carefully  removed  from  the  tensioning  apparatus.  The 
excess  PVDF  material  is  then  cut  away  from  the  6  inch  ring  with  sharp  scissors.  Care  must  be  taken 
not  to  remove  any  of  the  control  leads  and  also  to  insure  that  an  appropriate  space  is  left  for  the 
ground  wire. 

The  membrane  is  placed  in  a  150°  F  oven  for  at  lease  24  hours  to  further  solidify  the  epoxy 
cure.  Once  the  M-bond  610  is  completely  cured,  the  control  lead  wires  are  attached  to  the  control 
leads  on  the  PVDF  membrane.  The  control  lead  wires  are  covered  with  electrical  tape  to  insure  a 
solid,  isolated  connection  and  protect  against  unintentional  grounding  during  testing. 

The  membrane  is  now  appropriately  tensioned,  epoxied  to  the  six  inch  aluminum  ring  and 
ready  for  coating  with  silicone  substrate. 
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3. 6  Silicone  Coating 

After  the  etching  and  epoxy  processes  are  complete,  at  52  /im  thick,  the  material  cannot 
support  a  bending  moment  and  behaves  like  a  membrane.  However,  the  surface  of  the  mirror  is 
not  optical  quality.  Some  surface  roughness  and  inevitable  foreign  particles  have  accumulated  on 
the  mirror  surface.  The  mirror  surface  is  thoroughly  cleaned  and  coated  with  silicone  substrate  to 
reduce  the  surface  impurities  and  increase  the  optical  quality  of  the  surface.  A  smooth  coating  of 
gold  is  applied  to  the  finished  substrate  to  create  the  needed  reflective  properties  on  the  mirror. 

Before  each  mirror  is  coated,  a  few  measurements  are  taken.  In  order  to  estimate  the  thickness 
of  the  silicone  layer,  a  weight  measurement  before  and  after  coating  is  taken.  Also,  vibration  testing 
can  be  used  to  predict  the  tension  that  was  applied  to  the  PVDF  during  the  epoxy  phase  (see 
Section  3.5).  This  information  is  valuable  to  the  characterization  of  the  membrane  structure  and  is 
easiest  to  collect  before  coating. 

Coating  the  mirror  surface  with  the  silicone  substrate  requires  a  great  deal  of  care  and  expe¬ 
rience.  The  process  used  was  highly  iterative  and  experimental.  An  outline  of  the  procedure  as  well 
as  comments  about  some  of  the  problems  that  were  faced  are  discussed  below. 

The  silicone  substrate  that  was  selected  for  the  coating  process  is  GE  Silicones  RTV615.  This 
product  is  a  two  part  compound  which  requires  thorough  mixing.  The  mixture  requires  precisely  10 
parts  (by  weight)  RTV615A  resin  to  one  part  RTV615B  curing  agent.  RTV615  was  selected  because 
of  its  low  shrink  rate  and  room  temperature  cure.  It  is  clear  and  colorless,  pours  easily  after  mixing 
and  cures  in  one  week. 

The  flexibility  of  RTV615  is  similar  to  rubber  and  therefore  will  certainly  have  an  effect  on 
the  dynamics  of  the  membrane.  Silicone  substrate  can  carry  a  bending  moment,  where  the  uncoated 
PVDF  membrane  could  not.  This  alters  the  dynamics  of  the  optical  structure  from  those  of  a 
membrane  to  something  closer  to  a  thin,  deformable  plate.  In  order  to  maximize  the  membrane-like 
behavior  of  the  plate,  the  thinnest  possible  layers  were  applied.  Each  layer  was  approximately  1  to 
2  mm  thick,  however,  two  to  three  layers  were  needed  to  create  an  optical  quality  surface.  Table  3 
provides  some  useful  material  properties  for  RTV615. 

Two  configurations  are  possible  for  the  curing  process.  The  first  is  to  lay  the  PVDF  membrane 
on  a  flat  sheet  of  glass,  pour  the  silicone  substrate  and  allow  it  to  cure  with  the  weight  of  the  resin 
supported  by  the  glass  sheet.  This  configuration  would  provide  a  coat  of  even  thickness  of  silicone 
when  supported  by  the  glass.  However,  during  suspension  by  the  aluminum  support  ring  for  testing 
the  PVDF  membrane  would  have  to  support  the  weight  of  the  resin  and  may  sag.  This  would 
certainly  affect  the  optical  quality  of  the  reflective  surface.  Figure  12  illustrates  this  deformation 
problem.  The  second  method  is  to  pour  the  resin  and  allow  it  to  cure  in  the  testing  configuration  (i.e. 
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Table  3:  GE  Silicones  RTV615  Material  Properties 


Parameter 

Value 

Specific  Gravity 

1.02 

Viscosity 

4000  cps 

Shrinkage 

0.2% 

Useful  Temperature  Range 

-60°  C  to  204°  C 

Work  Time  @  25°  C 

4  hr 

Cure  Time  @  25°  C 

7  days 

Cure  Time  @  100°  C 

1  hr 

suspended  by  the  aluminum  ring).  In  this  case,  the  weight  of  the  resin  is  supported  by  the  PVDF 
during  the  cure  process.  The  silicone  will  cure  unevenly  (thicker  in  the  middle)  to  compensate  for 
the  sag  in  the  PVDF.  Figure  13  shows  an  approximation  to  the  surface  shape.  The  result  is  a  flat 
surface  when  suspended,  but  an  uneven  distribution  of  the  silicone  substrate.  Ultimately,  the  second 
approach  was  utilized  because  it  will  provide  a  flatter  surface  in  the  testing  configuration. 


Figure  12:  Glass  Cure  Configuration 


Figure  13:  Suspended  Cure  Configuration 

The  silicone  substrate  is  very  susceptible  to  contamination  during  the  curing  process,  therefore 
the  uncoated  PVDF  mirror  surface  must  be  thoroughly  cleaned  before  the  silicone  substrate  is 
poured.  An  air  jet  sprayer  is  used  to  clean  the  surface.  Any  chemicals  or  oils  on  the  surface  may 
prevent  the  silicone  from  properly  bonding  to  the  PVDF  material.  Additionally,  all  tools  and  mixing 
containers  are  cleaned  before  working  with  the  substrate.  The  lubricating  powder  on  latex  gloves 
also  seems  to  spoil  the  cure  process,  so  no  gloves  are  used  (the  substrate  is  chemically  very  neutral, 
so  a  good  washing  is  enough  to  remove  any  silicone  from  the  hands). 
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The  uncoated  PVDF  mirror  is  weighed  to  serve  as  a  baseline.  This  will  enable  an  inference 
of  the  thickness  of  the  silicone  layer,  see  Section  3.7.  The  raw  PVDF  surface  is  coated  with  GE 
Silicones  SS4120  priming  compound.  This  will  help  to  ensure  that  the  surface  is  ready  to  accept  the 
resin  with  a  tight  bond  and  that  no  air  bubbles  are  formed  on  the  bond  surface.  Simply  brushing 
on  the  primer  with  the  provided  brush  was  enough  to  coat  the  surface.  The  mirror  is  then  placed 
on  4  support  blocks  which  contact  only  the  outer  aluminum  ring.  A  sheet  of  glass  is  placed  over  the 
top  of  the  mirror  structure  to  prevent  dust  and  other  contaminants  from  collecting  on  the  primed 
mirror  surface. 

Next,  25  g  of  RTV615A  are  poured  into  the  clean  stirring  container.  2.5  g  of  RTV615B  curing 
agent  is  then  added,  ensuring  a  ten  to  one  ratio  of  resin  to  cure  agent.  The  curing  agent  and  the  resin 
are  mixed  for  ten  minutes  with  a  clean  stirring  implement.  The  mixture  is  then  poured  into  a  second 
mixing  container  and  mixed  with  a  second  stirring  implement  for  another  ten  minutes.  Finally  the 
mixture  is  poured  into  clean  test  tubes  (two  to  four  test  tubes  are  needed  for  this  amount),  placed 
in  a  centrifuge  and  spun  for  ten  minutes.  The  centrifuge  will  pull  all  of  the  entrained  bubbles  out  of 
the  mixture. 

The  resin  mixture  is  carefully  poured  onto  the  PVDF  membrane.  Pouring  is  done  in  one 
puddle,  in  the  center  of  the  membrane  and  care  is  taken  to  place  the  lip  of  the  test  tube  as  close  as 
possible  to  the  surface  of  the  membrane.  These  precautions  will  help  to  prevent  entrained  bubbles, 
which  could  work  to  the  surface  during  cure  and  destroy  to  optical  quality  of  the  mirror. 

After  the  resin  mixture  is  poured  onto  the  PVDF,  the  mirror  is  carefully  tipped  at  a  slight 
angle  to  encourage  the  mixture  to  flow  over  all  parts  of  the  exposed  surface. 

At  this  point,  some  of  the  early  mirrors  were  placed  into  a  vacuum  chamber  as  recommended 
by  the  manufacturer  [7].  This  process  was  thought  to  help  remove  any  remaining  bubbles  from  the 
mixture,  however  both  mirrors  which  went  through  this  procedure  exhibited  visible  surface  wrinkles 
after  the  resin  cured.  The  visible  surface  wrinkles  prevented  the  surface  from  being  optical  quality 
and  necessitated  an  additional  resin  coat  for  smoothing.  This  vacuum  process  seems  to  have  resulted 
in  loss  of  time  and  material  and  was  not  utilized  in  the  fabrication  of  the  most  advanced  mirrors 
(7-Patch  A  and  B  mirrors). 

The  silicone  resin  is  shielded  from  dust  by  a  glass  sheet  and  allowed  to  cure  for  seven  days 
at  room  temperature.  The  manufacturer  specifies  that  this  time  could  be  greatly  reduced  if  a  heat 
treatment  is  used,  however  the  longer  cure  time  was  used  to  allow  the  entrained  bubbles  time  to 
work  themselves  out  of  the  resin.  Accelerating  the  cure  time  may  have  resulted  in  entrained  bubbles 
and  a  disrupted  optical  surface. 
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After  the  silicone  is  completely  cured,  the  mirror  is  coated  with  a  thin  layer  of  gold  to  increase 
the  reflectivity  of  the  primary  optical  surface.  The  specimen  is  placed  in  a  vacuum  chamber  and  a 
small  amount  of  solid  gold  is  evaporated.  This  process  evenly  coats  the  mirror  which  is  suspended 
in  the  vacuum. 

The  optical  surface  quality  is  then  evaluated.  In  some  instances,  visible  wrinkles  or  surface 
imperfections  could  be  seen  with  the  naked  eye.  These  specimen  were  selected  for  recoating  and  the 
above  procedure  was  repeated  to  obtain  a  smoother  surface.  If  no  imperfections  were  visible  on  the 
surface  with  the  naked  eye,  a  small  interferometer  was  used  to  obtain  a  rough  measurement  of  the 
surface  quality.  Finally,  if  the  surface  appeared  flat  using  these  methods,  the  subject  was  tested  in 
the  AFIT  test  facility  with  the  Shack-Hartman  system  and  wavefront  measurements  were  obtained 
(Section  V). 

The  final  optical  structure  consists  of  a  base  layer  of  etched  PVDF  material,  one  or  more  layers 
of  silicone  substrate  and  a  very  thin  surface  layer  of  gold.  Details  of  the  silicone  coating  thickness 
are  closely  related  to  the  areal  density  of  the  mirror  and  are  therefore  given  in  Section  3.7. 

3.7  Areal  Density  Calculation 

The  areal  density,  p areal,  is  useful  as  a  comparison  of  mass  per  surface  area  for  different  optical 
structures.  It  is  defined  in  Equation  3: 


Pareal  —  ph  — 

V  =  h*  A  =  hnr2  (8) 

Volume  of  a  cylinder  is  defined  in  Equation  8,  where  r  and  A  are  the  mirror  radius  and  surface 
area,  respectively.  The  equation  for  areal  density  can  now  be  reduced: 

to  to  .  . 

Pareal  ~  ^  yh  =  y  (9) 

U7rrz  7rrz 

It  is  now  clear  that  the  areal  density  relates  the  mirror  mass  to  surface  area.  For  a  space 
application,  where  payload  weight  is  at  a  premium,  a  lower  areal  density  is  desirable. 

Table  4  includes  mass  measurements  for  all  the  membranes  that  were  produced.  Measurements 
were  taken  after  the  silicone  substrate  and  gold  reflective  coating  were  applied,  as  outlined  in  Sec¬ 
tions  3.2.3  through  3.2.2.  Notice  that  mirror  5-PatchB  was  not  coated  and  therefore  does  not  have 
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Table  4:  Mirror  Mass  Data 


Mirror 

Uncoated  Mass  (kg) 

Coated  Mass  (kg) 

Difference  (kg) 

7-PatchA 

0.357 

0.400 

0.043 

5-PatchA 

0.354 

0.410 

0.056 

7-PatchB 

0.355 

0.375 

0.020 

5-PatchB 

0.355 

- 

- 

Solid  Mirror 

0.370 

- 

- 

coated  measurements.  The  radius  of  all  mirrors  was  2.5  inches  (0.0635  m).  Given  this  data  we  can 
calculate  the  areal  densities  for  the  three  membrane-like  mirrors  which  were  constructed,  Table  5 
gives  this  data. 


Mirror 

Table  5:  Calculated  Areal  Densities 

Areal  Density  Substrate  Thickness  (mm) 

5-PatchA 

4.42 

4.3 

7-PatchA 

3.39 

3.3 

7-PatchB 

1.57 

1.5 

Solid  Mirror 

45.63 

- 

The  density  of  GE  Silicones  RTV615  substrate  is  known  to  be  1200  The  estimated  sub- 

.7  771° 

strate  thickness  is  then  calculated  using  Equation  3.  These  values  are  also  included  in  Table  5. 

Previous  fabrication  efforts  performed  at  AFIT  have  produced  mirrors  with  thicknesses  of  4  to 
6  mm  [22].  The  use  of  these  improved  fabrication  methods  have  resulted  in  mirror  thicknesses  of  1.5 
to  3.3  mm,  which  is  a  reduction  of  more  than  50  percent  in  thickness.  Because  the  previous  research 
was  done  with  exactly  the  same  mirror  size  and  materials,  the  areal  density  will  also  be  improved. 
Results  indicate  that  areal  densities  of  1.57  to  3.39  (^ii)  can  be  obtained  using  these  methods. 
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IV.  Modal  Analysis  of  Membranes 

4-1  Overview 

Modal  analysis  gives  essential  information  about  the  dynamics  of  a  system.  A  vibrational 
analysis  of  a  membrane  optical  structure  can  be  used  as  a  comparison  to  the  theoretical  prediction. 
This  analysis  also  allows  for  a  comparison  between  test  articles,  which  enables  verification  that  the 
fabrication  process  used  for  all  test  articles  has  been  carried  out  in  a  similar  fashion.  Finally,  this 
information  can  be  used  to  calculate  the  in-plane  tension  in  the  PVDF  layer  of  the  optical  structure, 
which  is  valuable  information  to  numerical  researchers. 

This  chapter  first  provides  a  theoretical  development  of  circular  membrane  theory  and  gives  a 
description  of  what  results  are  expected  from  modal  analysis  based  on  that  theory.  Next  a  detailed 
description  of  the  test  setup  and  data  collection  methods  is  discussed.  Then  a  calculation  of  the 
tension  in  the  membrane  is  performed. 

The  modal  analysis  data  is  then  used  to  compare  many  aspects  of  the  vibration  tests.  Attributes 
such  as  actuation  method,  input  signal  and  orientation  of  the  test  article  on  the  test  stand  are 
compared.  The  silicone  coated  test  mirror  is  also  subjected  to  modal  analysis  and  the  results  are 
presented  for  completeness. 

4-2  Theoretical  Development 

4-2.1  Modes  of  Vibration  of  a  Circular  Membrane.  Consider  a  circular  membrane,  such 
as  Figure  14,  described  by  polar  coordinates  (r,  9).  Please  note  that  the  formulation  of  this  section 
closely  follows  that  given  in  Meirovitch  [13]. 


Figure  14:  Circular  Membrane  with  Polar  Coordinates 

The  differential  equation  that  describes  membrane  motion  is  as  follows: 


X72W(r,6)  +  /32W(r,6)  =  0  (10) 
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The  Laplacian  and  j3  are  defined  as  follows: 


2  d2  Id  id2 

v  = - 1 - 1 - 

dr 2  r  dr  r2  dd2 

a2  (  w mn  \  2 

Pmn  = 

Here  W (r,  9)  is  the  displacement  of  the  membrane  surface,  u>  is  the  natural  frequency  of  the  membrane 
and  c  is  the  wave  speed.  Ultimately  the  solution  is  attained  through  a  separation  of  variables.  In 
order  to  apply  this  technique,  assume  a  solution  of  the  form: 

W(r,  9)  =  R(r)Q(9) 

Make  a  substitution  for  W(r,9)  into  Equation  10  and  rearrange. 

V2R(r)0(9)  +  (32  R(r)0(9)  =  0 


(R(r)G(0))  +  ~  (R(r)0(9))  +  ^^(R(r)G(0))  +  /32R(r)0(9)  =  0 


Next,  carrying  out  the  partials  and  factoring,  one  attains: 


0(9) 


d2R(r)  1  dR(r) 
dr2  +  r  dr 


R(r)d20(9)  , 


d92 


(32R(r)0(9)  =  0 


(11) 


Separate  11  into  two  equations,  starting  from: 


0(9) 


d2R(r)  0(9)  dR(r)  R(r)  d20(9)  2 


dr2 


r  dr  r 2  d92 


f32  R(r)0(9)  =  0 


Multiplying  by  one  obtains 


0(9) 


r2  d2R(r)  :  rO(9)  dR(r)  :  d20(9)  :  2  Q2 


R(r)  dr2  R(r)  dr 


d92 


+  r2(3  0(0)  =  0 


r\  2  f\ 

Now  solve  for  ,  obtain 


d20(9) 

d92 


=  -m 


r2  d2R(r)  |  r  dR(r)  |  2  2 

p  r 


R(r)  dr2  R(r)  dr 
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Now,  let  m2  be  defined  as  follows: 


2  __  r2  d2R{r)  |  r  dR{r)  2  2 
m  R(r)  dr2  +  R(r)  dr  +  ^  ? 

This  leads  to  two  partial  differential  equations: 


d2R{r) 

dr2 


1  dR(r) 
r  dr 


R(r)  0 


(12) 


d2Q(6) 

do2 


+  Q(9)m 2  =  0 


(13) 


Equations  12  and  13,  taken  together,  are  equivalent  to  Equation  11.  The  value  of  m2  was 
chosen  to  give  a  harmonic  equation  in  6.  The  solution  of  (13)  must  be  continuous,  therefore 


9  =  eo  ~  9  =  do  +  2nj,  j  =  1,2,3...  (14) 

So  for  any  value  of  90,  m  must  be  an  integer,  or  relation  (14)  would  be  false  and  the  solution 
could  not  be  harmonic  (which  is  a  contradiction  and  cannot  occur).  Now  solve  (12)  and  (13) 
separately. 

As  previously  discussed,  the  solution  to  (13)  must  be  periodic  and  is  therefore  best  represented 
as  a  trigonometric  function  of  the  form: 


Om(9)  =  Cimsin(m9)  +  C2mCos(m9),  m  =  0,1,2,...  (15) 

Equation  12  is  of  a  form  that  has  a  known  solution  involving  Bessel  functions: 

Rm{r)  =  C3mJm(Pr)  +  CAmYm{(3r),  m  =  0, 1, 2  . . .  (16) 

Where  Jm  is  a  Bessel  function  of  the  first  kind  and  of  order  m  and  Ym  is  a  Bessel  function  of 
the  second  kind  and  of  order  m.  Finally,  write  the  general  solution  by  combining  15  and  16. 

W(r,  9)  =  R(r)Q{9) 


28 


Where. 


W  (r,  9)  =  AlmJm(/3r)sin(rnO)  +  A2mJm(Pr)cos(mO)+ 
AsmYm{Pr)sin{nnO)  +  A4rnYm(f3r)cos(m9),  m  =  0, 1, 2  . . . 


(17) 


A\m  ^lm^m 
^2m  =  C2mC3m 
A^rri  =  ^lm^m 
^4m  “  I'Vm  Cim 

Now  the  boundary  conditions  can  be  applied,  to  obtain  the  values  of  the  constants.  Assume 
that  a  fixed  boundary  condition  exists  at  the  rim,  therefore: 

W{a,0)  =  0,  to  =  0,1, 2,... 

Since  all  points  of  the  interior  of  the  membrane  must  remain  finite  and  Bessel  functions  of 
the  second  kind  tend  to  infinity  as  the  argument  approaches  zero  (r  «  0  at  center  of  membrane), 
Azm  =  A4m  =  0  is  the  only  condition  which  will  prevent  the  displacement  from  blowing  up.  Rewrite 
Equation  17  as: 

W(r,  6)  =  AlmJm  (pr)  sin(m0)  +  A2mJm((3r)  cos  (m9),  to  =  0,1,2,...  (18) 

Equation  18  can  be  written  at  the  boundary,  r  =  a,  as: 

W(a,  9)  =  AimJm(/3a)sin(m.9)  +  A2mJm(/3a)cos(m9),  m  =  0, 1, 2, . . . 

This  can  only  be  satisfied  for  all  0,  if 

Jm{(3a)  =  0,  to  =  0,1,2,...  (19) 

This  equation  represents  an  infinite  set  of  characteristic  equations,  for  any  given  m,  there  are 
infinite  solutions,  pmn,  which  correspond  to  the  zeros  of  the  Bessel  function  of  the  first  kind,  Jm. 
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The  subscript  m  indicates  the  number  of  divisions  through  the  center  of  the  membrane  for  a 
given  mode  shape,  while  n  indicates  the  number  of  circumferential  divisions.  There  are  two  modes 
for  each  frequency,  u)mn,  but  only  one  mode  for  m  =  0.  Therefore,  the  natural  modes  are  degenerate. 
The  modes  can  be  written  as: 


W0n{r,e)  =  A0nJ0(^ry  n  =  1,2,3...  (20) 

Wmnc(r,9)  =  Amnc Jm  cos {m0),  m,n=  1,2,3...  (21) 

Wmns(r,  9)  =  AmnsJm  sin(m0),  m,n=  1,2,3... 

Where  Wmnc  is  the  cosine  component  of  the  solution  and  Wmns  is  the  sine  component.  Notice 
that  this  problem  is  self-adjoint  (i.e.  the  complex  conjugate  transpose  of  a  matrix  is  equal  to  that 
matrix),  positive  definite  (i.e.  xT Ax  >  0  where  x  £  Cn)  and  has  orthogonal  natural  modes  (i.e. 
AAT  =  1).  From  the  properties  of  orthogonality: 


PW2ondD  =  1 


>D 


Similarly, 


PAlnJo  rdrM  =  1 

npa2A2onJf  =  1 


A 


2 

on 


1 

■npa2J'l  (^a) 


PW2mncdD  =  1 


id 


p2ir  pa 


/0  ^0 


PA2mncJh  (U-^Tr)  cos2(m0)rdrd0  =  1 
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and  finally, 


7T  Oa 

-pa  A 


2 

mnc 


J2 

Jm+ 1 


( 


^ mn  \ 

- a 

C  J 


=  1 


A2  = 

mnc 


2 

7i -pa2J^+1  {^a) 


id 


pWmnsdD  = 


1 


PAmnsJm  (  sin2  (m9)rdrd9  =  1 

(^«)  =  1 


42 

^mns 


2 


7T  pa2J2z+1  (^a) 


Having  obtained  the  necessary  A  coefficients,  substitute  these  values  into  20  and  21  to  form 
the  particular  solution. 


Won(r)=  v  n=  1,2,3... 

pJWpaJi  (-^a) 


(22) 


Jm  COSlmQ'\ 

Wmnc(r,  9)  =  yA  c  ;  y  (,  m,n=  1,2,3. 

y7rpnt/m_|_i  \  c  a) 


(23) 


.  Jm  ( "?■ '  r)  sinimO) 

Wmns(r,  9)  =  "1  c  7  V  f  m,n=  1,2,3. 
^/npaJm+ 1  (^«J 


(24) 


Undamped  Response  of  a  Circular  Membrane.  The  undamped  response  of  a  circular 
membrane  will  be  developed  in  this  section.  This  development  is  included  for  completeness  and  is  not 
strictly  necessary  for  the  modal  analysis  comparisons  which  will  be  used  to  support  the  experimental 
results  in  the  remainder  of  the  chapter. 
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Figure  15:  Circular  Membrane  with  Axisymmetric  Pressure  Distribution 

Consider  a  circular  membrane,  such  as  Figure  15.  The  membrane  is  assumed  to  be  uniform 
thickness  and  clamped  at  the  boundary.  The  membrane  is  subject  to  time-dependent  uniform  pres¬ 
sure  because  this  is  similar  to  that  which  a  loud-speaker  actuator  will  produce.  This  pressure  is 
described  by: 


p{r,0,t)  =  p0f(t),  0  <r<b 


p(r ,  6,t)  =  0,  b  <r  <  a 

The  general  solution  to  the  partial  differential  equation  describing  membrane  motion  can  be 
written  in  the  following  form  through  the  use  of  the  expansion  theorem: 


W(r,e,t)  =  EE  Wmn(r,e)Vmn(t)  (25) 

m— 0 n= 1 


using  the  previously  developed  membrane  solution  (Equation  22  to  24): 


OO  OO  OO  OO  OO 

W(r,9,t)  =  ^2  Won(r,0)r]on(t)  +  E  E  Wmnc{r,  0)r)mnc(t)  +  E  E  Wmns(r,0)r)mns(t) 

n—0  m= 0  n= 1  m— 0  n= 1 


where  W(r,  9)  are  displacements  and  ij(t)  are  generalized  temporal  coordinates: 


Von  (t) 


1 

^ on 


-  /  N0n(r)sin{uj0n{t-  r))dr 


(26) 
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V  rrmc(^)  — 


Nmnc(r)sin(iL ; 

mn  (i-r))dr 


(27) 


P  mns(t) 


LUn 


Nomns(T)sin(u 

mn  (t  —  r))d,T 


(28) 


Non,  Nmnc  and  Nrnns  are  generalized  forces.  Recall  that  the  natural  frequencies  have  the  form: 


^ mn 


=  C~ 


(3n 


(29) 


Here,  T  is  the  tension  per  unit  length,  p  is  the  mass  per  unit  area  and  c  is  the  membrane  wave 
speed.  Also,  recall  that  (3mn  are  the  roots  of  the  Bessel  function,  Jm{(imn)= 0-  The  generalized  forces 
are  given  by, 


nr(t) 


(30) 


Here,  l  is  the  number  of  constrained  forces,  N{t)  are  generalized  forces  associated  with  the 
generalized  coordinates,  pit).  The  first  term  in  Equation  30  is  associated  with  the  distributed  forces, 
while  the  second  term  is  associated  with  the  concentrated  forces.  The  current  problem  as  posed 
above  involves  no  concentrated  forces,  therefore  we  are  left  with: 


na 

Won{r ,  9)p{r191  t)rdrd6 


2np0f(t)  be 

y/WpaJi  (^a)  to  on 


J ■  (t111) 


na 

Wmnc(r,  9)p(r ,  9,  t)rdrd9  =  0 


(31) 


na 

Wmna{r,  9)p(r,  9 ,  t)rdrd9  =  0 


(32) 


/  o  Jo 

It  is  indeed  reasonable  that  Equations  31  and  32  are  equal  to  zero  because  of  the  axisymmetric 
distributed  force  that  is  assumed  in  the  problem  statement. 

Substitute  back  into  26  to  28: 


Von(t) 


2np0bcJi  (^b) 
y-f>a^LJ\  (^a) 


f(t)sin(u>on(t 


r))dr 
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Vmnc{t)  —  'Hmnsif)  —  0 


Finally,  substitute  back  into  the  displacement  Equation  (25): 


t Jrt  a  ^  2pabc  ^ 

W(r,e,t)  =  — r  > 

paz  ^ 


Ji(^b)J0(^r) 


u; 


J? 


f(t)sin{won{t  -  t ))dr 


onu  1 


(33) 


So  Equation  33  is  the  solution  for  the  displacement  of  the  membrane  under  an  axisymmetric 
distributed  load  with  a  clamped  boundary  condition. 


4-2.3  Predicted  Surface  Deflection.  Equation  33  demonstrates  the  fact  that  the  roots  of 
the  Bessel  functions  play  an  important  part  in  predicting  the  the  natural  frequencies  and  surface 
deflections  of  a  circular  membrane  structure.  Figure  16  shows  the  first  few  Bessel  functions,  with  the 
zeros  clearly  labeled.  Note  that  m  corresponds  to  the  order  of  the  Bessel  function  and  n  corresponds 
to  the  nth  zeros  crossing  of  the  Bessel  function.  Table  6  provides  the  roots  of  these  functions  for 
reference. 


Figure  16:  Example  of  Bessel  Functions  Jm(x ),  for  m=0:2 
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Table  6:  Table  of  Bessel  Zeros 


n=l 

n=2 

n=3 

n=4 

n=5 

m=0 

2.4048 

5.5201 

8.654 

11.791 

14.931 

m=l 

3.8317 

7.0153 

10.173 

13.324 

16.471 

m=2 

5.1355 

8.417 

11.620 

14.796 

17.960 

m=3 

6.3799 

9.761 

13.015 

16.224 

19.409 

m=4 

7.5884 

11.065 

14.373 

17.616 

20.827 

m=5 

8.7716 

12.339 

15.700 

18.98 

22.218 

The  deformation  patterns  given  in  Figure  17  are  generated  with  a  Matlab®  code  (see  Appen¬ 
dix  B.6),  which  uses  Equation  33.  The  x,  y  and  z-axis  have  all  been  scaled  from  -1  to  1  because 
Figure  17  is  only  intended  to  illustrate  the  primary  modes  of  surface  deformation.  It  is  important 
to  note  that  the  theory  predicts  that  these  modes  will  appear  in  the  order  of  the  coefficient  values 
given  in  Table  6.  So  m  =  0,  n  =  1  is  expected  to  appear  first  in  a  frequency  profile,  followed  by 


to  =  1,  n  =  1,  then  to  =  2,  n 

=  1  and  so  on. 
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NOTE:  Plots  are  scaled  -1  <x<1 ,  -1  <y<1 ,  -1  <z<1 


Figure  17:  Predicted  Surface  Deflections  given  by  Bessel  Functions 

It  is  worth  noting  that  the  use  of  thin  membrane  theory  as  a  model  for  the  test  subjects  is  only 
a  rough  approximation  to  the  actual  dynamics  of  the  system.  There  are  two  primary  reasons  that 
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this  standard  theory  only  approximates  the  model.  First,  the  thickness  of  the  final  optical  structure 
(with  the  silicone  coating)  is  typically  a  few  millimeters.  This  structure  will  support  some  small 
bending  stress,  due  to  this  thickness  and  is  therefore  a  slight  violation  of  thin  membrane  theory  as  it 
was  developed  in  this  section.  The  second  reason  involves  the  actuation  method.  When  the  optical 
structure  is  actuated  with  the  speaker,  on  average  there  is  a  normal  force  on  the  surface.  This  is 
the  same  as  the  actuation  direction  which  was  assumed  in  the  development  of  the  thin  membrane 
theory.  However,  actuation  with  the  PVDF  patches  causes  the  patch  to  expand  in  the  plane,  the 
deformation  out  of  the  plane  is  a  result  of  bending  moments  due  to  the  thickness  of  the  optical 
structure.  This  actuation  direction  is  different  than  that  assumed  in  thin  membrane  theory  and 
may  result  in  additional  discrepancies  from  predicted  results  for  membrane  actuated  in  this  manner. 
Since  the  structure  is  overall  relatively  thin  and  flexible,  membrane  theory  was  still  utilized  to  model 
the  dynamics  of  this  system. 

4-3  Test  Setup 

The  test  setup  which  was  used  for  all  vibrations  tests  is  shown  in  Figures  18  and  19.  The 
scanning  laser  vibrometer  was  positioned  directly  above  the  test  membrane.  The  membrane  support 
ring  was  then  placed  onto  the  support  structure  with  the  test  article  suspended  in  the  middle  by 
heavy  iron  bars.  The  actuating  speaker  was  then  placed  on  supports  below  the  membrane  and  was 
rested  on  rubber  isolators,  which  helped  to  damp  vibration  in  the  test  structure.  Figure  19  illustrates 
the  setup  of  the  test  structure  with  a  close  side  view  of  the  equipment. 

The  iron  support  bars  were  carefully  wrapped  in  electrical  tape  to  isolate  the  membrane  from 
any  electrical  connection  and  vibration  in  the  structure.  Heavy  iron  bars  were  used  to  help  isolate 
the  excitation  of  the  membrane  support  structure  from  that  of  the  membrane  surface.  The  wrapped 
bars  can  be  seen  in  the  lower  image  of  Figure  18. 

This  test  setup  was  not  changed  for  any  membrane.  The  exact  location  and  orientation  of 
every  piece  of  equipment  was  marked  before  testing  and  this  orientation  was  verified  before  every 
test  was  performed.  These  efforts  helped  to  ensure  that  the  test  setup  and  vibration  in  the  support 
structure  were  similar  for  every  membrane  which  was  tested.  The  speaker  was  left  in  place  on  the 
test  setup  even  if  acoustic  actuation  methods  were  not  used. 

4-4  Data  Collection 

A  scanning  laser  vibrometer  manufactured  by  Polytec  Inc.  was  used  for  modal  analysis  of 
the  test  membranes.  A  signal  generator  integrated  with  the  laser  vibrometer  software  was  used  to 
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Figure  18:  Vibration  Test  Setup:  Global  View 


provide  an  input  signal  to  the  system  (see  Figure  21).  The  power  amplifier  used  for  all  tests  was 
manufactured  by  Active  Control  Experts  (ACX).  A  characterization  of  this  amplifier  is  given  in 
Appendix  C.l.  The  signal  was  then  used  to  power  either  a  speaker  or  the  PVDF  material  directly. 
The  scanning  laser  sampled  541  points  on  the  membrane  surface.  Figure  20  illustrates  the  pattern 
of  data  points  which  were  scanned  by  the  laser  vibrometer.  Software  provided  by  the  manufacturer 
was  used  to  view  the  data.  Frequency  response  curves,  max  deformation  plots  and  real  time  movies 
were  generated  using  this  software.  The  data  was  then  post-processed  using  Matlab®  code. 
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Figure  19: 


Vibration  Test  Setup:  Side  View 
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Figure  20: 


Data  Collection  Pattern  for  the  Scanning  Laser  Vibrometer 
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Figure  21:  Flow  Chart  of  Vibrometer  Experiment 
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4-4- 1  Test  Parameters.  Table  7  provides  an  overview  of  the  test  parameters  which  were 
used  for  data  collection  with  the  scanning  laser  vibrometer.  These  parameters  remained  the  same 
for  every  test  which  was  run. 


Table  7:  Table  of  Test  Parameters 


Parameter 

Setting 

Scan  Points 

541 

Optimal  Reading 

100  % 

Scanning  Head 

PSV-I-400  LR 

Front  End  Software 

PSV-E-400-3D  (ID) 

Acquisition  Board 

National  Instruments  PCI-4452 

Acquisition  Mode 

FFT 

Averaging 

Complex 

Average  Count 

10 

Cosine  Correction 

Active 

Bandwidth 

50  Hz  to  1.25  kHz 

Sampling  Lines 

3200 

Sampling  Frequency 

3.2  kHz 

Sample  Time 

2.56  s 

Resolution 

390.6  mHz 

Window 

Rectangular 

^.f.2  Description  of  Test  Runs  Performed.  A  summary  of  the  tests  performed  with  the 
vibrometer  is  given  below: 

•  12  Inch  Test  Article 

—  Speaker  actuation,  chirp,  standard  orientation,  optical  surface  up 

—  1-patch  actuation,  chirp,  standard  orientation,  optical  surface  up 

•  5  Inch  Uncoatecl  Test  Article 

—  5-Patch  Control  Pattern 

*  Speaker  actuation,  chirp,  standard  orientation,  optical  surface  up 

*  Speaker  actuation,  noise,  standard  orientation,  optical  surface  up 

*  PVDF  actuation  (all  patches),  chirp,  standard  orientation,  optical  surface  up 

*  PVDF  actuation  (all  patches),  chirp,  rotated  180°,  optical  surface  up 

*  PVDF  actuation  (all  patches),  chirp,  rotated  180°,  optical  surface  down 

—  7-Patch  Control  Pattern 

*  Speaker  actuation,  chirp,  standard  orientation,  optical  surface  up 

*  PVDF  actuation  (center  patch),  chirp,  standard  orientation,  optical  surface  up 
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*  PVDF  actuation  (all  patches),  chirp,  standard  orientation,  optical  surface  up 
•  5  Inch  Coated  Test  Article 

—  7-Patch  Control  Pattern  (membrane  A) 

*  Speaker  actuation,  chirp,  standard  orientation,  optical  surface  down 

—  7-Patch  Control  Pattern  (membrane  B) 

*  Speaker  actuation,  chirp,  standard  orientation,  optical  surface  down 


4-5  Tension  Calculation 

One  of  the  primary  goals  of  this  research  is  to  obtain  an  estimation  of  the  tension  entrained  in 
the  PVDF  material  during  fabrication.  The  results  of  the  vibration  tests  can  be  utilized  to  calculate 
this  tension. 

Recall  from  Section  4.2.2  that  Equation  29  gives  the  theoretical  natural  frequencies  of  the 
circular  test  article.  This  equation  can  be  rearranged  for  the  tension  as  given  in  Equation  34: 

T  =  p(UJprX  (34) 

The  zeros  of  the  Bessel  Function,  (3mn  are  known  quantities,  given  in  Table  6.  The  radius,  r, 
and  density  per  unit  area,  p,  are  known  quantities.  Therefore,  final  piece  of  the  puzzle  is  the  natural 
modes,  wmn,  for  the  structure. 

The  natural  modes  are  easily  identified  through  the  FFT’s  and  surface  deflection  plots  produced 
by  the  laser  vibrometer  testing  (see  Section  4.6.4  for  these  results).  The  task  is  to  simply  pick  a 
mode,  examine  the  surface  deflection  plot,  match  it  to  the  theoretical  surface  deflection  plots  given 
in  Figure  17  to  find  /3mn  for  that  mode  and  then  to  use  Equation  34  to  calculate  the  tension. 

This  procedure  has  been  done  for  selected  uncoatecl  test  articles  and  the  results  are  tabulated 
in  Table  8.  The  tension  calculation  presented  is  an  average  of  the  calculated  tension  for  the  first  five 
identifiable  modes  for  each  test  article. 


Table  8:  Calculated  Average  Tension 


Mirror 

Avg  Tension  — 

Std  Deviation  — 

m 

Modes  Identified  (m,n) 

12  Inch  Test  Article 

1.9128 

0.56359 

(0,1)  (1,1)  (2,1)  (1,2)  (3,1) 

5-Patch  A 

3.0049 

0.50113 

(0,1)  (1,1)  (2,1)  (1,2)  (2,2) 

7-Patch  A 

4.5874 

1.0303 

(0,1)  (1,1)  (2,1)  (0,2)  (3,1) 
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4-6  Comparison  of  Testing  Attributes 

Several  tests  were  conducted  using  membranes  of  different  orientations,  actuation  methods  and 
input  signals.  Membranes  were  also  tested  with  and  without  the  silicone  surface  coatings.  Results 
and  comparisons  that  illustrate  these  aspects  will  be  described  within  this  section.  Presentation  of 
all  results  and  data  is  simply  not  economical,  therefore  only  the  best  illustrative  examples  have  been 
chosen. 

4-6.1  Actuation  Method.  Two  different  actuation  methods  are  possible  for  the  test  articles. 
The  first  is  actuation  with  the  speaker,  mounted  in  close  proximity  to  the  membrane.  The  second 
are  the  PVDF  patches  etched  on  the  membrane  surface.  Actuation  with  the  center  patch  alone  and 
actuation  with  seven  patches  simultaneously  were  examined. 

In  order  to  test  the  dependence  of  actuation  method,  three  cases  were  tested  using  an  uncoated 
5  inch,  7-patch  membrane.  All  actuation  tests  were  run  with  the  chirp  input  signal. 


Figure  22:  FFT  and  Coherence  Plot  for  Speaker  Actuated  Case 

Figure  22  shows  the  Fast  Fourier  transform  (FFT)  and  coherence  data  for  the  speaker  actuated 
case.  The  FFT  data  represents  the  average  FFT  for  every  point  on  the  surface,  sampled  by  the 
vibrometer.  The  coherence  data  is  of  particular  interest.  Notice  that  below  150  Hz,  the  coherence 
drops  off  rapidly,  but  for  the  rest  of  the  spectrum  the  coherence  is  near  to  one.  This  indicates  that 
the  data  is  well  correlated,  therefore  it  is  more  likely  to  give  accurate  results. 

The  coherence  plot  for  the  1-patch  actuated  case  (Figure  23)  is  well  below  80%  over  the  entire 
domain  of  the  spectrum.  This  indicates  that  the  data  has  low  correlation.  This  is  expected  because 
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Figure  23:  FFT  and  Coherence  Plot  for  1-Patch  Actuated  Case 

only  a  small  patch  in  the  center  of  the  membrane  is  being  actuated.  The  time  delay  that  it  takes 
the  waves  to  propogate  to  the  exterior  of  the  membrane  may  reduce  the  correlation  of  the  data. 

Figure  24  gives  the  FFT  and  coherence  data  for  the  7-patch  actuated  case.  Here  the  coherence 
is  solidly  above  80%  after  350  Hz,  but  relatively  shakey  prior  to  that.  This  indicates  that  at  lower 
frequencies,  the  data  is  still  not  correlated  very  well,  despite  the  additional  control  patch  actuation. 
The  data  indicates  that  the  use  of  more  control  patches  gives  better  correlated  data  over  the  entire 
spectrum.  This  is  due  to  a  larger  area  of  the  membrane  being  actuated,  which  gives  a  reduced  time 
delay  and  larger  displacement. 

Table  9:  Comparison  of  Speaker,  1-Patch  and  7-Patch  Frequencies  with  Circular  Membrane  The¬ 

ory 


Theory 

Avg  Freq 

Speaker  Freq 

1-Patch  Freq 

7-Patch  Freq 

m=0,  n=l 

213.6  Hz 

213.6  Hz 

215  Hz 

218  Hz 

208  Hz 

m=l,  n=l 

331.4  Hz 

429.6  Hz 

433  Hz 

433  Hz 

423  Hz 

The  observed  frequencies  of  the  first  two  modes  are  given  in  Table  9.  These  frequencies  are 
averaged  and  this  value  is  given  in  the  ‘Avg  Freq’  column.  The  theoretical  value  is  obtained  based 
on  the  theory  presented  in  Section  4.2.3. 

Equation  29  gives  the  relationship  between  the  theoretical  natural  frequency  and  the  physical 
parameters  associated  with  the  membrane  test  subject  (i.e.  T,  p,  r  =  a).  The  zeros  of  the  Bessel 
functions  are  given  in  Table  6.  Since  the  tension,  radius  and  density  are  constant  during  this  test, 
the  values  in  this  table  can  simply  be  used  to  predict  the  subsequent  modes  of  the  structure.  Note 
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Figure  24:  FFT  and  Coherence  Plot  for  7-Patch  Actuated  Case 

that  the  first  theoretical  frequency  must  be  used  to  find  the  constant  multiplier  and  therefore  is  set 
equal  to  the  frequency  of  the  first  average  mode.  Equation  35  is  an  illustration  of  how  this  method 
can  be  used  to  predict  the  m  =  1,  n  =  1  natural  frequency. 


^ mn 


Kfin 


LOn  =  I<pu 


(35) 


J.6.2  Input  Signal.  The  input  signal  typically  used  for  vibrations  testing  is  random  noise. 
This  is  a  signal  that  contains  the  same  power  at  every  frequency.  Since  the  natural  frequencies  of 
the  optical  structures  was  predicted  to  be  in  the  hundreds  of  Hertz,  exciting  the  entire  frequency 
spectrum  could  prove  inefficient.  For  this  reason,  both  a  pseudo  random  noise  signal  and  a  periodic 
chirp  signal  were  used  to  excite  the  test  articles. 
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A  pseudo  random  noise  signal  has  equal  power  only  over  a  specified  frequency  range.  A  periodic 
chirp  signal  emits  sinusoidal  oscillations  over  a  specified  frequency  range.  The  chirp  signal  begins 
with  the  low  end  of  the  frequency  range  and  ends  at  the  high  end. 

In  order  to  test  the  dependence  of  the  vibration  data  on  the  type  of  input  signal,  these  signals 
were  run  over  the  frequency  range  from  50  to  1400  Hz. 


m  -65 
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Figure  25:  FFT  and  Coherence  Plot  for  Pseudo  Random  Noise  Input  Case 

In  total,  the  coherence  plots  for  both  test  cases  are  relatively  good,  averaging  values  above  80%. 
The  noise  generated  FFT  of  Figure  25  compares  closely  with  that  of  the  chirp  generated  FFT  in 
Figure  26  with  respect  to  location  and  magnitude  of  the  first  mode.  At  higher  frequencies,  the  noise 
generated  FFT  suffers  a  loss  of  detail  and  the  peaks  are  not  as  pronounced.  The  chirp  generated 
FFT  has  the  most  detailed  information  about  these  primary  modes. 

Due  to  the  clean  data  pattern  relative  to  a  pseudo  random  noise  input  signal,  the  periodic 
chirp  input  signal  is  determined  to  be  slightly  better  for  acquisition.  For  this  reason,  all  further  tests 
use  the  chirp  signal  generator. 


4-6.3  Orientation.  In  order  to  test  the  orientation  dependence  of  the  membrane  test 
article,  tests  were  performed  in  the  ‘standard’  orientation  and  then  the  specimen  was  rotated  180° 
and  tested  again.  In  order  to  protect  the  coated  mirrors,  they  were  tested  in  the  upside-down 
position.  Therefore,  the  test  membrane  was  also  flipped  upside-down  and  tested. 

The  ‘standard’  position  was  chosen  by  convenience  and  is  arbitrary.  For  this  discussion  it  will 
be  defined  as  the  control  lead  for  the  center  patch  running  parallel  to  the  iron  bar  supports  and 
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Figure  26:  FFT  and  Coherence  Plot  for  Periodic  Chirp  Input  Case 
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Figure  27:  Definition  of  ‘Standard’  Orientation 
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toward  the  nearest  table  edge.  This  orientation  is  clearly  illustrated  in  Figure  27.  Notice  the  PVDF 
flap  where  the  center  patch  control  lead  runs  and  the  grounding  wire  (covered  by  electrical  tape) 
can  be  seen  on  the  right  of  the  membrane. 


Magnitude  (dB)  Plot 


x-dir 


Figure  28:  Magnitude  and  Phase  for  ‘Standard’  Orientation  Case  (Freq=467Hz,  m=2,  n=l) 

Figure  28  shows  the  m=l,  n=2  mode  for  the  standard  orientation.  It  is  clear  from  Figure  29  that 
the  rotation  changed  the  vibration  of  the  membrane.  The  darkened  lines  indicating  zero  crossings 
on  the  membrane  surface  have  rotated  in  accordance  with  the  membrane  support  structure,  as  is 
expected. 

After  the  membrane  is  flipped  up-side  down  (Figure  30)  the  zero  crossing  lines  on  the  magnitude 
plot  return  to  the  orientation  of  the  original  standard  position.  This  is  also  as  expected  since  the 
membrane  is  now  in  an  upside-down,  180°  rotated  orientation. 

The  dark  bands  on  the  phase  plots,  which  indicate  180°  out  of  phase  are  also  important  to 
notice.  The  right  side-up  test  articles  show  that  the  center  portion  is  out  of  phase,  while  the  upside- 
down  membrane  shows  the  sides  out  of  phase,  as  expected. 

The  mean  for  the  first  and  second  natural  frequencies  were  calculated  and  compared  to  the 
theoretical  value  (obtained  as  in  Table  9).  This  information  is  given  in  Table  10. 


Table  10:  Comparison  of  Standard,  Rotated  and  Flipped  Frequencies  Data  with  Circular  Mem¬ 

brane  Theory 


Theory 

Avg  Freq 

Standard  Freq 

Rotated  Freq 

Flipped  Freq 

m=0,  n=l 

191.3  Hz 

191.3  Hz 

191  Hz 

192  Hz 

191  Hz 

m=2,  n=l 

410  Hz 

476.3  Hz 

467  Hz 

478  Hz 

484  Hz 
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Figure  29:  Magnitude  and  Phase  for  Rotated  Case  (Freq=478Hz,  m=2,  n=l) 
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Figure  30:  Magnitude  and  Phase  for  Upside-down  Case  (Freq=484Hz,  m=2,  n=l) 
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Since  the  results  of  these  tests  show  almost  no  variation  in  the  fundamental  frequencies  of  the 
test  article,  the  orientation  of  the  test  article  is  shown  to  be  largely  arbitrary  for  a  given  test.  There¬ 
fore,  testing  the  coated  membrane  optical  structures  in  an  inverted  position  to  protect  the  optical 
surface  are  assumed  to  yield  very  similar  results  to  tests  run  in  the  ‘standard’  test  configuration. 

4-6-4  Membrane  Coating.  Each  membrane  must  be  coated  with  silicone  to  improve  the 
optical  quality  of  the  surface.  This  coating  drastically  changes  the  dynamics  of  the  membrane.  For 
this  reason,  a  comparison  between  an  uncoated  mirror  and  a  coated  mirror  is  justified.  The  7-patch 
mirror  was  used  for  the  comparison. 


Figure  31:  FFT  and  Coherence  Plot  for  Uncoated  Case 

Figure  31  shows  the  Fast  Fourier  transform  (FFT)  and  coherence  data  for  the  uncoated  case. 
Figure  32  gives  the  same  information  for  the  silicone  coated  membrane.  Notice  that  the  high  fre¬ 
quency  vibrations  (above  400  Hz  or  so)  are  severely  damped  out  by  the  silicone  coated  membrane. 

The  first  three  peaks  in  both  the  coated  and  uncoated  cases  represent  the  same  modes  of 
vibration.  These  frequencies  can  be  compared  to  the  theoretical  prediction  of  the  fundamental 
frequencies  for  these  modes.  This  comparison  will  enable  an  estimation  of  how  the  presence  of  the 
silicone  substrate  affects  the  dynamics  of  the  membrane.  Table  11  provides  a  comparison  of  the  first 
three  modes  for  the  uncoated  case  (obtained  as  in  Table  9).  The  first  three  peaks  for  the  uncoated 
case  a 
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Table  11: 


Table  12: 


Figure  32:  FFT  and  Coherence  Plot  for  Coated  Case 


Uncoatecl  Comparison  of  Membrane  Frequencies  with  Circular  Membrane  Theory 


Theory 

Experimental 

%  Difference 

m=0,  n=l 

218  Hz 

218  Hz 

0.0% 

m=2,  n=l 

465.5  Hz 

547  Hz 

17.5% 

m=0,  n=2 

500.4  Hz 

656  Hz 

31.1% 

Coated  Comparison  of  Membrane  Frequencies  with  Circular  Membrane  Theory 


Theory 

Experimental 

%  Difference 

m=0,  n=l 

64.5  Hz 

76  Hz 

17.8% 

m=2,  n=l 

137.85  Hz 

194  Hz 

40.0% 

m=0,  n=2 

148.1  Hz 

228  Hz 

53.8% 
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Table  12  gives  the  experimentally  observed  fundamental  frequencies  for  the  coated  case.  The 
theoretical  prediction  is  obtained  by  the  use  of  Equation  29,  where  tension,  mass  and  the  zeros  of 
the  Bessel  function  are  known  quantities,  addressed  in  Table  8,  Table  4  and  Table  6,  respectively. 

The  results  indicate  that  a  frequency  shift  of  approximately  142  Hz  occurs  after  the  surface  is 
coated  with  silicone.  Theory  predicts  that  a  shift  of  153.5  Hz  should  occur  due  to  the  added  mass  on 
the  membrane.  There  is  approximately  11.5%  error,  which  is  well  within  the  error  deviations  that 
have  been  recorded  for  these  tests.  This  indicates  that  even  with  the  silicone  substrate,  the  system 
is  still  modeled  well  by  the  membrane  equations,  as  long  as  the  additional  mass  is  acounted  for  in 
the  system. 

4-7  Summary 

In  this  chapter,  a  modal  analysis  of  the  test  articles  has  been  carried  out  and  the  data  has 
provided  information  regarding  the  dynamics  of  the  optical  structure.  The  theoretical  development 
was  relied  upon  for  the  prediction  and  calculation  of  fundamental  frequencies  and  mode  shapes. 
These  theoretical  values  were  compared  to  the  experimentally  observed  values.  The  modal  analysis 
data  was  also  utilized  to  calculate  the  tension  in  the  membrane. 

The  tension  in  the  membrane  was  determined  to  be  3.0  to  4.6  —  with  a  standard  deviation 

m 

of  0.5  to  1.0  ^ .  The  orientation  of  the  test  article  was  determined  to  have  little  to  no  impact  on 
the  frequency  analysis  data.  The  silicone  coated  test  membrane  was  determined  to  match  the  values 
obtained  by  pure  membrane  theory  to  about  11.5%,  which  indicates  that  the  coated  membrane  can 
still  be  well  approximated  by  membrane  theory. 

The  insight  gained  through  this  modal  analysis  is  useful  to  the  development  of  understanding 
about  thin  deformable  structures  and  numerical  researchers  studying  this  membrane  system.  An 
improved  data  acquisition  system  would  allow  some  of  this  insight  to  be  applied  to  the  control 
problems  associated  with  this  optical  system. 
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V.  Wavefront  Imaging  of  Membrane  Surface 

5. 1  Overview 

This  chapter  provides  an  explanation  of  the  testing  procedure  and  results  for  the  wavefront 
imaging  experiments.  Included  sections  describe  the  test  setup,  data  collection  and  data  processing 
tools  that  are  developed  and  utilized.  Static  deformation  results  are  given  a  statistical  treatment 
and  compared  to  previous  research.  Dynamic  measurements  are  also  made  of  the  surface  and  the 
results  are  used  to  characterize  the  data  acquisition  system. 

5.2  Surface  Deflection  in  Terms  of  Zernike  Polynomials 

The  classical  formulation  for  the  expression  of  the  displacement  of  the  surface  in  terms  of  Bessel 
functions  (as  explained  in  Section  4.2.3)  is  not  the  only  way  to  mathematically  express  the  surface 
deformation.  Another  orthogonal  set  of  polynomials  which  are  often  used  to  express  the  surface  are 
called  Zernike  Polynomials.  This  set  is  often  used  in  the  optics  community. 

Appendix  C.3  gives  a  list  of  the  first  42  Zernike  Polynomials.  Since  the  surface  is  described 
by  a  linear  combination  of  the  polynomials,  additional  terms  in  the  sum  increases  the  accuracy  of 
the  surface  match  to  the  data.  42  polynomials  were  chosen  in  order  to  capture  the  most  significant 
deflection  of  the  surface. 

The  Zernike  polynomials  are  functions  of  radius,  r,  and  angle,  9.  Given  a  polar  coordinate, 
Equation  36  will  return  a  surface  deflection.  The  data  taken  during  the  wavefront  experiments  was 
Zernike  coefficients,  this  data  is  scaled  by  the  WaveScope®  software  to  give  a  solution  in  terms  of 
microns  of  surface  deflection  (see  Section  5.4.1). 


n 

W(r,0)=^2<:izi(r,0)  (36) 

»= i 

Here  £  is  the  Zernike  coefficient  and  Z  is  the  Zernike  polynomial.  Figure  33  provides  a  surface 
plot  of  the  first  25  Zernike  polynomials  for  comparison  with  Figure  17.  These  plots  illustrate  the 
Zflr,  6)  polynomials,  which  are  eventually  scaled  by  the  wavefront  data,  The  x,  y  and  z-axis  have 
all  been  scaled  from  -1  to  1.  The  mathematics  required  to  obtain  a  mapping  from  the  set  of  Bessel 
functions  to  the  set  of  Zernike  polynomials  is  beyond  the  scope  of  this  thesis.  For  the  purpose  of 
this  research,  data  will  be  transformed  into  terms  of  surface  displacements  before  analysis,  through 
the  use  of  one  of  the  methods  described  above. 
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Figure  33:  Component  Breakdown  of  Zernike  Polynomials 


5.3  Test  Setup 

The  test  setup  on  the  optical  table  is  an  extension  of  a  prior  arrangement  used  by  previous 
researchers  at  AFIT  [21].  For  completeness,  the  setup  of  the  optical  table  will  be  described  again 
here.  Any  significant  deviations  from  this  setup  will  also  be  noted. 

A  20  mWatt  helium-neon  laser  (A  =  633  mi)  was  used  to  illuminate  the  test  subjects.  The 
optical  path  of  the  laser  is  outlined  in  Figure  34  through  36. 

Light  from  the  laser  first  passes  through  a  set  of  neutral  density  filters,  which  alter  the  intensity 
of  the  beam  to  a  usable  level  for  the  wavefront  sensor.  The  filter  set  includes  a  primary  filter  wheel, 
with  lenses  ranging  from  10%  to  80%  transmission  and  two  variable  gradient  filter  wheels  which  can 
be  used  for  fine  adjustments.  These  filters  are  shown  in  Figure  35.  Due  to  differences  in  optical 
quality  and  reflectivity,  the  filter  wheels  are  tuned  for  every  test  article. 

The  optical  path  is  then  turned  90°  by  a  ^  flat  mirror  (Small  Pivot  Mirror  A,  34).  The  beam 
passes  through  a  spatial  filter,  which  produces  a  clean  Gaussian  wavefront  to  illuminate  the  test 
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Figure  34:  Optical  Table  Test:  Global  View 


article.  At  this  point,  the  laser  light  is  expanding.  A  1  inch  achromatic  doublet  lens  (Doublet  Lens 
A)  is  utilized  to  collimate  the  expanding  beam.  This  lens  has  a  focal  length  of  200  mm. 

A  Shack-Hartmann  type  wavefront  sensor,  requires  two  beams  for  wavefront  measurement. 
The  first  beam  is  the  image  from  the  optical  surface  under  investigation  and  will  be  referred  to  as 
the  test  beam.  The  second  beam  is  a  reference  beam  from  a  flat  optical  surface  and  will  be  referred 
to  as  the  reference  beam.  In  order  to  generate  these  two  beams  from  the  same  source,  a  beam  splitter 
is  utilized. 

The  collimated  laser  light  now  meets  a  3  inch  |  wedge  beam  splitter,  which  separates  the  beam 
into  two  equal  intensity  beams.  The  reference  beam  is  turned  90°,  reflected  off  of  a  ^  flat  mirror 
(Reference  Beam  Flat  Mirror,  36),  back  through  the  beam  splitter  and  into  the  wavefront  sensor. 

The  test  beam  is  allowed  to  pass  through  the  beam  splitter.  This  beam  is  focused  using  a  1 
inch  doublet  lens,  with  a  focal  length  of  250  mm.  The  focal  point  of  the  test  beam  is  located  ona-jj 
flat  mirror  (Small  Pivot  Mirror  B,  36),  which  turns  the  expanding  beam  toward  the  end  of  the  table. 
A  12.5  inch  parabolic  reflector,  with  a  focal  length  of  75.125  inches,  is  used  to  stop  the  expansion  of 
the  beam  and  to  direct  the  beam  path  onto  the  test  article.  Finally,  a  6  inch  flat  mirror,  mounted 
at  45°  from  vertical,  was  used  to  turn  the  test  beam  90°  downward  onto  the  test  article. 

The  test  beam  is  reflected  off  of  the  test  article  and  allowed  to  follow  the  same  optical  path 
back  to  the  beam  splitter.  The  beam  splitter  turns  the  test  beam  90°  into  the  wavefront  sensor. 
The  test  article  is  suspended  over  the  optics  table  by  spacers,  which  allows  the  membrane  to  vibrate 
freely. 
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Figure  35:  Optical  Table  Test:  Filter  View 


5-4  Data  Collection 

Data  collection  from  the  wavefront  sensor  can  be  completed  with  the  WaveScope®  electronics 
box  and  acquisition  software,  however  for  situations  where  greater  versatility  is  needed  a  data  acquisi¬ 
tion  methodology  was  developed  to  use  Matlab®  running  on  a  separate  computer.  This  methodology 
lays  the  foundation  for  future  development  of  the  system  for  use  with  Simulink®  controls  software 
and  dSpace®  control  system  prototyping  hardware.  This  section  provides  an  outline  of  the  new  data 
acquisition  setup.  The  system  is  illustrated  in  the  form  of  a  flow  chart  in  Figure  37. 

All  wavefront  measurements  utilized  a  Shack-Hartmann  type,  WaveScope®  wavefront  sensor, 
manufactured  by  Adaptive  Optics  Associates.  Data  is  generated  from  the  test  and  reference  beams 
(see  Section  5.3)  and  processed  by  the  WaveScope®  electronics  box  (Figure  38). 

The  data  is  then  used  to  calculate  various  coefficients  and  images  by  the  WaveScope®  software. 
For  the  present  system,  it  is  desirable  to  disable  all  unnecessary  calculations  and  image  displays  on 
the  WaveScope®  computer.  This  will  save  time  and  allow  the  data  to  stream  as  quickly  as  possible 
to  the  Matlab®  acquisition  computer. 

Software  engineers  at  Adaptive  Optics  Associates  graciously  provided  a  script  to  enable  Zernike 
coefficients  generated  by  the  WaveScope®  computer  to  be  exported  via  an  ethernet  connection.  This 
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Figure  36:  Optical  Table  Test:  Test  Specimen  View 

data  stream  from  the  WaveScope®  computer’s  ethernet  port  was  converted  into  a  serial  data  stream 
using  the  Nport®  Serial  to  ethernet  conversion  box. 

This  was  a  necessary  step  because  the  dSpace®  hardware  does  not  currently  support  data 
acquisition  through  an  ethernet  connection.  In  order  to  insure  future  compatibility  with  the  dSpace 
system,  the  conversion  to  the  serial  data  was  made. 

The  input  signal  for  all  tests  was  generated  by  an  Hewlett-Packard  1.5  MHz  function  generator 
and  amplified  by  a  stack  of  4  Trek  PZD  Duel  Channel  Amplifier  (see  Figure  39).  Transfer  functions 
for  these  amplifiers  are  given  in  Appendix  C.l.  The  output  of  the  amplifier  is  verified  using  an  Agilent 
6.5  digit  multimeter  before  being  used  to  actuate  the  control  patches.  Data  from  the  wavefront  is 
received  and  processed  using  the  acquisition  and  control  hardware  pictured  in  Figure  40. 

5-4-1  Data  Processing.  The  WaveScope®  software  allows  the  user  to  choose  the  number  of 
Zernike  coefficients  to  be  calculated,  with  a  maximum  of  42  values.  To  obtain  the  most  information 
possible  with  the  given  equipment,  42  Zernike  coefficients  are  exported  at  a  rate  of  approximately 
10  Hz.  No  averaging  is  performed  by  the  WaveScope®  software. 

Once  the  WaveScope®  software  has  been  configured  to  export  the  Zernike  Coefficients  (see 
Appendix  A. 6),  the  Seriallnit.m  (see  Appendix  B.ll)  script  file  is  run  on  the  Matlab®  computer. 
This  script  initializes  the  communication  to  the  serial  port  and  defines  some  characteristics  of  the 
incoming  data  stream. 
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Figure  37:  Flow  Chart  of  Wavescope  Experiment 


The  basic  command  to  grab  a  set  of  coefficients  from  the  serial  port  is  given  in  the  script 
SerialGrab.m  (see  Appendix  B.10).  This  script  simply  looks  for  the  start  of  a  set  of  data  and 
takes  42  lines  after  that  from  the  buffer.  A  provisional  output  has  been  included  to  tell  the  user 
if  noncontiguous  frames  have  been  taken.  This  lets  the  user  know  if  sets  of  Zernike  coefficients  are 
being  corrupted  or  overwritten  due  to  an  improper  communication  setting. 

The  data  processing  is  done  with  the  SerialPlot.m  script  (see  Appendix  B.12).  This  script 
serves  as  the  user  interface  to  the  data  acquisition  system.  The  program  first  calls  a  function  (see 
Appendix  B.18  to  calculate  a  surface  mesh  of  Zernike  Polynomials  (see  Appendix  C.3).  The  program 
then  asks  the  user  to  shut  off  all  inputs  to  the  system  so  that  a  static  reference  frame  of  the  mirror 
surface  can  be  taken.  This  reference  frame  is  an  average  of  30  sets  of  Zernike  polynomials  and  is 
subsequently  subtracted  from  all  measurements  of  the  surface. 

Next  the  SerialPlot.m  will  prompt  the  user  to  turn  on  equipment  again,  if  needed.  A  specified 
number  of  Zernike  coefficient,  n,  sets  are  then  read  from  the  serial  port.  Each  set  used  to  obtain  the 
deviation  of  the  surface  from  the  static  reference  frame  with  Equation  37. 
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Figure  38:  Electronic  Control  Boxes 
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Figure  39:  Signal  Generation  and  Data  Processing  Hardware 


n 

W{r,e)  =  Y,CiZi{r,e)  (37) 

i=l 

The  set  of  n  surface  shapes  is  then  stored  in  a  *.MAT  file  for  later  analysis  and  a  movie  of  the 
deformation  is  created.  Additional  post-processing  is  accomplished  using  the  StaticPlotter.m  and 
StaticFramePlot.m  script  files  (see  Appendix  B.14  and  B.15).  The  data  manipulation  used  for  these 
files  will  be  explained  as  the  results  are  discussed  in  the  remaining  sections  in  this  chapter. 
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Figure  40:  Data  Acquisition  and  Control  Hardware 

5.5  Data  Analysis:  Static  Imaging 

The  behavior  of  the  membrane-like  optical  structure  under  static  voltage  can  be  investigated 
using  the  dynamic  sampling  acquisition  system  to  generate  enough  data  for  statistical  analysis.  Mean 
Static  deformation  plots  will  be  discussed  for  each  of  the  patches.  A  calculation  of  the  standard 
deviation  in  the  data  is  provided.  Finally,  the  dependence  of  the  deformation  with  respect  to  applied 
voltage  is  discussed. 

5.5.1  Mean  Surface  Deformation.  The  mean  deformation  of  the  surface  with  respect  to 
the  static  reference  is  an  important  measurement.  This  information  is  obtained  by  averaging  the 
surface  deformation  data  calculated  with  the  SerialPlot.m  script.  The  code  which  was  developed  to 
perform  this  analysis  is  called  StaticPlotter.m  (see  Appendix  B.14). 

This  script  file  first  prompts  the  user  for  the  *.MAT  file  to  be  loaded.  It  discards  the  first  and 
last  surface  profiles  and  takes  the  average  and  standard  deviation  of  the  remaining  sets.  For  these 
experiments  30  sets  of  data  were  taken,  therefore  28  sets  of  data  are  averaged  to  obtain  the  mean 
surface  deformation. 

The  laser  could  not  be  expanded  enough  to  sample  the  entire  surface  of  the  mirror,  therefore 
area  that  could  be  reliably  illuminated  was  centered  on  the  middle  of  the  membrane  and  sampled  with 
the  wavefront  sensor.  In  order  to  match  the  scale  of  the  mirror  and  give  an  accurate  representation 
of  the  area  where  data  was  acquired  from  the  optical  surface,  the  sampled  data  is  padded  with 
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zeros  around  the  outside.  The  zero  displacement  region  in  the  images  represents  the  area  where  no 
wavefront  measurements  were  taken,  however,  the  surface  is  surely  deformed  in  these  regions. 

The  control  patch  design  was  then  overlaid  on  the  surface  deformation  plots.  This  information 
is  calculated  based  on  the  scale  of  the  original  computer  design  for  the  control  patch  and  is  not 
simply  a  drawing  of  what  appears  to  fit  the  data.  The  data  set  is  rotated  clockwise  until  a  close 
match  with  the  overlaid  control  pattern  is  obtained.  Typically  the  data  needed  rotation  of  about 
15°  to  match  the  control  pattern.  This  rotation  is  simply  a  matter  of  how  the  mirror  was  oriented 
on  the  test  stand  during  testing  and  has  no  effect  on  the  actual  data. 

The  maximum  and  minimum  deformation  is  calculated  over  the  entire  surface  and  this  infor¬ 
mation  is  provided  in  Table  13. 


Table  13:  Maximum  and  Minimum  Surface  Displacements 


Patch 

Min  Disp  ( fim) 

Max  Disp  ( /mi) 

PV  Value  ( fim) 

PV  Value  (A) 

1 

-1.0121 

0.79105 

1.80314 

2.84857 

2 

-0.70536 

0.67029 

1.37565 

2.17322 

3 

-1.4521 

0.98718 

2.43926 

3.85349 

4 

-1.6321 

0.98338 

2.61553 

4.13196 

5 

-1.7879 

0.83931 

2.62723 

4.15044 

6 

-2.0969 

1.0568 

3.15369 

4.98214 

7 

-1.1939 

0.73168 

1.92555 

3.04194 

Peak  to  valley  (PV)  calculations  can  be  calculated  for  each  data  sample  by  simply  adding  the 
absolute  value  of  the  minimum  deflection  to  the  maximum  deflection.  This  information  is  useful  for 
comparison  and  gives  a  good  idea  of  the  total  variation  of  the  mirror  surface.  The  PV  value  can 
also  be  used  to  describe  the  surface  quality  of  the  membrane-like  structure.  This  information  is  also 
provided  in  Table  13. 

Figure  41  is  the  ‘static  reference  frame’  which  was  subtracted  from  every  measurement  which 
was  taken.  This  is  due  to  the  surface  roughness  present  on  the  membrane’s  silicone  coating.  The 
frame  is  taken  prior  to  every  measurement  and  simply  subtracted  from  the  measurement  to  remove 
the  bias  due  to  this  roughness.  The  PV  value  of  the  resting  surface  is  determined  to  be  1.7  /r?n. 

The  actuated  center  control  patch  did  not  show  very  good  agreement  with  the  predicted  dis¬ 
placement  area.  The  actuated  area  appeared  between  control  patches  4  and  5,  with  maximum 
deflection  of  -1.01  microns.  One  possible  cause  of  this  displacement  is  that  the  test  beam  is  not 
centered  on  the  mirror  surface.  This,  however,  is  unlikely  to  be  the  case  because  the  centering  was 
verified  with  measurement  equipment  at  the  time  of  testing  and  the  other  control  patch  actuation 
results  all  appear  centered  on  the  surface.  This  patch  may  have  been  shorting  to  one  of  the  other 
patches  through  some  unforeseen  mechanism  in  the  test  setup. 
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Figure  41: 


Mean  Surface  Deformation  of  Zero  Voltage  Reference 


Figure  42:  Mean  Surface  Deformation  with  Patch  1,  Actuation  at  0  Hz,  400  V 
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Figure  43:  Mean  Surface  Deformation  with  Patch  2,  Actuation  at  0  Hz,  400  V 


Figure  44: 


Mean  Surface  Deformation  with  Patch  3,  Actuation  at  0  Hz,  400  V 


62 


Figure  45: 


Figure  46: 


Mean  Surface  Deformation  with  Patch  4,  Actuation  at  0  Hz,  400  V 


Mean  Surface  Deformation  with  Patch  5,  Actuation  at  0  Hz,  400  V 
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Figure  47: 


Figure  48: 


Mean  Surface  Deformation  with  Patch  6,  Actuation  at  0  Hz,  400  V 


Mean  Surface  Deformation  with  Patch  7,  Actuation  at  0  Hz,  400  V 
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Figures  42  through  48  represent  the  surface  deformation  patterns  for  patches  1  through  7, 
respectively.  Each  patch  is  actuated  at  a  static  voltage  of  400  V.  All  exterior  control  patch  actuation 
results  appear  to  fall  in  the  predicted  region.  The  low  area  in  each  plot  falls  within  the  lines  of  the 
control  patch  overlay,  which  is  expected.  The  displacements  in  the  membranes  range  from  to  -0.71 
to  -2.1  microns. 

5.5.2  RadialCuts.  Since  the  out  of  plane  deformation  is  most  interesting  in  the  area  of 
each  patch,  a  diagram  has  been  designed  to  show  this  deformation  for  a  radial  cut  pattern  around 
the  mirror.  There  are  6  radial  cuts  made,  one  through  each  of  the  outer  control  patches.  Each  cut  is 
located  on  the  line  connecting  the  of  the  control  patch  and  the  center  of  the  mirror.  By  examining 
these  patterns,  details  regarding  the  relative  displacements  of  the  surface  when  a  given  patch  is 
actuated  are  brought  out.  The  post-processed  data  used  in  Section  5.4.1  is  used  to  create  these  plots 
and  no  additional  processing  was  required.  The  data  points  represent  the  mean  deformation  from 
the  static  reference  frame  and  the  error  bars  are  figured  by  a  calculation  of  the  standard  deviation 
from  that  mean. 
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Figure  49:  Radial  Cut  of  Patch  2,  Actuation  at  0  Hz,  400  V 

Figures  49  through  54  represent  the  radial  cut  patterns  for  patches  2  through  7,  respectively. 
These  plots  are  simply  glimpses  of  the  data  contained  in  the  surface  deformation  plots  and  therefore 
exhibit  similar  statistical  characteristics. 

These  figures  highlight  a  unique  behavior.  The  actuated  control  patch  has  a  negative  deflection, 
as  expected,  but  the  neighboring  control  patches  often  exhibit  a  positive  compensative  deflection.  A 
great  example  is  the  plot  of  the  actuation  of  control  patch  3.  This  surface  exhibited  a  large  positive 
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Figure  50:  Radial  Cut  of  Patch  3,  Actuation  at  0  Hz,  400  V 
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Figure  51:  Radial  Cut  of  Patch  4,  Actuation  at  0  Hz,  400  V 
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Figure  52:  Radial  Cut  of  Patch  5,  Actuation  at  0  Hz,  400  V 
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Figure  53:  Radial  Cut  of  Patch  6,  Actuation  at  0  Hz,  400  V 
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Figure  54:  Radial  Cut  of  Patch  7,  Actuation  at  0  Hz,  400  V 

deformation  in  the  region  of  control  patch  4.  This  behavior  is  most  likely  tied  to  the  fact  that  the 
membrane  coating  cannot  return  to  zero  deformation  very  quickly  in  the  interface  regions  of  the 
membrane  (in  fact  the  reason  it  is  on  the  surface  at  all  is  for  the  smoothing  effect). 

5.5.3  Standard  Deviation.  The  standard  deviation  from  the  mean  surface  deflection  is 
calculated  using  the  Equation  38. 
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This  information  gives  an  idea  of  what  kind  of  variance  can  be  expected  in  the  measurement. 
A  low  variance  translates  to  a  higher  certainty  that  deformed  surface  will  be  close  to  the  stated 
mean. 


Actuation  on  the  center,  or  right  side  of  the  membrane  (Figures  55  to  58  showed  a  maximum 
standard  deviation  toward  the  outer  edge  of  the  right  hemisphere.  The  cases  where  the  actuation 
was  applied  to  the  left  patches  (Figures  59  to  61)  show  a  shift  in  the  maximum  standard  deviation 
to  the  center  section  of  the  membrane.  These  shifts  may  be  caused  by  slightly  uneven  tensioning  or 
non-uniformities  in  the  optical  system.  This  data  indicates  a  combination  of  two  results.  First,  the 
area  is  actuated  seems  to  have  an  increased  standard  deviation.  Second,  the  right  hemisphere  of  this 
mirror  seems  to  have  an  increased  variance.  The  increased  standard  deviation  for  the  patch  5  to  7 
in  the  center  region  may  be  a  result  of  a  combination  of  these  two  effects. 
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Figure  55:  Standard  Deviation  from  Mean  Deformation  of  Patch  1,  Actuation  at  0  Hz,  400  V 


Figure  56:  Standard  Deviation  from  Mean  Deformation  of  Patch  2,  Actuation  at  0  Hz,  400  V 


Table  14:  Maximum  Standard  Deviation  from  Mean  Surface  Deflection 


Patch 

Max  Standard  Deviation  from  Mean  (jirn) 

1 

0.47799 

2 

0.15178 

3 

0.1921 

4 

0.39651 

5 

0.27748 

6 

0.40928 

7 

0.38451 
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Figure  57:  Standard  Deviation  from  Mean  Deformation  of  Patch  3,  Actuation  at  0  Hz,  400  V 
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Figure  58:  Standard  Deviation  from  Mean  Deformation  of  Patch  4,  Actuation  at  0  Hz,  400  V 
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Figure  59:  Standard  Deviation  from  Mean  Deformation  of  Patch  5,  Actuation  at  0  Hz,  400  V 


'.5 
.45 
.4 
.35 
.3 
.25 
.2 
.15 
.1 

.05 

-BO  -40  -20  0  20  40  60 

X-dir  (mm) 

Figure  60:  Standard  Deviation  from  Mean  Deformation  of  Patch  6,  Actuation  at  0  Hz,  400  V 
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Figure  61:  Standard  Deviation  from  Mean  Deformation  of  Patch  7,  Actuation  at  0  Hz,  400  V 

The  standard  deviation  has  an  average  value  of  0.3271  pm.  The  maximum  values  for  the 
standard  deviation  of  each  plot  are  given  in  Table  14. 

5.5.4  Voltage  Dependence.  Figure  62  is  a  plot  of  the  observed  surface  deformation  as  the 
voltage  is  increased  from  100  to  600  volts.  Each  measurement  is  an  average  of  the  max  or  min  value 
of  the  seven  patches.  The  negative  values  are  the  displacement  of  the  surface  in  the  control  region. 
The  positive  values  represent  the  compensatory  effect  of  the  surrounding  region.  The  trend  is  very 
linear  in  the  region  from  200  to  600  V  and  an  increase  in  voltage  gives  a  larger  deflection. 

5.6  Data  Analysis:  Active  Imaging 

The  primary  objective  of  this  section  is  to  characterize  the  sampling  rate  of  the  new  data 
acquisition  system.  The  Zernike  coefficients  describe  the  surface  deformation  and  can  be  examined  in 
time  to  characterize  the  vibration  of  the  surface.  Significant  surface  deflection  were  observed  during 
static  actuation  for  patch  5,  therefore  this  patch  was  selected  for  this  analysis.  Figure  63  shows 
the  values  for  the  first  42  Zernike  coefficients.  This  plot  illustrates  the  fact  that  larger  coefficients 
typically  have  less  impact  (lower  values)  when  compared  to  the  small  coefficients.  This  is  due  to  the 
larger  coefficients  being  defined  to  describe  more  and  more  complex  surface  deformation  patterns. 

The  4th  Zernike  coefficient  is  chosen  for  tracking  due  to  its  large  magnitude.  By  keeping  the 
relative  size  of  the  tracked  coefficient  large,  compared  to  the  measurement  error  and  background 
noise,  the  results  are  isolated  from  those  errors. 
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Figure  62: 


Average  Surface  Deformation  as  a  Function  of  Voltage 


Figure  63:  Bar  Chart  of  all  42  Zernike  Coefficients  for  Patch  5  at  400  V 
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Figure  64:  Time  Signature  of  4th  Zernike  Coefficient  with  |  Hz  Input  Signal 


Figure  65:  Time  Signature  of  4th  Zernike  Coefficient  with  1  Hz  Input  Signal 
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Figure  66:  Time  Signature  of  4th  Zernike  Coefficient  with  2  Hz  Input  Signal 

Figure  64,  Figure  65  and  Figure  66  show  the  experimental  data  as  a  function  of  time  for  the 
4th  Zernike  coefficient.  This  data  was  obtained  by  taking  a  time  signature  every  time  the  data 
acquisition  system  read  a  line  from  the  serial  port  data  stream.  This  signature  is  then  processed  to 
obtain  a  time  span  for  the  data  read,  which  is  a  measurement  of  the  sampling  rate  of  the  system. 
The  Matlab®  code  which  was  used  to  make  these  calculation  is  provided  in  Appendix  B.3. 

The  input  signal  is  also  plotted  on  Figures  64  to  66.  Since  no  knowledge  of  where  the  data  is 
in  the  sinusoidal  cycle  can  be  easily  obtained,  a  fit  of  the  sinusoidal  input  curve  is  made  based  on 
an  L2  norm  from  a  difference  with  the  data  sample.  This  effectively  finds  the  best  fit  of  a  sinusoidal 
curve  to  the  provided  data.  As  illustrated  by  the  figures,  all  data  matches  the  input  signals  very 
closely. 

The  \  Hz  and  1  Hz  signals  were  slow  enough  that  the  data  acquisition  system  could  sample 
then  more  than  one  time  per  cycle.  In  the  case  of  the  2  Hz  input  signal,  the  data  acquisition  is  not 
able  to  sample  fast  enough  to  hit  every  cycle,  which  is  why  the  oscillatory  pattern  in  the  data  is 
apparent.  This  behavior  is  called  aliasing  and  is  associated  with  any  data  acquisition  system. 

The  Nyquist  frequency  is  said  to  be  half  of  the  sample  frequency  and  provides  a  line  below 
which  aliasing  does  not  occur.  The  fact  that  a  2  Hz  input  is  exhibiting  this  aliasing  phenomenon 
indicates  that  the  Nyquist  frequency  for  our  system  is  slightly  less  than  2  Hz.  Therefore,  the  sample 
frequency  of  the  data  acquisition  system  is  approximately  4  Hz. 

The  sample  frequency  was  also  calculated  analytically  based  on  the  time  signatures  which  were 
obtained  during  testing.  The  average  sample  rate  was  determined  to  be  4.2406  Hz,  with  a  standard 
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deviation  of  0.029414  Hz.  This  sample  rate  is  slightly  higher  than  that  predicted  by  the  graphical 
method. 

The  system  acquisition  speed  was  also  tested  based  on  workload.  The  number  of  Zernike 
coefficients  calculated  by  the  WaveScope®  computer  and  processed  by  the  Matlab®  were  reduced 
to  investigate  is  a  sample  rate  speed  increase  was  possible  in  the  system. 
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Figure  67:  Sample  Rate  as  a  function  of  Zernike  Coefficients  Computed 

Figure  67  gives  the  relationship  between  the  Zernike  coefficients  sampled  and  the  sample  rate 
of  the  system.  For  convenience,  a  second  order  trend  line  has  been  fit  to  the  data  and  the  equation 
for  this  line  is  given  in  the  upper  right  corner  of  the  figure.  This  figure  illustrates  the  fact  that 
as  fewer  coefficients  are  computed,  the  sample  speed  in  increased.  A  maximum  sample  speed  of 
approximately  14  Hz  can  be  obtained,  however,  so  few  Zernike  coefficients  would  be  available  that 
the  surface  deformation  prediction  could  not  accurately  represent  the  actual  surface. 


5. 1  Summary 

In  this  chapter,  the  new  data  acquisition  was  described  in  detail.  This  new  system  is  capable 
of  taking  data  at  sample  rates  as  high  as  14  Hz,  though  4  Hz  was  average  for  the  computation  of 
42  Zernike  coefficients.  The  acquisition  system  was  used  to  perform  a  statistical  measurement  of 
the  static  surface  deformations.  The  maximum  static  deformation  was  determined  to  be  -2.1  /im 
and  occurred  for  the  actuation  of  patch  6.  All  exterior  patches  appeared  to  actuate  in  the  region 
predicted  by  the  placement  of  the  control  patches.  The  center  patch  did  not  actuate  as  expected 
and  a  possible  short  circuit  in  the  control  lead  is  suspected  as  the  cause. 
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VI.  Conclusion 


6. 1  Overview  of  Experiments 

This  project  contained  two  distinct  experimental  test  setups.  A  scanning  laser  vibrometer  was 
used  to  characterize  the  optical  structures  over  a  wide  range  of  frequencies  and  to  calculate  the 
tension  in  the  membrane.  An  active  wavefront  sensor  was  used  to  characterize  the  optical  surface 
subject  to  static  and  periodic  applied  voltages.  In  addition  to  the  outline  of  these  two  experiments, 
an  overview  of  the  fabrication  and  design  process  is  provided. 

6.1.1  Fabrication  and  Design.  The  fabrication  and  design  process  included  the  design  of 
the  control  pattern  etched  onto  the  membrane,  the  selection  and  description  of  the  PVDF  actuator 
material,  mirror  etching  and  tensioning  and  coating  of  the  optical  surfaces.  The  areal  densities 
of  each  mirror  produced  is  also  provided  as  a  benchmark  for  the  fabrication  methods  used  in  this 
project. 

6.1.2  Modal  Analysis.  Next  a  modal  analysis  of  each  membrane  was  performed.  An 
extensive  theoretical  treatment  of  circular  membrane  theory  is  provided  along  with  an  overview  of 
the  test  equipment  and  setup,  data  collection  and  data  analysis  methods.  Frequency  response  curves 
and  surface  deformation  plots  are  provided  to  characterize  the  primary  modes  of  each  surface.  A 
estimation  of  the  tension  in  the  PVDF  membrane  is  provided.  Comparisons  are  also  made  based  on 
actuation  method,  input  signal,  orientation  and  membrane  coating. 

6.1.3  Wavefront  Analysis.  Finally  measurements  of  the  optical  wavefront  were  performed. 
A  description  of  the  test  setup  an  data  acquisition  system  is  provided  and  static  and  dynamic 
experimental  test  results  are  presented  and  discussed. 

6.2  Conclusions  Drawn 

6.2.1  Fabrication  and  Design.  While  an  established  fabrication  methodology  was  followed, 
some  improvements  to  the  technique  were  noted.  Thin  layers  of  silicone  are  not  necessarily  the 
best.  A  minimum  thickness  seems  to  exist  around  1.5  mm,  beyond  which  a  print  through  of  the 
control  pattern  could  be  observed  and  the  surface  was  not  considered  testable.  Exposing  the  curing 
membrane  to  vacuum  pressure  to  siphon  out  the  bubbles  is  not  a  valid  technique.  This  seemed 
to  create  a  myriad  of  surface  wrinkles  which  rendered  the  surface  untestable.  It  was  verified  that 
additional  layers  of  silicone  can  be  applied  directly  to  previous  gold  coated  silicon  layers,  with  no 
apparent  problem  of  wrinkling  or  bubble  formation  at  the  interface.  Finally,  structures  with  optical 
quality  reflecting  surfaces  can  routinely  be  produced  using  this  methodology.  These  structures  have 
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been  calculated  to  have  an  areal  density  of  1.5  -  4  Mr  and  thicknesses  of  1.5  to  3.3  mm,  which 
represents  a  50%  reduction  in  thichness  over  previous  research  completed  at  AFIT. 

6.2.2  Modal  Analysis.  The  tension  is  a  valuable  parameter  for  numerical  modeling  of  mem¬ 
brane  structures.  Modal  analysis  results  were  used  to  estimate  the  tension  in  the  PVDF  membrane. 
This  tension  was  estimated  to  be  3.0  to  4.6  —  with  a  standard  deviation  of  0.5  to  1.0  — . 

m  m 

Various  other  attributes  of  the  test  articles  were  also  investigated.  The  results  indicate  that 
the  silicone  coated  membrane  is  well  approximated  by  thin  membrane  theory.  Orientation  of  the 
test  article  on  the  test  stand  was  also  determined  to  have  little  to  no  impact  on  the  modal  analysis 
results. 

6.2.3  Wavefront  Analysis.  A  new  data  acquisition  system  was  developed  for  use  with 
AFIT’s  wavefront  imaging  hardware.  This  system  enables  static  data  to  be  taken  at  sample  rates 
up  to  14  Hz.  The  number  of  Zernike  coefficients  calculated  is  found  to  have  a  large  impact  on  the 
sample  rate  of  the  system.  42  Zernike  coefficients  give  a  reasonable  approximation  to  the  surface 
deformation  and  can  support  a  sample  rate  of  approximately  4  Hz. 

Static  deformations  at  400  volts  were  used  for  comparison.  This  actuation  method  produced 
surface  deformations  of  -0.71  to  -2.1  microns.  Each  exterior  control  patch  deformed  in  the  region 
which  was  expected.  The  redial  cut  diagrams  showed  that  often  times  when  a  particular  patch  is 
actuated,  its  neighboring  patches  are  adversely  affected.  A  calculated  average  standard  deviation 
was  determined  to  be  0.3271  microns.  The  largest  deviations  were  also  noted  to  be  in  the  patch 
which  was  being  actuated  or  toward  the  outside  of  the  observed  region  on  the  membrane.  Finally, 
the  surface  deformation  was  determined  to  go  nearly  linearly  with  increased  voltage  in  the  range  of 
200  to  600  volts. 

6.3  Areas  for  Further  Development 

It  is  the  author’s  hope  that  this  research  will  serve  as  an  appropriate  stepping  stone  for  the 
future.  There  are  many  possible  areas  which  seem  to  show  promise  for  possible  lines  of  research,  a 
few  are  discussed  below. 

6.3.1  Fabrication  and  Design.  The  silicone  substrate  which  was  utilized  for  this  experiment 
is  not  rated  for  use  in  space.  This  temperature  range  would  severely  change  the  dynamics  of  this 
coating  and  another  material  must  be  found  as  a  replacement.  Also  the  PVDF  etching  pattern  is 
time  consuming  and  relatively  inefficient  to  do  by  hand.  Other  methods  for  masking  the  etch  pattern 
onto  the  PVDF  material  would  allow  for  more  intricate  control  designs  and  speed  up  the  fabrication 
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process.  Finally,  better  methods  of  approximating  the  tension  in  the  membrane  could  be  developed. 
The  use  of  force  transducers  or  other  methods  could  be  used  to  verify  the  results  obtained  from 
modal  analysis. 

6.3.2  Wavefront  Analysis.  The  largest  area  of  possible  advancement  is  in  the  area  of  con¬ 
trol.  Now  that  a  way  of  collecting  data  in  ‘real  time’  has  been  demonstrated,  a  control  system  to 
damp  out  vibrations  on  the  surface  or  otherwise  shape  the  structure  are  within  reach.  Careful  atten¬ 
tion  has  been  paid  to  making  the  system  compatible  with  the  most  obvious  control  implementation 
tools,  namely  Simulink®  and  dSpace®. 


6  A  Summary 

This  project  consists  of  three  major  parts:  Fabrication  and  Design,  Modal  Analysis  and  Wave- 
front  Analysis.  The  relevant  information  has  been  presented  and  discussed  regarding  each  of  these 
three  areas.  Significant  advances  have  been  made  in  the  area  of  data  acquisition  and  processing 
using  AFIT’s  wavefront  analysis  hardware. 

The  development  of  this  methodology  has  brought  within  reach  the  seemingly  distant  goal 
of  ‘real  time’  investigation  of  the  wavefront  for  a  thin  deformable  membrane-like  structure.  This 
research  is  critical  to  the  overall  goal  of  placing  large  aperture  optical  structures  in  space.  Such 
systems  could  provide  valuable  information  for  the  mission  planners  and  war-fighters  of  tomorrow. 
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Appendix  A.  Laboratory  Notes 


A.l  Etching  the  Membrane 
A.  1.1  Supplies. 

•  Razor  blade,  sharpie,  straight  edge,  tape,  cotton  balls,  cutting  board,  paper  towel 

•  computer  generated  membrane  template 

•  rubber  safety  gloves 

•  Isopropyl  Alcohol  (91%  by  volume) 

•  Ferric  Chloride  Etchant  Solution 

•  14  by  14  inch  piece  of  PVDF  material  (52 pm,  with  Cu  Ni  coating) 

A.  1.2  Procedures  and  Comments. 

1.  Mark  the  positive  side  of  the  PVDF  with  a  (+)  in  the  corner  of  the  piece  with  the  black  marker. 
Packaging  of  PVDF  indicates  that  ‘positive  inside’  when  on  the  roll. 

2.  Measure  and  cut  a  14  inch  piece  from  the  roll  with  the  straight  edge  and  knife. 

3.  Use  the  knife  to  cut  the  computer  generated  membrane  template  out  as  cleanly  as  possible. 
The  area  that  is  to  be  masked  with  the  permanent  marker  must  be  completely  removed. 

4.  Mark  the  orientation  of  the  PVDF  sheet  with  3  parallel  lines,  parallel  to  the  direction  that  the 
PVDF  was  rolled  up. 

5.  Tape  down  the  template  in  the  appropriate  position. 

6.  Mask  the  PVDF  by  tracing  the  template  and  shading  in  with  permanent  marker. 

7.  Remove  the  template. 

8.  Carefully  clean  up  the  edges  of  the  mask  with  the  straight  edge  and  marker.  Note:  It  is  easy 
to  smudge  the  marker  at  the  point,  be  sure  to  remove  unwanted  marks  with  alcohol  before 
etching. 

9.  Put  on  the  rubber  gloves  and  ventilate  the  area. 

10.  Dip  cotton  ball  in  etchant  solution  and  rub  the  PVDF  to  remove  metallic  coating.  You  will 
need  only  a  little  solution.  Do  not  touch  skin  clothes  or  furniture. 

11.  After  the  PVDF  is  etched  and  metal  coating  is  dissolved,  rinse  in  water  to  stop  the  reaction. 

12.  Pat  with  paper  towel  to  dry  and  hang  to  dry  completely. 

13.  Remove  the  permanent  marker  with  the  alcohol  and  allow  to  evaporate. 

14.  Store  loosely  rolled  (do  not  fold)  and  away  from  dust  and  contamination. 
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A. 2  Tensioning  the  Membrane 

1.  Place  the  membrane  over  the  large  ring  in  the  following  order  (bottom  to  top):  bottom  of 
clamp,  12  inch  aluminum  ring,  rubber  o-ring,  PVDF  membrane,  aluminum  top  ring,  top  of 
clamp. 

2.  Space  clamps  at  equal  increments  around  the  circumference  of  the  12  inch  ring.  Do  not  tighten 
yet. 

3.  Gently  pull  all  wrinkles  out  of  the  membrane.  Pull  in  all  directions,  multiple  times.  Be  sure 
not  to  pull  the  PVDF  material  out  from  between  the  o-ring  and  the  aluminum  top  ring  (the 
PVDF  must  rest  on  the  o-ring  at  every  point  around  the  circumference). 

4.  Slowly  tighten  the  clamps  in  opposing  corners  until  the  slack  is  out  of  the  membrane.  Alternate 
the  corners  after  every  eighth  turn  of  the  clamp  handle. 

5.  Tighten  until  the  PVDF  is  stretched  tightly  in  all  directions.  Typically,  until  it  is  hard  to 
tighten  the  clamp  handles  anymore  with  only  two  fingers. 
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A. 3  Epoxy  the  PVDF  membrane  to  the  6  inch  ring 

1.  Begin  with  an  etched  piece  of  PVDF  stretched  over  the  12  inch  ring  and  under  appropriate 
tension. 

2.  Be  sure  that  the  top  of  the  membrane  has  windows  etched  in,  so  that  you  can  easily  center  the 
aluminum  ring. 

3.  Clean  and  prep  the  grooved  area. 

4.  Unscrew  the  wire  support  plates  on  the  top  of  the  6  inch  rings  (2  or  3  plates,  2  hex  screws 
each).  This  will  allow  the  ring  to  sit  more  level. 

5.  Be  sure  to  wear  gloves  and  be  in  a  well  ventilated  area. 

6.  M-bond  610  epoxy  was  used. 

7.  Use  the  plastic  dropper  to  carefully  fill  the  inner  groove  on  the  6  inch  ring. 

8.  Fill  until  the  shallow  groove  at  he  center  of  the  other  two  grooves  is  completely  covered  with 
epoxy. 

9.  Any  excess  epoxy  will  run  around  the  outer  groove  and  drain  out  the  side,  so  have  a  towel 
there  to  catch  the  runoff  (keep  the  ring  level  though). 

10.  Allow  the  epoxy  to  set  up  for  about  half  an  hour,  so  that  the  ring  can  be  turned  up  side  down 
without  the  epoxy  running  all  over  everything. 

11.  Use  the  epoxy  brush  to  coat  the  PVDF  with  epoxy  in  the  area  that  will  contact  the  ring.  Do 
not  coat  any  PVDF  that  will  be  on  the  interior  of  the  ring  (the  membrane  dynamics  could  be 
affected). 

12.  Flip  the  aluminum  ring  over  and  gently  place  it  on  the  stretched  PVDF.  Line  the  sides  of  the 
ring  up  with  the  pattern  shown  in  the  PVDF  window. 

13.  Allow  the  aluminum  ring  to  cure  to  the  membrane  for  two  days. 

14.  Gently  loosen  the  clamps  and  lower  the  membrane  onto  the  table. 

15.  Remove  all  tensioning  hardware  and  examine  the  epoxy  bond  of  the  membrane 

16.  Carefully  cut  the  excess  PVDF  material  away  from  the  edge  of  the  aluminum  ring.  Be  careful 
to  leave  all  of  the  electrical  leads  in  tact  and  to  leave  a  space  for  the  ground  wire  on  the  top 
of  the  PVDF  sheet.  Sharp  scissors  work  best. 

17.  Place  the  membrane  upside  down  in  the  oven  for  24  hours  at  150°  F  to  cure  the  epoxy. 

18.  Adhere  the  copper  lead  wires  to  the  PVDF  material.  Cover  with  electrical  tape  for  safety  and 
don’t  forget  the  ground  lead. 
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19.  Membrane  is  ready  for  coating. 
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A. 4  Coating  the  Membrane  with  Silicon  Rubber  Resin 

1.  Thoroughly  clean  the  surface  with  an  air  jet  cluster.  Silicon  rubber  is  very  susceptible  to 
cont  animation . 

2.  Suspend  the  membrane  on  supports  around  the  circumference  of  the  6  inch  ring.  This  method 
should  create  a  flat  surface  that  compensates  for  the  sag  in  the  membrane  due  to  gravity.  The 
other  method  is  to  set  the  membrane  on  a  very  smooth  glass  surface,  which  will  create  a  more 
evenly  distributed  coating,  but  will  deform  once  lifted  for  testing. 

3.  Prime  the  surface  of  the  PVDF  with  GE  Silicones  (SS4120  01P).  Simply  paint  this  substance 
over  the  entire  surface  of  the  membrane.  This  will  prevent  the  silicon  from  separating  from 
the  surface  due  to  air  bubble  or  impurities  on  the  surface. 

4.  Mix  two  test  tubes  of  the  GE  Silicon  (  Ten  parts  RTV615A  resin  (by  weight)  and  one  part 
RTV615B  curing  agent).  Stir  for  ten  minutes,  transfer  into  another  mixing  dish  and  stir  for 
another  ten  minutes.  Pour  into  the  test  tubes  and  place  in  a  centrifuge  for  ten  minutes  to 
remove  all  bubbles. 

5.  Carefully  pour  the  mixture  onto  the  clean  membrane  surface.  Pour  in  the  center  of  the  mem¬ 
brane,  with  the  lip  of  the  test  tube  as  close  as  possible  to  the  membrane  surface.  Any  distur¬ 
bance  will  cause  bubble  that  may  be  hard  to  remove  later.  Pour  only  the  minimum  amount  of 
resin  needed  to  coat  the  surface. 

6.  Tilt  the  membrane  smoothly  to  completely  coat  the  membrane  surface. 

7.  Pour  any  excess  silicon  off  of  the  membrane  surface.  This  will  ensure  the  thinnest  coating  that 
is  possible. 

8.  Carefully  wipe  the  excess  silicon  rubber  off  of  the  walls  of  the  6  inch  ring. 

9.  Optional:  Use  a  vacuum  chamber  to  suck  out  the  majority  of  the  bubbles.  This  technique 
was  attempted  for  the  first  7-patch  membrane,  however  the  longer  the  silicon  is  under  suction, 
the  more  its  composition  seemed  to  change.  This  seems  to  be  a  very  fragile  step,  to  much 
suction  seemed  to  create  wrinkles  in  the  surface,  but  no  suction  may  leave  entrained  bubbles 
in  the  resin.  This  step  is  assumed  responsible  for  the  wrinkled  membrane  surface  and  should 
be  exercised  with  caution. 

10.  Place  the  membrane  back  on  the  support  post. 

11.  Cover  the  silicon  with  a  glass  shield.  This  will  keep  the  dust  and  other  contaminants  away 
from  the  mirror  when  it  is  curing. 

12.  Allow  to  cure  for  3  to  seven  days  at  room  temperature. 
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13.  Coat  with  gold  in  the  vacuum  coating  chamber.  This  will  put  a  shiny  reflective  surface  on  the 
membrane. 

Research  shows  that  a  250/im  surface  depression  can  be  reduced  to  approximately  1  /jin  with 
the  first  coat  [14].  Visible  depressions  on  the  surface  have  a  magnitude  of  tens  of  microns  at  least 
(minimum  that  the  human  eye  can  discern),  so  this  is  motivation  for  coating  the  surface  with  at  least 
one  coating.  A  reflective  layer  placed  directly  on  the  PVDF  would  have  to  have  surface  imperfections 
of  this  order  of  magnitude. 
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A.  5  Turning  on  the  Laser 

1.  Turn  on  the  pressurized  air  table.  Maintain  25  psi. 

2.  Turn  on  the  WaveScope®  computer. 

3.  Turn  on  the  WaveScope®  electronics  box. 

4.  Start  WaveScope®  software  (link  on  desktop). 

5.  Remove  all  reflective  objects  (jewelry,  watches  and  hard  rank)  from  the  room. 

6.  Verify  that  the  door  is  locked  and  closed. 

7.  Put  on  protective  eye  wear  (correct  for  the  wavelength  of  the  laser). 

8.  Turn  off  the  lights. 

9.  Turn  on  the  laser. 
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A. 6  ‘ Real  Time'  Data  Collection  with  the  MATLAB  Computer 

1.  Turn  on  WaveScope®  computer. 

2.  Turn  on  WaveScope®  electronics  box. 

3.  Turn  on  Matlab®  computer. 

4.  Turn  on  all  clSpace®  hardware. 

5.  Make  sure  that  the  Nport®  ethernet  to  serial  adapter  is  plugged  in  and  properly  configured.1 

6.  Start  the  WaveScope®  software  on  the  WaveScope®  computer.  You  should  hear  the  sensing 
unit  cycle  through  different  lenses  at  startup. 

7.  Start  Matlab®  on  the  clSpace®  computer.2 

8.  Setup  the  test  subject,  turn  on  the  laser  and  calibrate  the  WaveScope®  as  normal. 

9.  Open  test  with  at  least  a  Zernike  display. 

10.  Open  the  TCL  command  prompt,  from  the  Misc  menu  on  the  WaveScope®  toolbar. 

11.  Start  a  live  display.3  Data  should  now  be  streaming  to  the  dSpace®  computer.  This  is  evident 
by  the  rapidly  flashing  lights  on  the  Nport®  ethernet  to  serial  connector. 

12.  On  the  Matlab®  computer,  initialize  the  serial  port  in  Matlab®.  The  Seriallnit.m  (Sec¬ 
tion  B.ll)  script  can  be  used  for  this. 

13.  Grab  and  process  the  streaming  data  with  the  Matlab®  computer.  The  SerialGrab.m  (Sec¬ 
tion  B.10)  and  SerialPlot.m  (Section  B.12)  scripts  can  be  used.  The  WavefrontZernikes.m 
(Section  B.18)  script  is  needed  to  make  the  Zernike’s  for  multiplication  with  the  coefficients 
grabbed  from  the  WaveScope®  computer. 

1Proper  configuration  for  the  Nport®  has  all  ‘dip’  switches  off,  the  device  is  in  TCP  Client  mode  and  the  data 
flow  control  is  turned  off. 

2 dSpace®  was  never  utilized  to  close  a  control  loop,  however  it  is  definitely  a  tool  that  warrants  further  study. 

3 The  data  collection  variables,  such  as  frame  rate  and  number  of  Zernike’s  calculated,  can  still  be  configured  from 
within  the  WaveScope®  software. 
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Appendix  B.  MATLAB  Code 


B.l  AvgRead 

This  code  is  used  to  process  the  data  from  the  laser  vibrometer.  It  reads  in  data  text  files  and 
was  used  to  create  Figures  22  to  24. 


I’ll’llVlTl 

\ l  Name:  AvgRead. m 

in 

1  This  program  opens  up  the  Average  FFT  and  Coherence  data 
1  files  collecetd  from  the  scanning  laser  vibrometer  and 
1  creates  a  plot.  Data  must  be  exported  using  ascii  format. 

Ill 

1  Author:  Eric  M.  Trad 

111 


clc ; 

clear  all; 
close  all; 
format  short  g; 

%Get  file  to  read  the  FFT  Files  for  Laser  Vibrometer 

[File, Path]  =  uigetfile(’*.txt ’, ’Choose  File  for  Average  Spectrum’); 

[Freq, Mag]  =textread (File ,  ’‘/0f  4/0f  ’ ,  ’headerlines’  ,5)  ; 

[CFile ,CPath]  =  uigetfile(’*.txt’ , ’Choose  File  for  Coherence’); 
[CFreq,CMag]=textread(CFile, ’%f  %f ’ , ’headerlines’ ,5) ; 

figure 

subplot (2, 1,1) 
plot(Freq,  Mag,  ’k’); 
xlabel( ’Frequency  (Hz)’) 
ylabel( ’Magnitude  (dB)’) 
subplot (2, 1,2) 
plot(CFreq,  CMag,  ’k’); 
xlabel( ’Frequency  (Hz)’) 
ylabel ( ’ Coherence ’ ) 

%Save  to  the  clipboard 
print  -dbitmap 
close  all 
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B.2  BesZeros 


This  code  is  used  to  find  the  zeros  of  the  Bessel  function. 


It  returns  the  values  in  Table  6. 


function  BesselTable=BesZeros (Mmax , Nmax) 


111111111 

1  BesZeros. m 

111 

1  This  function  will  return  a  matrix  of  values 
1  that  zero  the  bessel  function  of  the  first  kind 
1  of  size (BesO)= [Nmax, Mmax] . 

Ill 

1  VARIABLE  DEFINITIONS 
1  BesO 

1  Rows  are  n  =  1,2,3... 

1  Coulmns  are  m=0 ,1,2. . . 

Ill 

1  Author:  Eric  M.  Trad 

HI 


BesselTable=zeros (Nmax, Mmax) ; 

m=0 ; 
nu=0; 

while  m  <=  Mmax 
Starter=0; 
n=l ; 

%Fill  the  column  of  n  values  to  Nmax 
while  n  <=  Nmax 

TestString=linspace (Starter , Starter+30 , 50000) ; 
BesO=besselj (nu, TestString) ; 

BesselTable (n,m+l)=NextZero(TestString,BesO) ; 
Starter=BesselTable(n,m+l) ; 
n=n+l ; 

end 
m=m+l ; 
nu=nu+l ; 

end 
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B.3  DynamicDataFit 

This  code  is  used  to  process  the  dynamic  data  in  Section  5.6.  It  was  used  to  create  Figures  64 

to  66. 


ii’anni. 

0/oDynamicDataFit  .m 

'/  This  program  plots  a  bar  chart  of  a  Zernike  in  time.  This 
'/  provides  a  comparison  for  the  read  rate  of  the  system  and  other 
'/  dynamic  qualities. 

'/'/'/ 

l  Author:  Eric  M.  Trad 

Til 


clc 

clear  all 
close  all 
format  short  g 

File= ’ 7PatchB_Patch5_HalfHz_600V * ; 
load(File) 

Freq=.5; 

'/Create  the  time  hack 

Hack=Timer(l  :42:4200)/42;  °/0Pick  out  one  time  hack  per  Zcoef  Set 
Hack=Hack-min(Hack)  ;  '/Shift  to  the  origin 

Zc=Zcoef  (4,  : )  ;  '/Look  at  the  first  Zcoef  value  as  it  changes  in  time 
[ZcMax , Zclndex] =max (Zc) ; 

Zc=Zc/ZcMax;  '/Normalize  between  -1<Zc<1 

'/Find  a  Best  Fit  to  Data 
A=1 ; 

w=2*pi*Freq; 
t=linspace (0,25, 1000) ; 

Nlast=10; 

for  phi=0: .0001 :2*pi 

f=abs(A*cos(w*Hack+phi)-Zc) ; 

N=norm(f ) ; 
if  N<=Nlast 
Nlast=N ; 

FinalPhi=phi ; 

end 

end 

f=A*cos(w*t+FinalPhi) ; 


figure (1) 

plot ( t , f , ’b- ’ , Hack , Zc , ’ or : ’ ) 
xlabelC’Time  (s)’) 

ylabelC  ’4~{th}  Zernike  Value,  Normalized  by  L_\infty  Value’) 
legend( ’ Input  Function’ , ’Experimental  Values’, 4) 
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°/0Look  at  the  first  set  of  Zcoef 
figure (2) ; 
for  iii=l:16 

subplot (4,4, iii) 
bar (Zcoef (1 : 10, iii) ) 
axis( [.5  10.5  -1  1]) 
end 

text (-35, 9. 5, ’Bar  Charts  for  the  l“{st}  10  Zernike  Coefficients,  Frames  1:16’) 
figure (1) 
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B.J^  DynamicVisual 

This  code  is  used  to  process  the  dynamic  data  in  Section  5.6.  It  can  be  used  to  create  animations 
of  the  surface  deforming.  This  code  was  used  for  presentations,  but  not  for  any  figure  in  this 
document. 


‘/Dynamic  Visual  .m 

III 

l  This  program  plots  each  frame  of  the  surface  deflection. 

I  It  can  be  used  to  create  animations.  Requires  ’plotman’ 
l  function  to  plot  the  surface.  Use  movie2avi()  to 
l  create  the  avi’s. 

Ill 

l  Author:  Eric  M.  Trad 

HI 


niuiiii 

l  Initialize 

HI 

clc 

clear  all 
close  all 
format  short  g 

global  moviefile  b 


lllllllll 

°/0Load  File 

HI 

l  [File, Dir] ^igetfileC’*. mat’ , ’Choose  the  MAT  file  you  want  to  load’); 
l  load(strcat (Dir , File)) ; 
l  disp(File); 

load( ’ 7PatchB_Patch5_2Hz_600V . mat ’ ) 

[tmpl,  tmp2,  tmp3] =size(Zsurf ) ; 

Zs=Zsurf ( : , : , 2 : tmp3-l) ; 


lllllllll 

°/0Expand  the  data  with  zeros  to  fill  the  membrane 

HI 

R=28-16; 

P=48; 

radius  =  linspace(l ,2.5/1 .5,R+1) ; 
theta  =  [-pi: (2*pi)/P:pi] ; 
for  i  =  1:R+1; 
for  j  =  1:P+1; 

r  =  radius (i); 
t  =  theta(j )  ; 
x(17+i, j)  =  r*cos(t) ; 
y(17+i , j)  =  r*sin(t) ; 

Zs (17+i , j , :)  =0; 
end 
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end 


mxiiixi 

“/.Multiply  by  in2m()  and  1.5  to  fit  to  actual  membrane  size 

111 

x=in2m(x*l .5)*le3;  °/0Want  the  x  axis  in  mm 
y=in2m(y*l . 5)*le3;  “/.Want  the  y  axis  in  mm 


%%  T/:/:/:/:a 

“/.Shift  the  data  one  or  more  points,  rotate  to  align  with  control  outline 

“/.“/.“/. 

S=0;  °/0Shifter,  clockwise,  0<S 
Zs_shift=zeros(size(Zs)) ; 
for  kkk=S+l :49 

Zs_shif t ( : , kkk-S , : ) =Zs ( : , kkk , : ) ; 

end 

for  kkkk=l:S 

Zs.shif t ( : , 49- (S-kkkk) , : ) =Zs ( : , kkkk+1 , : ) ; 

end 

Zs=Zs_shift ; 

'/.PLOT  THE  ABSOLUTE  MAX  DEFORMATION 

for  b=l : 98 

plotman (x , y , Zs , b)  ; 

movief il(b)=getframe(l)  ; 

end 
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B.5  MagPhaseRead, 

This  code  is  used  to  process  the  data  from  the  laser  vibrometer.  It  reads  in  data  text  files  and 
was  used  to  create  Figures  28  to  30. 


i  Name:  MagPhaseRead.m 

Hi 

i  This  program  opens  up  the  magnitude  and  phase  surface  plot 
i  data  files  collected  from  the  scanning  laser  vibrometer 
i  and  creates  a  plot .  Data  must  be  exported  using  ascii  format . 

ill 

i  Author:  Eric  M.  Trad 

Hi 


clc ; 

clear  all; 
close  all; 
format  short  g; 

'/.Get  file  to  read  the  FFT  Files  for  Laser  Vibrometer 

[File, Path]  =  uigetfile(’*.txt’ , ’Choose  File  for  Magnitude  and  Phase’); 
[u,  v,  Mag,  Phase] =textread (File, ’%f  °/.f  if  if  ’,  ’headerlines  ’,  11)  ; 

'/.Shift  the  center  of  the  circle  to  x=0,  y=0 
Ushif  t=min  (u)  +  (max  (u)  -min  (u)  )  /2 ; 

Vshif  t=min  (v)  +  (max  (v)  -min  (v)  )  /2 ; 
u=u-Ushift;  '/.Center  the  data 
v=v-Vshift;  '/.Center  the  data 

'/Scale  the  axes  between  -1  and  1 
u=u. /max(u) ; 
v=v./max(v) ; 

tmp  =  linspace(min(min(u,v) ) ,max(max(u, v) ) ,30) ; 

[U,V]  =  meshgrid(tmp, tmp) ; 

figure (1) 
subplot (2, 1,1) 

M  =  griddata(u, v,Mag,U, V) ; 

surf (U,V,M) 

xlabel( ’x-dir ’ ) 

ylabel ( ’ y-dir ’ ) 

title ( ’Magnitude  (dB)  Plot’) 

colormap  bone 

colorbar 

shading  interp 

axis([-l  1  -1  1]) 

subplot (2, 1,2) 

P  =  griddata(u, v, Phase ,U,V) ; 
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surf (U,V,P) 

xlabel( ’x-dir ’ ) 

ylabel ( ’ y-dir ’ ) 

title (’Phase  (deg)  Plot’) 

colormap  bone 

colorbar 

shading  interp 

axis([-l  1  -1  1]) 

'/Save  to  the  clipboard 
print  -dbitmap 

'/Close  the  plot 
close  1 
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B.  6  MembraneSurface 

This  code  is  used  to  calculate  the  surface  deflection.  It  is  based  on  Equations  22  to  24. 


function  [x , y , W] =MembraneSurf ace (m , n , Radius , rho , Tension , Frequency) 


111111111 

1  MembraneSurface  .m 

z:/:i 

1  This  function  will  return  the  surface  deflection  of  a  circular 
1  membrane  made  of  PVDF,  given  a  tension  and  a  radius. 

Ill 

1  VARIABLE  DEFINITIONS  AND  UNITS 
1  x=x-coordinate  (m) 

1  y=y-coordinate  (m) 

1  W=vertical  surface  deflection  (m) 

1  Frequency  is  in  radian 
1  tta=theta  (rad) 

111 

1  Author:  Eric  M.  Trad 

Til 


rn  =  30; 
pn  =  48; 

R  =  [0:Radius/rn: Radius] ; 

P  =  [0:2*pi/pn:2*pi] ; 

for  iii=l : length (R)  lx  loop 

for  j jj=l : length (P)  %y  loop 
r=R(iii) ; 
tta=P(jjj) ; 

W(iii , j j j )=inv(sqrt (pi*rho) *Radius* . . . 

besselj (m+1, Frequency/sqrt (Tension/ rho) *Radius))  . . . 

*bessel j (m , Frequency/ sqrt (Tension/ rho) *r  ) * (cos (m*tta) +sin (m*tt a) ) ; 
if  abs(r)>Radius 
W(iii , j j j)=0; 

end 

x(iii,jjj)  =  r*cos(tta); 
y(iii,jjj)  =  r*sin(tta) ; 
end 

end 
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B.7  ModeShape 


This  code  is  used  to  predict  mode  shapes  of  fundemental  natural  frequencies. 


I'a'aTITl 

1  This  program  first  calculates  a  tension  from  an 
1  experimentally  observed  mode  shape.  Next  it  predicts 
’/  the  shape  of  any  other  mode  of  vibration  given 
’/  an  m  and  n  value . 

IXl 

*1  Author:  Eric  M.  Trad 

111 

1  Required  Functions  Include: 

’/  BesZeros.m 

’/  FreqCalc.m 

’/  MembraneSurf ace .m 

’/  NextZero.m 

’/  Tens i onCal c . m 

111 


llllllTll 

1  Initialize 

HI 

clc 

clear  all 
close  all 
format  short  g 
tic 


1:1:1:1:1:1X11 

’/Define  Parameters 

111 

Radius=in2m(6) ;  %(m) 

rho=l . 78e3*52e-6;  %(kg/m~3)*(m)=(kg/nr2)  multiply  by  thickness 


1:1:1:1:11111 

’/Normalize  by  the  observed  mode  shape 

111 

iii=l 

FFF=  [218,434, 547 , 656 , 875 ,979,11 14] ; 

MMM= [0 , 1 , 2 , 0 , 4 , 2 , 0] ; 

NNN= [1 , 1 , 1 , 2 , 1 , 2 , 3] ; 

BBB= [2.405,3. 8314 , 5 . 1365 , 5 . 5206 , 7 . 5878 , 8 . 4158 , 8 . 6529] ; 
Frequency_obs_Hz=FFF(iii) ; 

Frequency_obs=BBB(iii) ;  ’/Frequency  observed  in  the  experiment (rad/s) 
m_obs=MMM(iii)  ;  ’/m  value  for  the  observed  mode  shape 
n_obs=NNN(iii)  ;  ‘Zn  value  for  the  observed  mode  shape 
Tension=TensionCalc (Frequency_obs ,m_obs ,n_obs , Radius ,rho) ; 


1:11:1:1:1111 

’/Find  the  max  of  the  first  mode  to  scale  by 

111 

m=0 ; 
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n=l ; 

BesselTable=BesZeros(m,n) ; 

Omega_mn=FreqCalc (BesselTable , Frequency_obs_Hz , Frequency_obs) ;  / (rad/s) 
[x,y,W_on] =MembraneSurf ace (m,n, Radius, rho, Tension,  . . . 

BesselTable (n,m+l)  )  ;  °/0(m) 

Scaler_x=max (max (abs (x) ) ) ; 

Scaler_y=max(max(abs(y))) ; 

Scaler_z=max (max (abs (W_on) ) ) ; 


ix/ix/xa 

°/0Ex  trap  date  to  Surface  Prediction 

IXi 

clear  x  y  W_on 
m=0 ; 
n=l ; 

BesselTable=BesZeros(m,n) ; 

Omega_mn=FreqCalc (BesselTable ,Frequency_obs_Hz  . . . 
,Frequency_obs) ;  / (rad/s) 

disp(sprintf ( ’The  predicted  frequency  for  this  mode 
is:  °/00.2f  Hz’  ,Omega_mn(n,m+l) ) ) 

[x,y ,W_on] =MembraneSurf ace (m,n, Radius , rho ,  . . . 

Tension, BesselTable (n,m+l)) ;  %(m) 
x=x/Scaler_x; 
y=y/Scaler_y ; 

W_on=W_on/Scaler_z ; 


ix/x/xm 

°/,  PLOT 

IX/ 

x_in=x; 
y_in=y ; 

W_on_in=W_on; 

figure 

subplot (2, 1,1) 

contour (x_in,y_in,W_on_in, 10) 
shading  interp 
colorbar 
grid  on 

xlabel( ’X-dir ’ ) 
ylabel ( ’ Y-dir ’ ) 

subplot (2, 1,2) 

surf (x_in , y_in , W_on_in) 

shading  interp 

colorbar 

xlabel( ’X-dir’ ) 

ylabel ( ’ Y-dir ’ ) 

zlabel( ’Surface  Deflection’) 

axis( [min(min(x_in) ) ,max(max(x_in) ) ,min(min(y_in) ) ,  . . . 
max (max (y_ in) ) ,  min(min(W_on_in) ) ,max(max(W_on_in) )] ) 
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%  Stop  the  clock 

disp(sprintf ( ’ \n  It  took  °/0f  seconds  to  run’ 
’  this  code... end  transmission. \n’ , toe) ) 


B.8  NextZero 


This  code  is  used  to  find  the  next  zero  of  a  function  and  is  called  by  the  BesZero.m  code. 


function  min=NextZero(t ,x) 


111111111 

\ l  NextZero.m 

in 

1  This  function  will  return  the  next  [MIN]  zero  of  a  function. 
1  Note:  Developed  for  use  with  Bessel  Fn  and  not 
1  thouroughly  tested  for  other  functions 
111 

1  Author:  Eric  M.  Trad 

HI 


go=l; 
iii=l ; 
min=0 ; 
x=abs(x) ; 
while  go==l 

iii=iii+l ; 
if  x(iii)<=x(iii-l) 
min=t (iii) ; 
elseif  min==0 
x=x*-l ; 
go=l; 

elseif  abs(x(iii))>=. 1 
x=x*-l ; 
go=l; 

else 

go=0; 

end 

end 


too 


B.9  plotman 

This  function  is  used  to  plot  the  deformation  of  the  membrane  surface  and  is  called  by  Dy- 
n  amic  Visual,  m. 


llllilil/ 

‘/plotman.m 

1  This  function  plots  the  deformation  of  the  membrane  surface 
'/  and  is  called  by  Dynamic  Visual  .m  to  create  movies. 

in 

1  Author:  Eric  M.  Trad 

111 

function  plotman (x,y,Zs, frame) 


lllll'll'll 

'/PLOT  the  Surface  Shape 

111 

hh=f igure(l) ; 

1  surf  (x,y ,Zs(  :  , :  , frame) )  ;  '/Plot  whole  surface 

surf  (x(l :  17,  :  )  ,y (1 : 17,  :  )  ,Zs  (1 : 17,  :  ,  frame)  )  ‘/Plot  only  interior  region 
xlabel( ’X-dir  (mm)’); 
ylabel( ’Y-dir  (mm)’); 

'/  axis (  [-in2m(2 . 6) *le3  in2m(2.6)*le3  -in2m(2 . 6) *le3  in2m(2 . 6) *le3  -2.7  2.7]); 
axis ([-30  30  -30  30  -2.7  2.7]); 
caxis([-2.7  2.7]); 
colorbar ; 

'/  view(0,90)  ; 
colormap(getmap) ; 
shading  interp; 
grid  off; 


llllllTll 

1  Draw  the  control  pattern 

111 

1  hold  on; 

'/  plot3(x(7,  :  )  ,y(7,  : )  ,Zs(7,  :  ,  frame)  ,  ’k.’)  ‘/Plot  the  center  circle 

'/  plot3(x(:  ,37)  ,y(:  ,37)  ,Zs(:  ,37, frame)  ,  ’k.  ’)  ‘/Plot  the  center  patch  control  lead 

1  rad= [1:4,6:12,14:20,22:28,30:36,38:44,46:48]; 

'/.  rad2=[4, 6, 12, 14,20,22,28,30,36,38,44,46]  ; 

1  plot3(x(10,rad) ,y(10,rad) ,Zs( 10, rad, frame) , ’k. ’)  ‘/Draw  r=0.75  inch  lines 
1  plot3(x(20,rad) ,y(20,rad) ,Zs (20, rad, frame) , ’k. ’)  ‘/Draw  r=1.75  inch  lines 
1  plot3(x(10:20,rad2) ,y(10:20,rad2) ,Zs(10:20,rad2, frame) ,  ’k.’)  /Draw  lower  radials 
1  plot3(x(10:20,rad2) ,y(10:20,rad2) ,Zs(10:20,rad2, frame) ,  ’k.’)  ‘/Draw  upper  radials 
1  hold  off; 
drawnow 
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B.10  Serial  Grab 


This  code  is  used  to  grab  the  wavefront  data  from  the  serial  port  on  the  Matlab®  computer. 

function  A=SerialGrab2() 


II 

1  Eric  M.  Trad 
•/•/•/ 

1  Date :  14  Jan  05 

III 

1  This  function  grabs  Zernike  Coefficients  that  are  streaming 
1  through  a  serial  connection.  The  serial  connection  must  be 
1  defined  as  type  SERIAL  in  the  variable  ’port’  (use  the 
1  Seriallnit.m  file  to  accomplish  this).  Function  returns  a 

1  vector  of  data  equal  in  length  to  the  number  set  in  ’NumZern’ 

1  variable  and  should  match  the  number  of  coefficients  sent  from 
1  the  Wavescope  side. 

HI 


clear  A  data 


global  NumZern  port  Frames  Timer 


111111111 

1  GRAB  THE  DATA  FROM  THE  SERIAL  PORT 

111 

count=l ; 
iii=0; 
corrupt=0 ; 

WholeSet=0; 
c=0 ; 

A=zeros (NumZern, 2, Frames) ; 

cnt=l ; 

tic 

while  WholeSet  ==  0  ‘/.Check  if  the  SAMPLE  set  is  complete 
for  jjj=l: Frames  ‘/Takes  a  whole  SAMPLE  set 

while  corrupt  ==  0  ‘/Check  if  the  ZERN  set  is  complete 
while  count <=NumZern  ‘/Grab  all  the  ZERN  coef’s 
data=f gets (port) ; 

T imer ( cnt ) =t oc ; 
cnt=cnt+l ; 

if  length (data) <=28  |  length (data) >=35 
iii=iii+l ; 
if  iii>2 

1  disp(’Frame  Grab  Disrupted...’) 

end 

else 

A (count , 1 , j j j )=str2num(data(7 : 10) ) ; 
if  A(1 , 1 , j j j)==l 

tmp=length (data) -2 ; 

A (count ,2, j j j)=str2num(data(16 : tmp) ) ; 


102 


count=count+l ; 


end 

end  ‘/.Check  the  ZERN  set 
end  ‘/.End  Grab 

‘/.Check  for  a  corrupted  set 
for  tester=l :NumZern 

if  A(tester,l, jjj)~=tester 
%  disp( ’Corrupt  Set...’) 

c=c+l ; 
corrupt=0; 
count=l ; 
iii=0; 
break 

else 

corrupt=l ; 

end 

end  ‘/Corruption  Test 
end  ‘/.While  Corrupted 
corrupt=0 ; 
if  c>0 

break 

end 

end  ‘/.Take  the  whole  SAMPLE  set 
if  c  >0  ‘/.Check  the  SAMPLE  set 
%  disp(’Not  a  Whole  Set...’) 

WholeSet=0; 

c=0; 

else 

WholeSet=l ; 

end 

end  ‘/.While  SAMPLE  set  is  not  complete 

disp(sprintf  ( ’Grabbed  °/»i  consecutive  sets  of  Zernike  Coeff  icients ’  .Frames) ) 
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B.ll  Seriallnit 

This  code  initializes  communication  with  the  serial  port  on  the  acquisition  computer. 

function  port=Serial!nit () 


n 

1  Eric  M.  Trad 

in 

1  Date :  14  Jan  05 

HI 

1  This  function  initializes  the  SERIAL  type  variable  port  for  use 
1  with  the  SerialPlot  and  SerialGrab  M-files.  Note:  Occasionally 
1  there  are  instances  where  the  COM  PORT  will  not  be  properly 
1  initialized,  usually  running  this  script  again  will  take  care  of 
1  the  problem . 

Ill 

°/0Note:  Be  sure  to  set  port=SerialInit ()  at  the  command  line  when  running, 
1  because  all  other  scripts  assume  this  name  for  the  connection. 

HI 


111111111 
1  INITIALIZE 
HI 

clear  all; 
close  all; 
format  short  g; 
clc ; 


111111111 

1  INITIALIZE  AND  OPEN  COM  PORT 

HI 

1  Sets  the  COM  PORT  to  C0M1 

1  Sets  the  Baud  Rate  to  115200,  which  is  the  max  on  the  Nport 
1  Sets  the  ByteOrder  to  Little  End  In 
1  Sets  the  Timeout  to  10  seconds 
1  Sets  the  Terminator  to  Carriage  return 

1  Pauses  3  seconds  to  allow  time  to  establish  the  connection  before  opening 
1  the  port 

HI 

port  =  serial ( ’C0M1 ’ , ’BaudRate’ ,115200, ’ByteOrder’ , ’LittleEndian’  ... 

, ’TimeOut ’ , 10, ’Terminator’ , ’CR’ ) ; 

pause (3) ; 
f open (port) ; 
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B.12  SerialPlot 


This  code  is  used  to  acquisition  the  wavefront  data  on  the  Matlab®  computer. 


rarara 

%  Eric  M.  Trad 

ZZZ 

%  Date :  14  Jan  05 

%  Update:  15  Feb  05  -  Include  Matlab  v6  export  and 
%  coefficient  grab. 

zzz 

l  This  script  sets  up  environment  varibles  and  runs  a  number 
%  of  functions  to  grab  Zernike  data  from  the  AFIT  Wavescope 

%  box,  process  it  and  display  it  in  ’real  time’. 

ZZZ 

%  NOTE:  This  script  uses  the  make_zernikes .m  function 
%  to  make  the  Zernike  definitions. 

zzz 


%clear  all; 
close  all; 
format  short  g; 

clear  A  Movie  Zcoef  Zpoly  Zsurf  filename  h  hh  i  j  k  x  y  Timer; 
clc ; 
tic ; 

global  NumZern  port  Frames  Timer 


ZZZZZZZZZ 

%  Initialize 

ZZZ 

NumZern=42 ; 

Framelnterval= . 1 ;  ^Seconds  from  Wavescope 

%  ShowZern=questdlg( ’Would  you  like  to  show  the  Zernike  Component  Plots?’,  ... 
I  ’Question’ , ’No’) ; 

ShowZern=’No’ ; 

%  NewRef=questdlg( ’Would  you  like  to  take  a  static  reference?’,  ... 

%  ’Question’ , ’No’) ; 

NewRef=’No’ ; 

%  NewRef=’ Yes ’ ; 

%  MakeMovie=questdlg( ’Would  you  like  to  make  a  movie  and  save  the  run?’,  ... 

%  ’ Question’ , ’Yes ’) ; 

MakeMovie=’ Yes ’ ; 

Frames=100; 

Pnum=3 ; 


ZZZZZZZZZ 

Z  Call  the  fn  to  put  the  polynomials  together 

ZZZ 

R  =  16; 

P  =  48; 

radius_area  =  1.0; 

[x,y, Zpoly]  =  Wavefront_Zernikes(R,P,radius_area) ; 
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if  strcmp(ShowZern, ’Yes ’ )==1 ; 
figure (1) 
for  i=l:NumZern; 

subplot (5 , ceil (NumZern/5) , i) ; 
mesh(x,y,Zpoly( : , : , i) ) ; 

end 

end 


\ l  Grab  a  reference  image 

m 

if  strcmp(NewRef , ’ Yes ’ )==1 

h=warndlg( ’Prepare  to  grab  static  reference.  Press  OK  to  continue. 
waitfor(h) ; 

Reasonable= ’  No  ’ ; 

while  str cmp (Reasonable, ’No’) ==1 
clear  Atmp  Aflat  Zsurf 
Atmp=SerialGrab2 ; 

Af lat=mean(Atmp( : ,2, : ) ,3) ; 
figure (3) 

Zsurf =Zpoly( : , : , 1)*0; 
for  i=l:NumZern; 

Zsurf =Zsurf+Aflat (i)*Zpoly( : , : ,i) ; 

end 

surf (x , y , Zsurf ) 
shading  flat 
myclmp=getmap ; 
colormap(myclmp) ; 
colorbar  southoutside 
°/,axis([-l  1-11-55]) 
drawnow 

Reasonable=questdlg( ’Does  the  reference  frame  look  reasonable?’ 
’Question’ , ’No’) ; 

end 

close  3 

end 

if  strcmp(NewRef , ’ Yes ’ )==1 

h=warndlg( ’Prepare  to  grab  dynamic  frames.  Press  OK  to  continue.’) 
waitfor(h) ; 

end 


mmm 

%  Grab  moving  frames 

m 

hh=f igure(2) ; 

A=SerialGrab2 ; 
for  j=l:Frames 

Zcoef ( : , j)=A( : ,2, j)-Af lat ;  %  Zernike’s  from  Wavescope 
Zsurf ( : , : , j)=Zpoly(: , : ,1)*0; 
for  i=l :NumZern; 

Zsurf ( : , : , j )=Zsurf ( : , : , j)+Zcoef (i , j )*Zpoly ( : , : , i) ; 
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end 


end 


iiT/:/:/:/:a 

%  Display  the  moving  frames  and  make  the  movie 

m 

for  k=l: Frames 

surf (x , y , Zsurf ( : , :  ,  k)  ) 
shading  interp 
colormap  bone 
colorbar  southoutside 
axis( [-1  1-11-22]) 

%  view (0,90); 

drawnow 

%  Movie (k) =getframe (hh) ; 

end 

close 

Timer =42*Timer (length(Timer)-Frames*NumZern: length(Timer) ) ; 
FRate=mean(Timer(2:Frames)-Timer(l :Frames-l)) ; 

FRateStd=std (Timer  (2 :  Frames) -Timer  ( 1 :  Frames-1)  )  ; 
disp(sprintf ( ’ \nMean  Time  Between  Samples:  %0.2g  s’,FRate)) 
disp(sprintf (’Std  Dev  of  Time  Between  Samples:  °/00.2g  s ’  ,FRateStd) ) 
disp(sprintf ( ’Frame  Rate:  %0.2g  Hz  \n’ , 1/FRate) ) 

if  str cmp (MakeMovie , ’ Yes ’ ) ==1 

filename=inputdlg( ’Please  enter  a  file  name: ’, ’Question’) ; 

%  disp( ’Creating  Movie...’); 

°/0  movie2avi  (Movie ,  char  (filename)  ,  ’fps’  ,5)  ; 

disp(’Saving  Workspace  Variables  in  Matlab  v6  Format...’); 
save (char (filename) , ’-v6’ ) ; 

end 

disp(’ Script  is  done. . .END  TRANSMISSION’) 
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B.13  StaticBarGenerate 


This  code  is  used  to  process  the  static  data  in  Section  5.5.4.  It  was  used  to  create  Figure  33. 


i:i:i:aTU:i 

%St at icBar Generate .m 

111 

1  This  program  plots  a  bar  chart  of  the  Zernike’s.  For  comparison 
1  any  number  of  data  sets  can  be  called  up  for  ploting  in  a  group. 
1  Also  plots  a  spiff y  Zernike  component  breakdown. 

HI 

1  Author:  Eric  M.  Trad 

HI 


clc 

clear  all 
close  all 
format  short  g 

repeat=’Yes’ ; 
ii=l ; 

while  strcmp (repeat , ’ Yes ’ ) 

°/0Load  File 

[File ,Dir] =uigetfile( ’* .mat * , ’Choose  the  MAT  file  you  want  to  load’); 
load(strcat (Dir, File)) ; 
disp(File) ; 

Zc( : ,ii)=Zcoef ; 

repeat=questdlg( ’Would  you  like  to  add  another  data  set’); 
ii=ii+l ; 

end 


IHIHIH 

°/0PLOT  the  bar  chart  of  the  Zernike ’s 

111 

figure (1) ; 

h=bar (Zc (1:42,:),’ group ’ ) ; 

colormap  bone 

axis ([0,42, -.75, .75]) 

xlabel( ’Zernike  Coeficient  Number’) 

ylabel( ’Coefficient  Value’) 

print  -dbitmap 


iininn 

1  Plot  the  Zernike  Polynomial  Components 

HI 

figure (2) 

NumZern=25;  °/0Best  if  multiple  of  5 
for  i=l:NumZern; 

subplot (5 , ceil (NumZern/5) , i) ; 
surf (x,y ,Zpoly( : , : ,i)) ; 
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shading  interp 
view(0,90) 
colormap (getmap) 
grid  off 
axis  off 

text (-1  ,-l , sprintf  (’“/oi’  ,i)) 

end 
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B.IJ^  StaticPlotter 


This  code  is  used  to  process  the  static  data  in  Section  5.5.4. 

to  61. 


It  was  used  to  create  Figures  42 


•av/iim. 

ZStaticPlotter .m 

111 

1  This  program  plots  the  surface  deflection,  radial  cut  pattern 
1  and  standard  deviation  plots  for  the  wavescope  data.  Various 
1  max,  min  and  other  info  is  also  calculated. 

HI 

1  Author:  Eric  M.  Trad 

HI 


111111111 

1  Initialize 

HI 

clc 

clear  all 
close  all 
format  short  g 

Flag=0;  °/00=Def lection  Plot,  l=Radial  Cut,  2=Std  Plot 


111111111 

°/0Load  File 

HI 

1  [File, Dir] ^igetfileC’*. mat’ , ’Choose  the  MAT  file  you  want  to  load’); 
1  load(strcat (Dir , File)) ; 

1  disp(File); 

load  7PatchB_Patch7_0Hz_100V_V6 ; 


111111111 

°/0Average  and  STD  the  Surface 

HI 

[tmpl,  tmp2,  tmp3] =size(Zsurf ) ; 
Zavg=mean(Zsurf ( : , : ,2 : tmp3-l) ,3) ; 
Zstd=std (Zsurf ( : , : , 2 : tmp3- 1 ) ,0,3) ; 


111111111 

°/0Expand  the  data  with  zeros  to  fill  the  membrane 

111 

R=28-16; 

P=48; 

radius  =  linspace(l ,2.5/1 .5,R+1) ; 
theta  =  [-pi: (2*pi)/P:pi] ; 
for  i  =  1:R+1; 
for  j  =  1 :P+1 ; 

r  =  radius (i); 
t  =  theta(j ) ; 
x(17+i,j)  =  r*cos(t); 


no 


y(17+i , j)  =  r*sin(t) ; 
Zavg(17+i,j)  =  0; 
Zstd(17+i, j)  =  0; 
end 
end 


111111111 

°/0Multiply  by  in2m()  and  1.5  to  fit  to  actual  membrane  size 

nt 


x=in2m(x*l .5)*le3;  °/0Want  the  x  axis  in  mm 
y=in2m(y*l . 5)*le3;  %Want  the  y  axis  in  mm 


111111111 

°/0Shift  the  data  one  or  more  points,  rotate  to  align  with  control  outline 

111 

S=l;  °/0Shifter,  clockwise,  0<S 
Zavg_shift=zeros(size(Zavg) ) ; 
for  kkk=S+l :49 

Zavg_shift(: ,kkk-S)=Zavg( : ,kkk) ; 

end 

for  kkkk=l:S 

Zavg_shift(: ,49-(S-kkkk) )=Zavg( : ,kkkk+l) ; 

end 

Zavg=Zavg_shif t ; 


111111111 

\ l  Alter  the  filename  for  display  in  title 

i:a 

Fname=regexprep (File , ’ _ ’ , ’  ’); 

Fname=regexprep(Fname , ’ .mat’ , ’ ’ ) ; 

if  Flag==0 

°/0PLOT  the  Surface  Shape 

111 

figure (1) 
surf (x,y ,Zavg) ; 
xlabel( }X-dir  (mm)’); 
ylabel( ’Y-dir  (mm)’); 

axis([-in2m(2.6)*le3  in2m(2 . 6) *le3  -in2m(2 . 6) *le3  in2m(2.6)*le3  -2.7  2.7]); 

caxis([-2.7  2.7]); 

colorbar ; 

view(0,90) ; 

myclmp=getmap ; 

colormap(myclmp) ; 

shading  interp; 

grid  off; 


mmm 

°/0  Draw  the  control  pattern 

17.1 

hold  on; 


in 


plot3(x(7,  :  )  ,y(7,  :  )  ,Zavg(7,  : )  ,  ’k.’)  '/Plot  the  center  circle 

plot3(x( : ,37) ,y( : ,37) ,Zavg( : ,37) , ’k. ’ )  '/Plot  the  center  patch  control  lead 

rad= [1:4, 6: 12, 14: 20, 22: 28, 30: 36, 38:44, 46: 48] ; 

rad2= [4, 6, 12, 14, 20, 22, 28, 30, 36, 38, 44, 46] ; 

plot3(x(10,rad)  ,y(10,rad)  ,Zavg(10,rad)  ,  ’k.  ’ )  '/Draw  r=0.75  inch  lines 
plot3(x(20 ,rad)  ,y(20,rad)  ,Zavg(20,rad)  ,  ’k.  ’)  '/Draw  r=1.75  inch  lines 
plot3(x(10:20,rad2) ,y(10:20,rad2) ,Zavg(10:20,rad2) ,  ’k.’)  “/Draw  lower  radials 
plot3  (x  (10 : 20,  rad2)  ,y  (10:20,  rad2)  ,Zavg(  10:20,  rad2)  ,  ’k.’)  “/Draw  upper  radials 


'/mmn 


'/  Find  the  min  and  max  of  the  surface 
'/'/'/ 


'/Find  the 
[Bminmin, 
[Bminmin, 
‘/Find  the 
[Iminmin, 
[Iminmin, 
‘/Find  the 
[Bmaxmax , 
[Bmaxmax , 
‘/Find  the 
[Imaxmax, 
[Imaxmax, 


absolute  minimum  of  the  surface 

lx]  =min(min(Zavg,  []  ,2))  ; 

ly]  =min(min(Zavg,  []  ,1))  ; 
interior  minimum  of  the  surface 

lx]  =min(min(Zavg(l :  15,  :  )  ,  []  ,2) )  ; 

ly]  =min(min(Zavg(l : 15, :),[], 1) ) ; 
absolute  maximum  of  the  surface 

lx]  =max(max(Zavg,  []  ,2))  ; 

ly]  =max(max(Zavg,  []  ,  1)  )  ; 
interior  maximum  of  the  surface 

lx]  =max(max(Zavg(l :  15,  :  )  ,  []  ,2))  ; 

ly]  =max(max(Zavg(l : 15 ,:),[], 1) ) ; 


disp(sprintf ( ’Min  Disp  =  °/0.5g  \nMax  Disp  =  %0.5g  \n’ , Bminmin, Bmaxmax) ) ; 
disp(sprintf ( ’Min  Interior  Disp  =  %0.5g  \nMax  Interior  Disp  =  °/0. 5g’ , Iminmin, Imaxmax) ) ; 


r/mr/n 

'/  Find  the  Peak-Valley  Displacements 

'/'/.'/ 

PV=abs (Bmaxmax) +abs (Bminmin) ; 

disp(sprintf  (  ’ \nPV  =  °/0.5f  microns  ('/0.5f  wavelengths)  ’  ,PV,PV/ . 633)  )  ; 

PV_ inter ior=abs ( Imaxmax) +abs ( Iminmin) ; 

disp (sprintf ( ’PV_ interior  =  °/0.5f  microns  (%0.5f  wavelengths) \n’ ,PV_interior,PV/ . 633) ) ; 

‘/Copy  to  the  clipboard 
'/  print  -dbitmap 
close  1 

elseif  Flag==l 

'/•/  ix/xm 

'/  PLOT  the  radial  cut  graphs 
'/'/'/ 

figure (2) 
hold  on; 

RRR=linspace (0 , in2m (2 . 5 ) * le3 , length (Zavg ( : , 1 ) ) ) ; 
for  jjj=l:8:48 

errorbar (RRR , Zavg ( : , j  j  j ) , Zstd (:,jjj),’k:’) 

end 

hh=plot (RRR, Zavg ( : ,33) ’ , ’k+’ ,RRR,Zavg( : ,25) ’ , ’ko’ ,RRR,Zavg( : , 17) ’ , ’kx’ , . . . 

RRR , Zavg ( : ,9) ’ , ’ks’ ,RRR,Zavg(: ,1) ’ , ’kd’ ,RRR,Zavg(: ,41) ’ , ’k~’) ; 
legend (hh, ’Patch  2’, ’Patch  3’,  ’Patch  4’,  ’Patch  5’,  ’Patch  6’,  ’Patch  7’,1) 
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set(hh, ’MarkerEdgeColor ’ , ’k’ , ’MarkerFaceColor ’ , ’none’ , ’MarkerSize’ ,6) 

axis([0  in2m(2 . 5) *le3  -3  3]) 

xlabel( ’Radius  (mm)’) 

ylabel( ’Displacement  (\mu  m)’) 

hold  off; 

'/Copy  to  the  clipboard 

print  -dbitmap 

else 

mmm 

°/,PL0T  the  Surface  Standard  Deviation 

7.7.7. 

figure (3) 
surf (x,y,Zstd) ; 
xlabel( ’X-dir  (mm)’); 
ylabelC ’Y-dir  (mm)’); 

axis(  [-in2m(2.6)*le3  in2m(2.6)*le3  -in2m(2 . 6) *le3  in2m(2.6)*le3  0  .5]); 
caxis(  [0  . 5]  )  ; 

'/  title (sprintf  (’“/.s’  ,Fname)) 
colorbar ; 
view(0,90) ; 

'/  colormap  bone; 
myclmp=getmap ; 
colormap (myclmp) ; 

'/  colormap  flag; 
shading  interp; 
grid  off; 


mxiixa 

'/  Draw  the  control  pattern 

7.7.7. 


hold  on; 

plot3(x(7,  :  )  ,y(7,  :  )  ,Zstd(7,  : )  ,  ’k.’)  '/Plot  the  center  circle 
plot3(x(: ,37) ,y(: ,37) ,Zstd(: ,37) , ’k. ’)  '/Plot  the  center  patch  control  lead 
rad= [1:4, 6: 12, 14: 20, 22: 28, 30: 36, 38: 44, 46: 48] ; 
rad2= [4, 6, 12, 14, 20, 22, 28, 30, 36, 38, 44, 46] ; 

plot3(x(10,rad)  ,y(10,rad)  ,Zstd(10,rad)  , ’k. ’)  '/Draw  r=0.75  inch  lines 
plot3(x(20,rad) ,y(20,rad) ,Zstd(20,rad) , ’k. ’ )  '/Draw  r=1.75  inch  lines 
plot3(x(10:20,rad2) ,y(10:20,rad2) ,Zstd(10:20,rad2) ,  ’k.’)  “/Draw  lower  radials 
plot3(x(10:20,rad2) ,y(10:20,rad2) ,Zstd(10:20,rad2) ,  ’k.’)  “/Draw  upper  radials 


nr/mr/ 

'/  Find  the  min  and  max  of  the  surface 

7.7.7. 

'/Find  the  absolute  maximum  of  the  surface 
[Bmaxmax ,  lx]  =max  (max (Zstd ,  []  ,  2)  )  ; 

[Bmaxmax,  Iy]  =max (max (Zstd,  []  ,1))  ; 

'/Find  the  interior  maximum  of  the  surface 
[Imaxmax,  lx]  =max (max (Zstd (1 : 15, :  )  ,  []  ,2) )  ; 
[Imaxmax,  Iy]  =max (max (Zstd (1 : 15, :),[],  1) )  ; 

disp (sprintf  ( ’Max  Disp  =  °/0 . 5g’ , Bmaxmax)  )  ; 
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disp(sprintf ( ’Max  Interior  Disp  =  %0.5g’ ,Imaxmax)) ; 


%Copy  to  the  clipboard 

print  -dbitmap 

end 
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B.15  StaticFramePlot 


This  code  is  used  to  process  the  static  data  in  Section  5.5.4.  It  was  used  to  create  Figure  41. 


i:i:i:aTU:i 

%StaticFramePlot .m 

HI 

1  This  program  plots  the  mean  surface  deflection  or  standard 
1  deviation  of  the  mean  deflection  for  the  static  reference 
1  frame  that  is  selected  in  a  data  file. 

Ill 

1  Author:  Eric  M.  Trad 

HI 


111111111 

1  Initialize 

HI 

clc 

clear  all 
close  all 
format  short  g 


IHlltltl 

°/0Load  File 

111 

[File ,Dir] =uigetfile( ’  * .mat ’ , ’Choose  the  MAT  file  you  want  to  load’); 
load(strcat (Dir, File)) ; 

Flag=0;  1  0=calculate  mean,  l=calculate  std 


111111111 

1  Average  Reference  frames 

111 

if  Flag==0 

Amean=mean(A,3) ; 

else 

Amean=std(A,0,3) ; 

end 

Amean=Amean( : ,2) ; 


111111111 

1  Create  the  Surface  Deformation 
1  from  coefficients  and  polynomials 

HI 

ZZZ=0 ; 

for  i=l:NumZern; 

ZZZ=ZZZ+Amean(i)*Zpoly ( : , : ,i) ; 

end 


111111111 

°/0Expand  the  data  with  zeros  to  fill  the  membrane 

HI 

R=28-16; 
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P=48; 

radius  =  linspace(l ,2.5/1 . 5,R+1) ; 
theta  =  [-pi: (2*pi)/P:pi] ; 
for  i  =  1:R+1; 
for  j  =  1:P+1; 

r  =  radius (i); 
t  =  theta(j ) ; 
x(17+i,j)  =  r*cos(t); 
y(17+i,j)  =  r*sin(t) ; 

ZZZ(17+i , j )  =  0; 

Zstd(17+i , j )  =  0; 

end 

end 

x=in2m(x*l .5)*le3;  %Want  the  x  axis  in  mm 
y=in2m(y*l . 5)*le3;  %Want  the  y  axis  in  mm 


II  mini 

°/0Shift  the  data  one  or  more  points,  rotate  to  align  with  control  outline 

III 

S=l;  “/Shifter,  clockwise,  0<S 
Zavg_shift=zeros(size(ZZZ) ) ; 
for  kkk=S+l :49 

Zavg_shift ( : ,kkk-S)=ZZZ( : ,kkk) ; 

end 

for  kkkk=l:S 

Zavg_shif t ( : , 49- (S-kkkk) ) =ZZZ ( : ,kkkk+l) ; 

end 

ZZZ=Zavg_shif t ; 


111111111 
1  Plot 
HI 

surf (x,y,ZZZ) 
shading  flat 
myclmp=getmap ; 
colormap(myclmp) ; 
shading  interp; 
grid  off; 

xlabel( ’X-dir  (mm)’); 
ylabel( ’Y-dir  (mm)’); 
if  Flag==0 

axis([-in2m(2.6)*le3  in2m(2 . 6) *le3  -in2m(2 . 6) *le3  in2m(2.6)*le3  -2.7  2.7]); 
caxis([-2.7  2.7]); 

else 

axis( [-in2m(2.6)*le3  in2m(2 . 6) *le3  -in2m(2 . 6) *le3  in2m(2.6)*le3  0  0.75]); 
caxis( [0  .75] ) ; 

end 

colorbar 
view(0,90) ; 


111111111 

1  Draw  the  control  pattern 
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Ill 

hold  on; 

plot3(x(7,  :  )  ,y(7,  :  )  ,ZZZ(7,  :  )  ,  ’k.’)  °/0Plot  the  center  circle 
plot3(x(:  ,37)  ,y(:  ,37)  ,ZZZ(:  ,37)  ,  ’k.  ’)  °/0Plot  the  center  patch  control  lead 
rad= [1:4, 6: 12, 14: 20, 22: 28, 30: 36, 38: 44, 46: 48] ; 
rad2= [4, 6, 12, 14, 20, 22, 28, 30, 36, 38, 44, 46] ; 

plot3(x(10,rad)  ,y(10,rad)  ,ZZZ(10,rad)  , ’k. ’)  °/0Draw  r=0.75  inch  lines 
plot3(x(20,rad)  ,y(20,rad)  ,ZZZ(20,rad)  , ’k. ’)  °/0Draw  r=1.75  inch  lines 
plot3(x(10:20,rad2)  ,y(10:20,rad2)  ,ZZZ(10:20,rad2)  ,  ’k.’)  '/.Draw  lower  radials 
plot3(x(10:20,rad2) ,y(10:20,rad2) ,ZZZ(10:20,rad2) ,  ’k.’)  %Draw  upper  radials 
hold  off; 

print  -dbitmap 
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B.16  StaticVoltageDependence 


This  code  is  used  to  process  the  static  data  in  Section  5.5.4. 


It  was  used  to  create  Figure  62. 


i:i:i:i:i:iTi:i 

%StaticVoltageDependence .m 

i:a 

1  This  program  plots  the  change  of  displacement  with  voltage. 
1  The  data  is  gathered  from  other  code  outpus. 

Ill 

1  Author:  Eric  M.  Trad 

111 


clc 

clear  all 
close  all 
format  short  g 

°/0Enter  Data  from  other  calculations  (low  tech) 


MinlOO= [  -0.82823 
-0.3501 
-0.36026 
-0.4821 
-0.55265 
-0.51593 
-0.43385] ; 

Min200= [-1.5915 
-0.28348 
-0.61241 
-0.78587 
-1.1828 
-0.95211 
-0.60629] ; 

Min300= [-1.0929 
-0.49212 
-1.0891 
-1.3505 
-1.3835 
-1.4651 
-0.77172] ; 

Min400= [  -1.0121 
-0.70536 
-1.4521 
-1.6321 
-1.7879 
-2.0969 
-1.1939] ; 

Min500= [-0.99005 
-1.0358 
-1.7944 
-1.8423 
-2.1288 


118 


-2.5875 
-1.4155] ; 

Min600= [-0.94603 
-1.1854 
-2.1269 
-1.7407 
-2.6766 
-2.6243 
-1.5429] ; 

Volt=ones(6,7) ; 

Volt(l, :)=100; 

Volt (2, : )=200; 

Volt (3, : )=300; 

Volt (4, : ) =400 ; 

Volt (5, : )=500; 

Volt (6, : )=600; 

figure 

plot (Volt, [mean(MinlOO) ;mean(Min200) ;mean(Min300) ; . . . 

mean(Min400) ;mean(Min500) ;mean(Min600)] , ’kx-’ ,100:600,0, ’k 

xlabel ( ’ Charge  (volts) ’ ) 

ylabel( ’Average  Surface  Deformation  (\mum)’) 

print  -dbitmap 
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B.17  TensionCalc 


This  code  is  used  to  calculate  the  estimated  tension  in  the  membrane.  It  is  related  to  Equa¬ 
tion  34. 


function  T=TensionCalc (Freq ,m , n , Radius , rho) 


iimiiii 

1  TensionCalc .m 

HI 

1  This  function  will  return  the  tension  as  calculated  by  circular 
"/  membrane  theory.  The  material  is  assumed  to  be  PVDF  and  the 
"/  density  is  set  accordingly 

m 

"/.  VARIABLE  DEFINITIONS  AND  UNITS 

"/  Freq  =  frequency  of  the  tension  calculation  [rad/s] 

"/  m,n  =  describe  the  mode  of  freq 

"/  Radius  =  radius  of  the  circular  membrane  [meters] 

%  T  =  claculated  tension  [newtons] 

111 

1  Author:  Eric  M.  Trad 

111 


111111111 

1  Find  the  Zeros  of  the  Bessel  Fn 

111 

Mmax=m ; 

Nmax=n ; 

BesTable=BesZeros(Mmax,Nmax) ; 
BesVal=BesTable(n,m+l) ; 


111111111 

1  Calculate  the  tension 

HI 

T=rho*(Freq*Radius/BesVal)''2;  "/Tension  [kg/nT2]  [rad~2/s~2]  [m~2]  =  [kg/s~2]  =  [N/m] 
disp(sprintf ( ’ \nThe  tension  is:  %0.5e  N/m  f or (%i,%i)\n’ ,T,m,n) ) 
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B  .1 8  WavefrontZernikes 

This  code  generates  the  Zernike  Polynomials  for  use  with  the  acquisition  code. 


function[x,y,Zpoly]  =  Wavefront _Zernikes(R,P,radius_area) ; 

^/********************************************************************** 

%  PURPOSE:  The  Wavefront_Zernikes  function  generates  a  set  of  basis  elements 
%  for  plotting  the  Zernike  polynomials  from  the  Zernike  Coefficients 
%  generated  by  the  Wavefront  sensor.  The  basis  elements  are  consistent  with 
%  only  the  Wavefront  software,  as  the  ordering  of  the  Zernike  polynomials 
%  by  Wavefront  is  non-standard. 

i 

%  Inputs: 

%  R  :  number  of  radial  cuts 
%  P  :  number  of  angular  cuts 

%  radius_area:  radius  of  cicular  area  (Optional) 

%  The  basis  elements,  Zpoly,  may  be  plotted  by  using  the  surf  command,  i.e. 

°i  surf  (x,y ,Zpoly(  :  ,  :  ,4)  prints  the  fourth  polynomial. 

% 

‘ l  ****  OR  ***** 

%  single  input  R 

%  R  :  [n  x  3]  array  of  coordinates  such  that 
%  R(l,:)  is  the  index 
%  R(2,:)  is  the  radial  coordinate 
%  R(3,:)  is  the  theta  coordinate 

l 

7.  The  basis  elements,  Zpoly,  may  be  plotted  by  using  the  plot3  command, 

%  i.e. 

%  plot3(x,y ,Zpoly( : , : ,4) , ’ . ’ )  prints  the  fourth  polynomial 

l 


7. 


%  WRITTEN 

7. 

°l 

1 

l 

l 


BY :  Maj  Shepherd 

Air  Force  Institute  of  Technology 
Aeronautics  and  Astronautics  Department 
AFIT/ENY 
2950  P  Street 

Wright -Patterson  AFB,  OH  45433-7765 


7. 


I  CREATED:  18  Jan  05 
l  REVISED:  n/a 
l 


if  nargin  >  1 
if  nargin  ==  2 
radius_area  =  1 ; 
end; 

radius  =  [0:radius_area/R:radius_area] ; 
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theta  =  [-pi: (2*pi)/P:pi] ; 
for  i  =  1:R+1; 
for  j  =  1:P+1; 

r  =  radius (i); 
t  =  theta(j) ; 
x(i, j)  =  r*cos(t) ; 
y(i, j)  =  r*sin(t) ; 

Zpoly(i,j,l)  =  r*cos(t);  '/.(X  Tilt) 

Zpoly(i,j,2)  =  r*sin(t) ;  '/»  (Y  Tilt) 

Zpoly(i,j,3)  =  2*r~2-l;  '/.(Focus) 

Zpoly(i,j,4)  =  r~2*cos(2*t) ;  %(0  Astigmatism) 

Zpoly(i,j,5)  =  r~2*sin(2*t) ;  %(45  Astigmatism) 

Zpoly(i,j,6)  =  (3*r~2-2)*r*cos(t) ;  '/.(X  Coma) 

Zpoly(i,j,7)  =  (3*r~2-2) *r*sin(t) ;  °/.(Y  Coma) 

Zpoly(i,j,8)  =  6*r~4-6*r~2+l ;  '/.(Spherical) 

Zpoly(i,j,9)  =  r~3*cos(3*t) ; 

Zpoly(i,  j  ,  10)  =  r~3*sin(3*t) ; 

Zpoly(i, j ,11)  =  (4*r~2-3)*r~2*cos(2*t) ; 

Zpoly(i, j ,12)  =  (4*r~2-3)*r~2*sin(2*t) ; 

Zpoly(i, j ,13)  =  (10*r~4-12*r~2+3)*r*cos(t) ; 

Zpoly(i, j ,14)  =  (10*r~4-12*r~2+3)*r*sin(t) ; 

Zpoly(i,  j  ,15)  =  20*r~6-30*r~4+12*r~2-l; 

Zpoly (i , j , 16)  =  r~4*cos(4*t) ; 

Zpoly(i, j ,17)  =  r~4*sin(4*t) ; 

Zpoly (i, j ,18)  =  (5*r~2-4)*r~3*cos(3*t) ; 

Zpoly (i, j , 19)  =  (5*r~2-4)*r~3*sin(3*t) ; 

Zpoly(i, j ,20)  =  (15*r~4-20*r~2+6)*r~2*cos(2*t) ; 

Zpoly(i, j ,21)  =  (15*r~4-20*r~2+6)*r~2*sin(2*t) ; 

Zpoly(i, j ,22)  =  (35*r~6-60*r~4+30*r~2-4) *r*cos (t) ; 

Zpoly(i, j ,23)  =  (35*r~6-60*r~4+30*r~2-4) *r*sin(t) ; 

Zpoly(i, j ,24)  =  70*r~8-140*r~6+90*r~4-20*r~2+l; 

Zpoly(i, j ,25)  =  r~5*cos(5*t) ; 

Zpoly(i ,  j  ,26)  =  r~5*sin(5*t) ; 

Zpoly(i, j ,27)  =  (6*r~2-5)*r~4*cos(4*t) ; 

Zpoly(i, j ,28)  =  (6*r~2-5)*r~4*sin(4*t) ; 

Zpoly(i,  j  ,29)  =  (21*r~4-30*r~2+10) *r~3*cos (3*t)  ; 

Zpoly(i,  j  ,30)  =  (21*r~4-30*r~2+10) *r‘'3*sin(3*t)  ; 

Zpoly(i, j ,31)  =  (56*r~6-105*r~4+60*r~2-10)*r~2*cos(2*t) ; 

Zpoly(i, j ,32)  =  (56*r~6-105*r~4+60*r~2-10) *r~2*sin(2*t) ; 

Zpoly(i,  j  ,33)  =  (126*r~8-280*r~6+210*r'‘4-60*r~2+5) *r*cos(t)  ; 

Zpoly(i, j ,34)  =  (126*r~8-280*r~6+210*r~4-60*r~2+5)*r*sin(t) ; 

Zpoly(i, j ,35)  =  252*r~10-630*r~8+560*r~6-210*r~4+30*r~2-l ; 

Zpoly(i , j ,36)  =  r~6*cos(6*t) ; 

Zpoly (i ,j ,37)  =  r~6*sin(6*t) ; 

Zpoly(i, j ,38)  =  (7*r~2-6) *r~5*cos (5*t) ; 

Zpoly(i, j ,39)  =  (7*r~2-6)*r~5*sin(5*t) ; 

Zpoly (i, j ,40)  =  924*r~12-2772*r~10+3150*r~8-1680*r~6+420*r~4-42*r~2+l; 
Zpoly(i, j ,41)  =  r~7*cos(7*t) ; 

Zpoly(i,  j  ,42)  =  r~7*sin(7*t)  ; 
end  '/»i 
end  °/.j 
else 
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"/  R  =  2  "/test  line 

"/  P  =  8  "/test  line 

"/  [X1X2,EL]  =  FEM_Membrane_cylindrical(R,P) ;  "/test  line 
X1X2  =  R; 

for  i  =  1: length (X1X2(: ,1)) ; 
r  =  XlX2(i,2) ; 

t  =  XlX2(i,3); 

x(i , 1)  =  r*cos (t) ; 
y(i , 1)  =  r*sin(t) ; 


Zpoly (i , 

,1,D 

= 

r*cos(t);  '/, (X  Tilt) 

Zpoly (i , 

,1,2) 

= 

r*sin(t) ;  i  (Y  Tilt) 

Zpoly (i , 

,1,3) 

= 

2*r~2-l; 

Zpoly (i , 

,1,4) 

= 

r~2*cos (2*t)  ;  °/o(0  Astigmatism) 

Zpoly (i . 

,1,5) 

= 

r~2*sin(2*t) ;  %(45  Astigmatism) 

Zpoly (i . 

,1,6) 

= 

(3*r~2-2)*r*cos(t)  ;  7„ (X  Coma) 

Zpoly (i , 

,1,7) 

= 

(3*r~2-2)*r*sin(t)  ;  7.(Y  Coma) 

Zpoly (i , 

,1,8) 

= 

6*r~4-6*r~2+l ;  '/.(Spherical) 

Zpoly (i . 

,1,9) 

= 

r~3*cos (3*t) ; 

Zpoly (i , 

,1,10) 

= 

r~3*sin(3*t) ; 

Zpoly (i . 

,1,11) 

= 

(4*r‘2-3)*r‘2*cos(2*t) ; 

Zpoly (i . 

,1,12) 

= 

(4*r~2-3)*r‘2*sin(2*t) ; 

Zpoly (i , 

,1,13) 

= 

(10*r‘4-12*r‘2+3)*r*cos(t) ; 

Zpoly (i , 

,1,14) 

= 

(10*r‘4-12*r‘2+3)*r*sin(t) ; 

Zpoly (i , 

,1,15) 

= 

20*r~6-30*r‘4+12*r‘2-l ; 

Zpoly (i , 

,1,16) 

= 

r“4*cos(4*t) ; 

Zpoly (i , 

,1,17) 

= 

r~4*sin(4*t) ; 

Zpoly (i . 

,1,18) 

= 

(5*r‘2-4)*r‘3*cos(3*t) ; 

Zpoly (i , 

,1,19) 

= 

(5*r‘2-4)*r‘3*sin(3*t) ; 

Zpoly (i . 

,1,20) 

= 

(15*r‘4-20*r‘2+6)*r“2*cos(2*t) ; 

Zpoly (i . 

,1,21) 

= 

(15*r‘4-20*r‘2+6)*r“2*sin(2*t) ; 

Zpoly (i . 

,1,22) 

= 

(35*r~6-60*r~4+30*r~2-4) *r*cos (t) ; 

Zpoly (i . 

,1,23) 

= 

(35*r‘6-60*r‘4+30*r“2-4) *r*sin(t) ; 

Zpoly (i , 

,1,24) 

= 

70*r‘8-140*r‘6+90*r“4-20*r“2+l ; 

Zpoly (i , 

,1,25) 

= 

r“5*cos(B*t) ; 

Zpoly (i , 

,1,26) 

= 

r~5*sin(5*t) ; 

Zpoly (i , 

,1,27) 

= 

(6*r‘2-5)*r‘4*cos(4*t) ; 

Zpoly (i . 

,1,28) 

= 

(6*r‘2-5)*r‘4*sin(4*t) ; 

Zpoly (i . 

,1,29) 

= 

(21*r‘4-30*r‘2+10)*r‘3*cos(3*t) ; 

Zpoly (i . 

,1,30) 

= 

(21*r‘4-30*r‘2+10)*r‘3*sin(3*t) ; 

Zpoly (i . 

,1,31) 

= 

(56*r‘6-106*r~4+60*r‘2-10)*r‘2*cos(2*t) ; 

Zpoly (i . 

,1,32) 

= 

(56*r‘6-105*r*4+60*r‘2-10)*r‘2*sin(2*t) ; 

Zpoly (i . 

,1,33) 

= 

(126*r~8-280*r‘6+210*r~4-60*r~2+5)*r*cos(t) ; 

Zpoly (i , 

,1,34) 

= 

(126*r*8-280*r‘6+210*r*4-60*r‘2+5)*r»sin(t) ; 

Zpoly (i , 

,1,35) 

= 

282*r‘  10-630*r‘8+560tr~6-210*r‘4+30*r‘2-l ; 

Zpoly (i , 

,1,36) 

= 

r‘6*cos (6*t) ; 

Zpoly (i . 

,1,37) 

= 

r~6*sin(6*t) ; 

Zpoly (i , 

,1,38) 

= 

(7*r‘2-6)*r‘5*cos(5*t) ; 

Zpoly (i , 

,1,39) 

= 

(7*r‘2-6)*r‘5*sin(5*t) ; 

Zpoly (i . 

,1,40) 

= 

924*r‘12-2772*r~10+3160*r‘8-1680*r‘6+420*r‘4-42*r"2+l; 

Zpoly (i , 

,1,41) 

= 

r“7*cos(7*t) ; 

Zpoly (i . 

,1,42) 

= 

r"7*sin(7*t) ; 

end;  "/  : 

L  loop 
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end;  °/0nargin  if  then 
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Appendix  C.  Other  Interesting  Notes 


C.l  Amplifier  Characterization  Data 


Table  15:  Test  Parameters  for  ACX  Amplifier  Characterization 


Parameter 

Setting 

Test  Range 

0  to  5  kHz 

Input 

Gaussian  White  Noise 

Amplitude 

10  mV 

Window 

Hanning 

Lines 

3200 

Averages 

10 

lifciM  SipnalCalc  620  ^^ansfe^unetion  rC:\Documcnts  and  Settingsmser\MvDocuments\Userfiles\AmDT^sts\BopenAmD.trf.  Currcnt:RunOOQ37.  Next:S...  I—  lPt?Hxl 


Figure  68:  ACX  Amplifier  Characterization  from  0  to  5  kHz  (S/N:  8H263) 
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Figure  69:  ACX  Amplifier  Characterization  from  0  to  500  Hz  (S/N:  8H263) 
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Table  16:  Test  Parameters  for  TREK  Amplifier  Characterization 


Parameter 

Setting 

Test  Range 

0  to  100  Hz 

Input 

Gaussian  White  Noise 

Amplitude 

10  mV 

Window 

Hanning 

Lines 

3200 

Averages 

10 
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Figure  70:  TREK  Amplifier  1  Characterization  from  0  to  100  Hz  (S/N:  118A) 
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Figure  71:  TREK  Amplifier  2  Characterization  from  0  to  100  Hz  (S/N:  118B) 


Figure  72:  TREK  Amplifier  3  Characterization  from  0  to  100  Hz  (S/N:  121A) 
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Figure  73:  TREK  Amplifier  4  Characterization  from  0  to  100  Hz  (S/N:  121B) 
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Figure  74:  TREK  Amplifier  5  Characterization  from  0  to  100  Hz  (S/N:  120A) 
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Figure  75:  TREK  Amplifier  6  Characterization  from  0  to  100  Hz  (S/N:  120B) 


Figure  76:  TREK  Amplifier  7  Characterization  from  0  to  100  Hz  (S/N:  115A) 
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Figure  77:  TREK  Amplifier  8  Characterization  from  0  to  100  Hz  (S/N:  115B) 


C.2  Relevant  Telephone  Numbers 


Table  17:  Phone  Numbers 


Name 


Phone  Number 


Adaptive  Optics  Assoc.  (AOA) 

Herb  DeSilva  (AOA  software) 

Jay  Anderson  (AFIT/ENY  Lab  Supervisor) 
Dr.  David  Mollenhauer  (ML  contact) 


Brian  Langley  (ML  coop  student) 
Joel  Johnson  (ML  contact) 


617-806-1400 
617-806-1816 
937-255-3636  x4865 
937-255-9727  (office) 
937-255-9059  (lab) 
937-255-9070  (branch  office) 
937-255-9059  (lab) 
937-255-4651 
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C.  3  Zernike  Polynomials 

Given  below  are  the  first  42  Zernike  Polynomials.  Each  term  has  a  coefficient  multiplier  which 
scales  the  set.  Upon  combination,  the  surface  deflection  is  given  in  terms  of  radius  (r)  and  theta  (0). 
These  coefficients  are  scaled  so  that  the  surface  deflection  is  in  terms  of  microns  (yim). 

Zx  =  rcos(6) 

Zi  =  rsin(9) 

Z3  =  2r2  -  1 
Z4  =  r2cos{29 ) 

Z5  =  r2sin{29) 

Zq  =  (3r2  —  2  )rcos{9) 

Z7  =  (3  r2  —  2  )rsin{9) 

Z&  =  6r4  -  6r2  +  1 
Zg  =  r3cos{39) 

Zio  =  r3sin(39) 

Zn  =  (4  r2  —  3)r2cos{29) 

Z12  =  (4r2  —  3)r2sin(29) 

ZXz  =  (10r4  —  12r2  +  3  )rcos(9) 

Z14  =  (10r4  —  12r2  +  3  )rsin{9) 

Z15  =  20 r6  -  30r4  +  12r2  -  1 
Ziq  =  r4cos(4:9) 

Zn  =  r4sin(49) 

Zia  =  (5r2  —  4)r3cos(30) 

Zw  =  (5  r2  —  4)r3sm(30) 

Z20  =  (15r4  -  20r2  +  6)r2cos(29) 
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Z21  =  (15 r4  -  20 r2  +  6)r2sin(29) 

Z22  =  (35r6  —  60r4  +  30r2  —  4  )rcos{9) 

Z23  =  (35  r6  —  60r4  +  30?’2  —  4  )rsin{9) 

Z24  =  70r8  -  140r6  +  90r4  -  20r2  +  1 
Z23  =  r5  cos  (59) 

Z26  =  r5sin(59) 

Z27  =  (6r2  —  5)r4cos(40) 

Z28  =  (6r2  —  5)r4sm(40) 

Z29  =  (21r4  -  30r2  +  10)r3cos(3(9) 

Z30  =  (21  r4  -  30r2  +  10)r3sm(36>) 

Z31  =  (56r6  -  105r4  +  60r2  -  VS)r2cos(29 ) 

Z32  =  (56r6  -  105r4  +  60r2  -  lO)r2sm(20) 

Z33  =  (126r8  -  280r6  +  210r4  -  60r2  +  5 )rcos(9) 

Z34  =  (126r8  -  280r6  +  210r4  -  60r2  +  5 )rsin(9) 

Z35  =  252r10  -  630r8  +  560r6  -  210r4  +  30r2  -  1 
Z36  =  r6cos(69) 

Z37  =  resin(69) 

Z38  =  (7r2  —  6)r5cos(59) 

Z39  =  (7  r2  —  6)r5sin(59) 

Z40  =  924r42  -  2772r40  +  3150r8  -  1680r6  +  420r4  -  42r2  +  1 
Z44  =  r7cos(79) 

Z42  =  r7  sin(79) 
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